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Abstract 

The superstring theory is now regarded as the most promising candidate for the unification of the 
Standard Model and gravity, and this field has been rigorously investigated. However, we have seen 
a setback of the perturbative analysis of the superstring theory, because it has so many vacua that 
we have no way to determine which are the true ones. In order to remedy this problem and thus 
for the superstring theory to have a power to predict our real four-dimensional world, we need the 
nonperturbative formulation of the superstring theory. 

In the late 1990's, our understanding of the nonperturbative aspects of the superstring theory 
have been greatly deepened. Especially, the large- A'" (A'^ is the size of the matrices) reduced models 
have been proposed as the nonperturbative formulation of the superstring theory. One of the most 
promising candidates is the IIB matrix model, which is defined by the dimensional reduction of 
the ten-dimensional A/" = 1 super- Yang-Mills theory to zero dimension. It has been conjectured 
that this model resurrects the behavior of the string theory at the large-A limit. There have been 
a lot of interesting discoveries of the IIB matrix model, such as the dynamical generation of the 
four-dimensional spacetime and the interpretation of the diffeomorphism invariance. 

On the other hand, there are a lot of problems to surmount, if a large-A reduced model is to be 
an eligible framework to unify the gravitational interaction. Firstly, it is still an enigma how we can 
realize the local Lorentz invariant matrix model. In addition, we need to understand how we can 
describe the curved spacetime more manifestly, in terms of a large-A reduced model. 

This thesis discusses several attempts to address these issues concerning the gravitational inter- 
action. This thesis is based on the following works [26,38,46,60]. 
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1 Introduction 



One of the main themes in the elementary particle physics is unification. There are four kinds of inter- 
action in the nature; the weak interaction, the strong interaction, the electromagnetic interaction and 
the gravitational interaction. In 1967, Glashow, Weinberg and Salam succeeded in the unification of the 
weak and electromagnetic interaction in terms of the SU{2) x U{1) gauge group. The ensuing success is 
the advent of the "Standard Model", described by the SU{3) x SU{2) x U{1) gauge group. The Stan- 
dard Model is striking in the sense that it describes all the experimental phenomena by setting the 18 
parameters in the model. 

However, there are two serious drawbacks to the Standard Model. First is that we cannot explain 
theoretically how these 18 parameters are fixed. Namely, we have to rely on the experimental data in 
setting these parameters. Second is that the Standard Model does not unify the gravitational interaction. 

Now, the superstring theory is regarded as the most promising candidate to resolve these two draw- 
backs. The superstring theory regards not a zero- dimensional point but a one-dimensional string with 
the length Ip = 10^'^'''cm as the fundamental object. Here, Ip is called "the Planck Scale", which is the 
fundamental scale in the superstring theory. The superstring theory incorporates not only the gauge 
particles but also the gravitons in its oscillation modes. There is an infinite tower of the mass level in 
the oscillation modes. The massless modes of the open string, which has two ends, arc regarded as the 
gauge particle. On the other hand, the massless modes of the closed string, which constitutes a closed 
circle like a rubber band, gives the spin-2 graviton. In this sense, the superstring theory is thought to 
include not only the Standard Model but also the gravitational theory in the low-energy limit, where 
the massive modes of the superstring theory are ignored. The superstring theory is fascinating in that it 
incorporates no free parameters. It is in contrast to the Standard Model, which suffers the problem of a 
plethora of free parameters. The superstring theory has a possibility to explain the gauge group of the 
Standard Model, the number of the generations of the quarks, the mass of the Higgs particle, quarks or 
leptons, the dimensionality of our spacetime • • •, from a theory without any parameter. 

We have seen the so-called "first string boom" in the early 1980's, in which the pcrturbative aspects 
of the superstring theory have been elucidated. The striking discovery in 1984 is that the superstring 
theory is free from the divergence of the gravitational energy. Physicists have been so far puzzled by the 
UV divergence in the quantization of the gravitational interaction. This conundrum has been resolved by 
the superstring theory which alleviates this divergence. Moreover, it has also been known that there are 
only five kinds of superstring theory well-defined in the ten dimensions. They are the type I, type IIA, 
type IIB, SO{32) heterotic and Es x Es heterotic superstring theory. Especially, the Es x Es heterotic 
superstring theory is propitious for the immersion of the Standard Model in the superstring theory. The 
exceptional Lie algebra Eg includes the Lie algebra SU{3) x SU{2) x U{1), which is the gauge symmetry 
of the Standard Model, as its subalgebra. 

The superstring theory is defined in the ten-dimensional spacetime, but this is not so problematic 
because it is a theory of gravity. The spacetime is regarded as being given not a priori but dynamically 
from the classical solution of the theory. It is now believed that the superfiuous six dimensions are 
compactified at an extremely small size at the incunabula of the universe. Mathematically, the six 
superfluous dimensions are compactified by the Calabi-Yau manifold. 

However, the superstring theory suffers a serious setback. Since the superstring theory is defined 
only perturbatively, there is no telling which is the true vacuum. While there are only five well-defined 
superstring theories, there are a myriad of ways to compactify the ten-dimensional superstring theory into 
the four dimensions and thus we are overwhelmed by the infinite vacua. This means that the siiperstring 
theory does not have the ability to predict our four-dimensional world, and we cannot tell whether the 
superstring theory dynamically generates the Standard Model. In order to surmount this difficulty, we 
definitely need the constructive definition (i.e. the definition without the perturbation). 

The so-called "first string boom" was over at the end of the 1980's. However, there were important 
discoveries at this period about the noncritical string theory and the bosonic matrix model (the extensive 
review can be foimd in [1,2]). Distler and Kawai [72] succeeded in the quantization of the non-critical 
string theory up to one dimension via the conformal gauge. In addition, there is a striking relation 
between the bosonic string theory and the matrix model, proposed by F. David [70]. For the simplest 
zero-dimensional non-critical string, the corresponding matrix model is a one-matrix model described by a 
simple NxN hermitian matrix. The correspondence has been found for the central charge c = 1 — 
(for m = 2, 3, 4, • • •), in which the string corresponds to the multi-matrix model. The path integral of the 
worldsheet for all genera has been approximated by the random triangulation of the worldsheet. This 
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has led us to identify the Feynman rule with that of the one-matrix model. In this way, we have seen 
the formal correspondence between the matrix model and the string theory. Even more crucial is that 
the matrix model can be solved nonperturbatively. Brezin and Kazakov [74] solved the matrix model by 
means of the orthogonal polynomial method. They extracted the nonperturbative aspects of the bosonic 
matrix model via the differential equation called Painleve equation. They calculated the parameter called 
the string susceptibility, and found that this agrees with the calculation of Distler and Kawai [72] for the 
noncritical string. 

We have failed in extending this idea to the " super" string, and this progress is limited to the bosonic 
string. This is due to the same difficulty as we face in putting the chiral fermion onto the lattice. Thus, 
these discoveries, per se, do not give a technical clue to the constructive definition of the superstring 
theory. The 'state-of-the-art' matrix models (such as the IIB matrix model) do not inherit the same 
techniques as the old matrix model. Nevertheless, these discoveries serve as an important touchstone for 
the legitimacy of the large- iV reduced models as the constructive definition of the string theory. 

In the late 1990's we have seen the so-called "second string boom", in which we grasped the nonper- 
turbative aspects of the superstring theory. During the second string boom, Polchinski discovered the 
soliton-like object called the D-brane. The discovery of the D-brane led to the idea of the two kinds 
of the duality of the string theory. First is the T-duality, which is the duality between the large and 
small radius of the compactification. Second is the S-duality, which relates the strong coupling and the 
weak coupling. It is now believed that the five ten-dimensional superstring theories are related with one 
another through the T/S duality. This is an illuminating discovery, in the sense that it is reminiscent of 
the more fundamental theory whose certain limit might reproduce these five superstring theories. 

Another big progress is the proposal for the constructive definition of the superstring theory via the 
large-iV reduced model. Here, the 'large-A^ reduced model' means the model defined by the N x N 
hermitian matrices. In 1996, Banks, Fischler, Shenker and Susskind proposed the matrix model defined 
by the dimensional reduction of the ten-dimensional super Yang-Mills theory to one dimension — the 
so-called BFSS model [4] . This matrix model is more related to the type IIA superstring theory, and the 
effect of the type IIA supergravity is induced by the one- loop effect. 

Other attempts are the IIB matrix model (the IKKT model) [5] and the matrix string theory [6], 
which are defined by the dimensional reduction of the ten-dimensional super Yang-Mills theory to the 
zero and two dimensions, respectively. Especially, the IIB matrix model is now regarded as one of the most 
promising candidates for the constructive definition of the superstring theory (this model is extensively 
reviewed in [16,75]). As its nomenclature shows, this model is deeply related to the type IIB superstring 
theory. The simplest and the most direct correspondence is that the IIB matrix model is defined by the 
matrix regularization of the Grcen-Schwarz action of the type IIB superstring theory; namely by replacing 
the Poisson bracket with the commutator of the N x N matrix. In addition, it has been shown [7] that 
the Wilson loops satisfy the string field equations of motion for the type IIB superstring in the light-cone 
gauge. The prominent feature of this model is that it does not include any free parameter. The overall 
coefficient which is an artifact of the coupling constant before being reduced is absorbed by a simple field 
redefinition. 

There liavc^ been hitherto a lot of progressc^s with respect to the IIB matrix model. First is the relation 
to the gravitational interaction. The most fundamental feature is the spacetime J\f = 2 supersymmctry [5] 
when we regard the eigenvalues of the bosonic matrices as the spacetime coordinate. This urges us to 
interpret the eigenvalues as the spacetime coordinate. In addition, the graviton and dilaton exchange 
has been calculated in [5]. The interpretation of the general coordinate invariance is proposed in [14], in 
which the permutation invariance of the eigenvalues is identified with the general coordinate invariance 
of the low-energy effective action. 

Second is the relation to the noncommutative geometry, which has been widely investigated since 
Seiberg and Witten elucidated the relationship with the superstring theory [18]. By expanding the IIB 
matrix model around the classical background, it reproduces the noncommutative Yang-Mills theory [17]. 

Third is the dynamical generation of the four-dimensional spacetime. In [11], it has been found that 
the Hausdorff dimension of the eigenvalue distribution is four. In [36], Nishimura and Sugino proposed the 
breakdown of the Lorentz symmetry from 50(10) to 50(4) by means of the third order of the Gaussian 
expansion. Their analysis has been extended to the seventh order in [39,48]. 

The discovery of the IIB matrix model implies the trilateral relationship of the three important notions 
in the elementary particle physics; the superstring theory, the quantum field theory and the matrix model. 
The discovery of the relationship between the one-matrix model and the bosonic string theory is another 
good example to relate these notions. In the course of the "second string boom" we have other important 
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discoveries for this trilateral relationships. In 1997, Maldacena advocated the duality between the type 
IIB superstring theory on the AdSz x spacetime and the AA = 4, four-dimensional super Yang-Mills 
theory. This is the so-called "AdS/CFT correspondence" [9,12,15]. This correspondence has been 
conjectured through the coincident symmetries of these two theories. Since the proposal of Maldacena, 
the AdS/CFT correspondence has been rigorously investigated by many authors (among them is the 
author's work [32], while we do not delve into this work in this thesis). Another important discovery is 
the so-called Dijkgraaf-Vafa duality [41,43,45]. This is a duality between the J\f = 1 four-dimensional 
super Yang-Mills theory and the one-matrix model, which used to attract much attention in relation 
to the bosonic string theory in the early 1990's. This discovery reinstated the one-matrix model into 
many authors' attention. At first [41,43,45], this relationship is understood with the intervention of 
the topological string theory. However, Cachazo, Douglas, Seiberg and Witten [47,49,53] explained this 
duality directly through the super Yang-Mills theory without topological string theory. They elucidated 
that these two theories comply with the same Schwinger-Dyson equation, employing the notion of the 
chiral ring and the Konishi anomaly [67,68]. Moreover, in [54,58], the direct relationship between the 
large- A'' reduction and the Dijkgraaf-Vafa duality has been found. In this way, in the course of the "second 
string boom", many interesting discoveries, including the IIB matrix model, have been unraveled. In the 
future, we will surely have a synergy in the studies of these three germane areas. 

While the IIB matrix model has a lot of interesting properties, there is a plenty of room to investigate. 
The more manifest correspondence with the gravitational interaction is definitely a $64,000 question. If 
the IIB matrix model or its extensions are to be the eligible frameworks to unify the gravitational 
interaction, we must take seriously the relation to the gravitational interaction. While several important 
evidences have been found, there are a lot of questions to clarify. It has been an enigma how we can 
describe the local Lorentz symmetry in terms of the large- A?^ reduced models. In addition, the definition of 
the IIB matrix model relies heavily on the flat background. Especially, the IIB matrix model incorporates 
the background of the flat space, not the curved space. This prohibits the perturbation around the curved 
space background. 

In this thesis, we discuss some of the attempts [26, 38, 46, 60] to address these issues. This thesis is 

organized as follows. In Section 2, we give a brief review of the IIB matrix model. Especially, we focus on 
the development of the IIB matrix model in relation to the gravitational interaction. Moreover, we have 
a careful look at the alterations of the IIB matrix model that incorporate a curved-space background, in 
relation to some of the author's works [46,60]. 

In Section 3, we review the author's works about the supermatrix model [26,46], based on the 
osp{l\32, R) super Lie algebra. The studies of the supermatrix models give a rich perspective for the 
generalization of the IIB matrix model in relation to the gravitational interaction. We start with in- 
vestigating the correspondence of the nongauged osp(l|32, R) model with the IIB matrix model, paying 
attention to the supersymmetry. We next consider the so-called "gauged version", namely the model 
whose Lorentz symmetry and the gauge symmetry are mixed together. This is an important idea in real- 
izing the matrix model equipped with the local Lorentz symmetry, because the eigenvalues of the bosonic 
matrices arc identified with the spacetime in the IIB matrix model. We next go back to the nongauged 
version of the osp(l|32,i?) supermatrix model including the mass term [46]. It has been known that the 
IIB matrix model with the tachyonic mass term, whose review we give in Section 2, incorporates the 
curved-space background because the trivial commutative background is destabilized. We focus on the 
similarity of such models and the massive supermatrix model. 

In Section 4, we elaborate on the idea of the "gauged matrix model" without the supermatrix model 
[38]. Especially, we allocate the odd-rank matrices of the ten dimensions for the matter field, and the 
even-rank matrices to the parameters of the local Lorentz transformation. When we mix the Lorentz 
symmetry and the gauge symmetry, the model must inevitably include the higher-rank fields in order for 
the action to be invariant and for the algebra of the symmetry to close with respect to the commutator. 
This is a good news in that the rank-3 fields are identified with the spin connection in the curved space. 
We consider the explicit identification of the matrices with the differential operator, and clarify that 
the bosonic part actually reduces to the Einstein gravity in the low-energy limit. We also discuss the 
structure of the Af = 2 supersymmetry. 

In Section 5, we study the stability of the fuzzy sphere background in the quantum sense. While it is 
easy to discuss the stability of the curved-space background (all we have to do is just to plug the solution 
into the action), it is by no means easy to discuss the stability in the quantum sense. In this section, 
we address this question via the heat bath algorithm of the Monte Carlo simulation. We focus on the 
simplest case; namely the bosonic three-dimensional IIB matrix model with the Chern-Simons term. We 
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find that this bosonic matrix model has a first-order phase transition between the Yang- Mills phase (in 
which the quantum effect is large and the fuzzy sphere classical solution is unstable) and the fuzzy sphere 
phase (in which the model is subject to the meager quantum effect and the fuzzy sphere is stable), as we 
change the radius of the fuzzy-sphere solution (equivalently, the coefficient of the Chern-Simons term). 
Moreover, we find that this model has a one-loop exactness in the fuzzy sphere phase in the large-TV 
limit. The latter result is especially exciting, in the sense that this helps the analysis of the dynamical 
generation of the gauge group. 

Section 6 is devoted to the conclusion and the future outlook of the work. 

In Appendix A, we summarize the notation of this thesis. We give the detailed notation of the gamma 
matrices and the supermatrices. In addition, we give a proof of the miscellaneous formulae of the gamma 
matrices. 

In Appendix B, we give the calculation of the Seeley-de-Witt coefficients, which plays an important 
role in Section 4. 

In Appendix C, we give the detailed recipe for the heat bath algorithm of the matrix model. We start 
with the quadratic one-matrix model in full detail, because this gives the key idea of the simulation of the 
IIB matrix model with/without the Chern-Simons term. We nextly discuss the simulation of the quartic 
one-matrix model. Finally, we discuss the heat bath algorithm for the IIB matrix model with/ without 
the Chern-Simons term. 



2 Brief review of the IIB matrix model 

In this section, we give a brief review of the IIB matrix model [5] . As we have mentioned in the introduc- 
tion, the IIB matrix model is regarded as the most promising candidate for the constructive definition of 
the superstring theory. This model is also called " the IKKT model" , where IKKT is the acronym for the 
authors of [5]; Ishibashi, Kawai, Kitazawa and Tsuchiya. 

2.1 Definition and symmetry of the IIB matrix model 

The IIB matrix model is defined by the dimensional reduction of the ten-dimensional M — \ super 
Yang-Mills theory to zero dimension as 

S^--Tr{- Y,[A,,A,][A^^,A'^] + -Y,4'T^[A,m]- (2-1) 

^ \ ti,u=0 fj.=0 I 

Here, A^ are the ten-dimensional bosonic vector, and ?A is the ten-dimensional Majorana-Weyl (hence 
sixteen-component) fermion. Both and are promoted to the N x N hermitian matrices. This model 
incorporates the 5*0(9, 1) Lorentz symmetry and the SU{N) gauge symmetry. This model is a totally 
reduced model in that it is reduced to zero dimension. The prominent feature of this model is that it 
has no free parameter. The overall coefficient g can be trivially absorbed into the fields by the following 
redefinition. 

A^^giA^, ^/j^gi^. (2.2) 

The coefficient g is an artifact of the coupling constant of the super Yang-Mills theory before being 
reduced, and is nothing but a scaling parameter. 

The path integral of the IIB matrix model is given by 

Z = J dAdije-^^''^ (2.3) 

where »S'(^) is defined in the ten-dimensional Euclidean space by the Wick rotation of Aq and in the 
action 1)2. The convergence of the path integral of the IIB matrix model is not trivial, because the 
gauge group SU (N) is not compact. However, the convergence of the path integral is discussed in the 
paper [24,28] (the review is found in the Ph.D. thesis of Austing [34]). The authors corroborated that 
the path integral converges for^ d — 4,6, 10. Only for the bosonic part, the path integral is shown to 
converge for d > dc- Here, dc = 5 for the SU{2) gauge group, dc — A for the S'C/(3) and dc — "i for all the 
other simple Lie group. 

is the dimension of the IIB matrix model, and the model Il2.ll is defined for d = 10. More explicitly, the d-dimcnsional 
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2.1.1 Relation to the type IIB superstring theory 

The IIB matrix model has a deep relation to the type IIB superstring. Here, we discern that the IIB matrix 
model is obtained by the matrix regularization of the Green-Schwarz action of the type IIB superstring 
theory. The Green-Schwarz formalism is defined so that the spacetime supersymmetry should be more 
manifest. The bosonic Nambu-Goto action is given by 

Sb = J (faV^h'^^d^X^'dpX^. (2.4) 

Here, hap is a metric on the worldsheet, and has a Lorentz signature ry"'^ = diag(— , +). A naive guess 
for a supersymmetric extension is to replace daX^ with 

H^ = daXf" - iie^vdaO^ - Pvdc^e^). (2.5) 

H(^ is trivially invariant under the supersymmetry transformation 

SsO'^^e'^, SsX" ^iie^Ve^ -e^r^e"^). (2.6) 

Here, A runs over 1, 2, and this is an index for the supersymmetry. However, the naive replacement of 
daXf^ with for the action (|2.4() is not enough, because it does not have a k symmetry and hence it 
has twice as many degrees of freedom as it should. Instead, we give the following action: 

Sgs = ^ J d^a (V^/i"^H^H;3m + 2ie"f'daX'^{FT^dp0^ + PT^df^e^) + 2e^^ {e^r'^da9^){S^T ^8^9^ 

In this section, we define the rank-2 epsilon tensor as — 1 (and hence eoi — —1). This action is 
defined for the d = 10 dimensional spacetime, and the index fi runs over 0, 1, 2, • • • , 9. The supersymmetry 
transformation of this action is 

SsSgs =2^1 d'ae''^[ie'T^do.9'WT,dp0') - ie^T^da.9^)iS^T^d09% (2.7) 

where we drop the surface terms. The term (|2.7f) is shown to vanish (i.e. the action Sqs is supersymmetry 
invariant) by noting the following rewriting: 

A = e'^'^lT^'^a99T^^(s9 ^eT''96T^9' -eT''e'9T^6. (2.8) 

Here, — t and ai ~ a. 9 — ^ and 9' — A can be further written as 

A = '^[eT''99Tf,9' + 1^9' 9T ,,9 + eT^'99'Tf,9] + l--^[er''99Tf,9'] - l--^[eTf'99Tf,9]. (2.9) 
3 3 OT 3 oa 

The first term vanishes due to the Fierz identity 

er^7/.[i7/i2r^V^3] = 0. (2.10) 

We delegate the proof of this Fierz identity to Appendix IA.1.61 The second and third terms vanish 
because they are nothing but surface terms. 

This model also incorporates the symmetry named the k symmetry: 

SJ"^ = 2r^H;jK^", S^X^' = i9^T^'5J^ - iPT''S^9'', 

S^iV^h"/^) = -16tV^{d^9^K^I^Py - d^pK^f^P^^). (2.11) 



supersymmetric IIB matrix model means the following action 

Cd-l d-l \ 

^1,1^=0 M=0 / 

which is well-defined for d = 3, 4, 6, 10. The explicit definition of the d-dimensional bosonic IIB matrix model is given by 

d-l 
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Here, is a projection operator defined by 

2' " ~ 



= 7;{h°''^ ±-T=)- (2-12) 



This satisfies the fohowing properties 

pfhp^P^^=Pf, Pfhp^Pl^ = Q, pfp^s^pfpf^ (2.13) 
which follow from the inversion rule 

h^P^U \" -J""' ) . (2.14) 



h \ —hio hoa 

Here, the parameters of the K-symmetry is subject to the following projection rule: 

P^^kI = k''', Pfnl^K^''. (2.15) 
The K-symmetry transformation varies the action as 

6.Sgs - ^ ydVe"'3([(9^0"i)r'^(a„0i)][^ir^5«0i] + 2[(a„0"i)r'^(5«0i][^ir^a^0i] 

-[idpP)r''{d^e')][PT,sj'] - 2[{dJ^)r^sjWT^dp9^]) , (2.16) 

up to the surface terms that vanish in the integral. We can show that the Green-Schwarz action is 
invariant under the K-symmetry with the help of the identity H2.10|l . 

Nextly, we explain how the matrix regularization of the Green-Schwarz action of the type HB super- 
string theory leads to the IIB matrix model. For the type IIB superstring theory, the chirality of the two 
spinors 9^ and 6^ is identical. Therefore, we set 9^ = 9^ = 9. This simplifies the term as H^ — daX^. 
Then, the action is simplified as 

Sgs = ^ J (h"f^V^^c.X^'^fiX^ + 4ie"'3a„x^^r^5/30) . (2.17) 

Integrating out haf}, we obtain 

Sgs = ^ J d^a (^^- dct{^^X^^^pX^) + 2ie''''daX^^9r^dp9^ 

-I /■ .o / n 



Jd^a ( ^J-^ie'-'^da^Xt^d^X^)^ + 2te'^f'do.Xf^9T^dp9] . (2.18) 



27ra' 

In the last equality, we note that this determinant can be rewritten as 



detid^x'^dpx^) - idoX''d„x^){diX''diX,)-{doX''diX^){diX''doX,) 

= ^{e^f'd^X'^dpX^)^ (2.19) 

We introduce the Schild action and discern that this reduces to the Green-Schwarz action. The Schild 
action is defined by 

^Schild ^-J"^'" (v^« (^i^^^'^'^'^^' + ^^r'^{X^,^}) + V^sbN^ . (2.20) 

Here, we introduce the scalar density \/—hs as an independent variable. The Poisson bracket {q,p} is 
defined by 

{q,p} = ^^e"f'daqdpp. (2.21) 
Varying this Schild action with respect to \/—hs, we obtain 

V^s = 2^ ^^-{e-0d^Xf^d^X,)^. (2.22) 
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This scalar density is clearly identical to the determinant in the Green-Schwarz action \/ —h up to a 
constant. Plugging this into the action, we obtain 

^Schiid 1^^"" [Vd>^~{e-Pd^x^^df3X^Y + ^e"^vir'^a„x^a^v) , (2.23) 

which reduces to the Green-Schwarz action by setting (a, 6) = {-^, wna' )• discuss the matrix 

regularization of the Schild action (|2.2()|l . We perform the following replacement: 

-ilj^Nll J d'a^^Tr. (2.24) 

By setting a = p- and dropping the term y/—hsb, the Schild action arrives at the IIB matrix model. 

We next discuss the supersymmetry of the Green-Schwarz action. We identify the two combinations 
of the parameter of the original supersymmetry and the k symmetry, to obtain a new supersymmetry 
transformation. To this end, we relate the parameters of the supersymmetry and the k symmetry as 

/t^" = -^r'^(e' - e')n0M^-^ = ^r^(e' - e^)'np^Pf- (2-25) 
We take the following gauge for the Schild action of the type IIB superstring theory: 

9^ = e^{=e), ^^^^o,x^', ^^e^'^n^m = (2.26) 

V-hs 

Thus, the fcrmionic supersymmetry transformation is written as 

= l(}- l^'^iX,^ ^^}) i^' ~ ^'). '5.^^' = ^ ^r^^i^^^. ^^}) - ^')- (2-27) 
Thus, the k, symmetry transformation for X^ is 

= -i^r^(ei - 62). (2.28) 
We introduce the following redefinition of the supersymmetry parameter: 

X = 2(ei-e2)^ y^(!i±f!). (2.29) 
The combination of the supersymmetry and the k symmetry is thus given by 

((5e + = ixr^fl, (<5, + 8^)9 = x' + ^r'^''{X^, X,}x- (2.30) 

2.1.2 Supersymmetry 

Now, let us have a careful look at the supersymmetry of the IIB matrix model. We have formulated the 
supersymmetry of the Schild action of the type IIB superstring. In the IIB matrix model we consider 
that the matrix-regularized version of that supersymmetry is inherited. The M ~ 2 supersymmetry of 
the IIB matrix model is then distinguished by 

• homogeneous: J.^^V = AJ\Vi''' b'f^ A^, = ieVf^ij}. 

The feature of the 'homogeneous supersymmetry' is that this supersymmetry transformation de- 
pends on the matter fields Afj, and "tp. And this supersymmetry transformation vanishes if there is 
no matter field. 

(2) (2) 

• inhomogeneous: ip = ^, A^ = 0. 

The feature of the 'inhomogeneous supersymmetry' is that the translation survives without the 
matter fields. 

The commutators of these supersymmetries give the following important results: 

(1) [jW,4i)]A^ = 0, [5W,4i)]V = 0, (2.31) 

(2) [4f,4f]A^ = 0, [4f,4f]=0, (2.32) 

(3) [5'^^\5f]A^^-ieT^i, [5'^^\6f]^ = 0. (2.33) 
These properties can be verified by taking the difference of the two supersymmetry transformations. 
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1. This is the most comphcated to compute. For the gauge field, we should consider the following 
transformation 

5(1) 5(1) I 

^ A^ + ie2T^iP ^ A^+i{ei+e2)T,,^^ -e2T,,[A,,Ap]r''Pei, (2.34) 
5(1) 5<i) 1 

^ A^ + ieiT^iJ; A^+i{ei+e2)T,,^- -eiT,,[A,,Ap]T''Pe2. (2.35) 
Then, the commutator is 

[6il\sil^]A^ - „ie2r^[^.,A,]r''''6i + ^e^T^[A,,Ap]n2. (2.36) 
Utihzing the formula Tf^VP = + ■qf'TP - ifPr"", we obtain 

[5il^,Sil^]A^ = 2e^rPe2[A„Ap]. (2.37) 
For the fermions, we have only to repeat the similar procedure: 



5(1) 5(1) 

i/) -I- 

2' 



^ 3 ^+-[A^^A,]TP''e2 ^ i> + -[A^,A,]T^'''iei+e2)-[A^,eiT,ij]T''''e2, (2.38) 



^ ^+'-[A^,A,]TP''e,'M i,+ '-[A^,A,]T^^ie,+e2)^[A^,e2T,iP]T'^''e,. (2.39) 

Using the formula of the Fierz transformation, 

-^eiF^esFpF'^V-TTr^,,- 
16 16 X 5! 

(2.40) 



eir,V'r'"'e2 = eiF^es^ - — eiF^esFpF'^V - ——eiTP^-P'e2Tp,...p,T''ij + (rank 3 term). 



^Tf^[A^,i;]=0, (2.42) 



whose proof we present in Appendix IA.1.51 we verify that the commutator of the supersymmetry 
transformation is 

[Sil\Sil^]i; = 2[^,eir''e2Ap]. (2.41) 
Here, we utilize the equation of motion 

d5 _ 1 
d'ip g 

in order to eliminate the second and the third terms of the Fierz identity 1)2.4011 . 

Next, we note that the commutators of the supersymmetry transformation (|2.37|) and (|2.41l) vanish 
up to the gauge transformation. The gauge transformation of the IIB model is to multiply the 
unitary matrix ia G SU{N). The gauge transformation is expressed in the infinitesimal form as 
follows : 

Af„i; ^ Af, + i[Af„a], i/' + i[V',a]. (2.43) 

The supersymmetry transformation H2.37|l and H2.41|l can be gauged away by the gauge parameter 
a — 2eiTPe2Ap. We now complete the proof of H2.31|l up to the gauge transformation. 

(21 

2. This is trivial because the supersymmetry involves only a constant. 

3. This can be proven by taking the difference of these two transformations 

5(1) _ 5(1) _ 5<"' 

A^ A,, ^ Ap + leTp^, whereas A^, A^ + leTpip -> + leTpii/j + C), 

V' + tp + ^+^ri"'[Af„A^]e, whereas 



^ ^P+'-r'^'^[Ap,A,]e% ^ + ^+'-r^^[A^,A,]e. (2.44) 
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This completes the proof of the commutation relation of the supersymmetry transformation. When we 
take the linear combination of the homogeneous and the inhomogeneous supersymmetry as 

=<5(i)+<5(2), (2.45) 

the commutator is written as 

[5("\^fV = 0, [5[''\~5f]A^,^~2ieT,,i5''^ (2.46) 

where a,/3 run over 1,2. The commutation relation (|2.46(l has a serious consequence. When we regard 
the eigenvalues of the bosonic matrices as the spacetime coordinate, the IIB matrix model carries a 
spacetime Af = 2 supersymmetry. Namely, the commutator of the supersymmetry transformations gives 
a translation of the spacetime by = — 2jer^^. This urges us to interpret the spacetime as emerging 
from the eigenvalues of A^. This is a crucial property in order for this model to include the gravitational 
interaction. If the IIB matrix model has a massless particle, it must have a spin-2 graviton. We have a 
good reason to reduce not the four or six-dimensional, but the ten-dimensional Af = 1 super Yang-Mills 
theory. This model has the 32 maximal supersymmetry only if we reduce the ten-dimensional super 
Yang-Mills theory. The maximal supersymmetry is an essential aspect for the theory of gravity. 



2.2 Interaction between the BPS object 

We next review the interpretation of the BPS objects (such as the D-brane) in terms of the IIB matrix 
model. We start with considering the classical equation of motion of the IIB matrix model: 

= 0. (2.47) 

Here, we set the fermion tp to zero. The corresponding equation of motion for the type IIB superstring 
theory is given by 

{X^,{X^',X''}} ^0. (2.48) 
In terms of the type IIB superstring theory, this has a solution representing the Dl-brane: 

Xo = Tt, Xi^^a, X2^X3 = --- = Xg^0, (2.49) 

where T and L are the compactification radii of X'^ and X^, respectively. The parameter r and a take 
values < T < 1 and < a < 2tt. The Poisson bracket is given by 

L TL 

{Xo, Xi} - fOi^o^oaiXi = r X — = — . (2.50) 

ZTT ZTT 

The corresponding commutator for the matrix regularization is 

Tr TT 

-i[A^.A^]^^^{X,,X,} = —. (2.51) 

This commutation relation is realized only for the infinite-size matrices (we can see that this is impossible 
for the finite-size matrices by using the cyclic rule of the trace Tr). Namely, we have the following solution: 

where q,p are the infinite-size canonical pair satisfying [q,p\ — i. 

We discern that the solution H2.52|) is a BPS saturated state in the following way. We substitute the 
solution H2.52|l . along with A2 = A3 = ■ ■ ■ = Ag = and ■(/; = 0, into the supersymmetry transformation 
in the previous section to obtain 

4'V = e, (2.53) 
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The homogeneous supersymmetry can be clearly canceled by the inhomogeneous supersymmetry, when 
we properly choose the parameter ^. In this sense, the Dl-brane solution of the IIB matrix model H2.52|l 
actually preserves half of the supersymmetry, and thus is a BPS saturated object. 

By the same token, we can establish the solutions of the IIB matrix model that correspond to the 
-D3, D5, ■ ■ ■ brane, by taking the similar canonical pairs as 



[^2,^3] OC 1 



NxN, 



[X4,X5] OC 1 



NxN, 



(2.54) 



This nicely corresponds to the fact that the type IIB supcrstring theory accepts the Dl,3,5,- • -branes. 

We next see that the bosonic matrices of the IIB matrix model describe the multi-body system. 
We can easily build the solution for the two Dl-branes located at X2 = ±| by taking the following 
block-diagonal matrices. 



Ao = 



_ ( V27rni 





Ai = 



V27rni-' 





A. = 



A3 = • • • = = 0, 



(2.55) 



where ni_2 is the size of each block, and satisfy = = i- 

By the same token, we can build the Dl — Dl brane system through the following classical solution. 



An 



V27rni 









Ai 



V27m 





-p 



A. 







A3 = • • • = Ag = 0, (2.56) 
In this way, the IIB matrix model describes not only the single D-objects , but also the multi-body 




Figure 1: The two parallel Dl-branes, represented by the classical solution of the IIB matrix model (|2.55|l . 

system by allocating the plural matters on the block-diagonal components. In this sense, there is a great 
difference from the action of the D-instanton action, which appears to be the same action. Unlike the 
D-instanton action, the IIB matrix model accepts the multi-body system and thus is a second-quantized 
action. 

Not only do the matrices of the IIB matrix model accommodate the multi-body system but also 
the interaction is embedded in the same matrices. Namely, the off-diagonal part can be interpreted as 
representing the interaction of each block. 

We next investigate the interaction between the D-objects by the one-loop calculation. To this end, 
we derive the one-loop effective action. We separate the matrices with the classical solution and the 
fluctuation as 

Af,=p^,+a^, tp = x + 'P- (2.57) 

Here, and x &re the classical solution, and and ip are the fluctuation. Since the IIB matrix 
model is invariant under the SU{N) gauge transformation H2.43|l . we need the gauge fixing. Here, we 
adopt the background gauge 

-z[p^,A^]=0. (2.58) 



10 




Figure 2: The ofF-diagonal parts represent the interaction in the IIB matrix model. 



The corresponding gauge-fixing term £ and the ghost term S'gj-^Qgi^ are respectively 

Sg_i=-lTr[p^,A^^]', Sg^^^^ = -Tr{[p^,bW\c\), (2.59) 

where c and b are the ghosts and the anti-ghosts respectively. In the following, we set the parameter g to 
1. Then, we extract the quadratic term of the fluctuation for the total action ^(total) _ 5 + S'„ f +'S'g}]^ost 



4^°*''^^ = -Tr (^[p^, a,]' + b^,p.]K, a"-] + l^T^ip^., ^] + [p^, b][p^, c 



= Tr Qa'^(p2^^^ _ 2tF^,)a'' + ^^T^P^^ - bP\j . (2.60) 

Here, we introduce the following adjoint operators 

P^X = (adp^)X = [p^,X], F^,X ^ {sid{^i[p^,p,]))X ^ -i[P^,P,]X, (2.61) 

where the last equality for follows from the Jacobi identity of the commutators. Integrating the 
fluctuation, we obtain the following effective action. 



W = —log J dadbdcdifi exp{— 82^^^^^^ ) 



T;Tr ( trioxiolog(P'5^. - 2iF^,) - \ogP^ ~ -tnexw^og P^ + -F^.r'^-^^ ) ) (2.62) 



1„ / . ^ . . „o 1 . f ■? i ^ „„„i + r« 

-^triexie\oglP' + - 

Here, the trace Tr is for the N x N matrices, while the trace trioxio and triQxw £^re respectively for 
the ten-dimensional vector indices and the ten-dimensional gamma matrices ( is the Weyl projector, 
and effectively the size is not 32 but 16). 

We evaluate the interaction between the diagonal blocks using the one-loop effective action (|2.62|) , and 
elucidate that the IIB matrix model includes the gravitational interaction. We consider the backgrounds 
having a block-diagonal form: 

A,=p^^dmg{p^^\p'^^\---). (2.63) 

Here, each block p), represents the separate D-object and is an rii x rii matrix. We decompose the 
classical background between the trace and the traceless part as 

=4'^ln.x„.+P«. (2.64) 

Since the eigenvalues of the bosonic matrices are regarded as the spacetime coordinate in the IIB matrix 
model, the trace d/j should be identified with the collective coordinate of the D-objects. We assume that 
the classical D-objects are separated with each other so that the distances d/^ — d// are sufficiently large. 
The adjoint operator P^ operates on the components of the matrix X^^'^'> as 

(P^xf^n ^ (rfW _ - x(-^)p(f) = ((dW - d(f)) -f p^f -f pgf )x(-^), (2.65) 
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where P^^'^-'x'*'^^ = p"jl^ X^"^'^^ and -Pr'^"* = —X^^'^^p^^'^ represent the operation of from the left and 
right, respectively. Similarly, we simplify the notation of the operation of the adjoint operator F^^. The 
background field strength is rewritten as 

f,, = diag(^[p«,rf)],^[p(2),p(2)],. . .) ^ ^fWJ^^.. .). (2.66) 
The operation of i^^^ can be likewise decomposed into the left and right operation as 

= f\^iX^'=) - ^ {f'^^II + f'^^^X^''^ (2.67) 
Since the left and right operation are independent, we have 



TrO= }^TrO'-l^"TrO'^''>. (2.68) 

We obtain the one-loop interaction for these D-objects. To this end, we expand the effective action H2.62|l 
with respect to the power of Ff^. The first bosonic term is expanded as 

^Tr (trioxio \og{P^S^, - 2iF^,)) 
= TrlSlogP +jp^F^y—F^^ — ——F^y—F^p—Fpf^ — 2—Ff^y—Fiyp—Fp^—Fy^^\+0{P ). 



Whereas, the third fermionic term of 12.62|l is expanded as 



rTr (tri6xi6 log (P' + l^ i^^.r''"- 



9 4il_l„l„ 1 



Tr ( 4 log P ^2 ^IJ.!^ p2 ^"P- + 3 p2 ^A'l' p2 ^"P p2 ^PP + p2 ^P'' p2 p2 ^PX p2 ^XP- 



11 1 1 1 11 1 1 1 

'^'p2^P''~p2^PXjy2^P''J}2^PX + 2J^^'"'J^^'^'^J}2^PXJ^^PX 

~'^Jy2^P'^j32^PP~p2^X'^~p2^XP ) +^{P °)- 



In this derivation, we utilize the following formulae of the gamma matrices: 
^^^YPit^2Y'^iu2^ = -2ir(?7['^i[''i7/^^l''^l) 

^^(pPlP2p'>ll'2pPlP2pXlX2-) _ ^^(4j^[Pl[fl jyP2]!^2]y^[pl[XljyP2]X2] _j_ gj^[Pl [Pl ^P2] P2] ^[I'l [Xl ^''2] X2] 

_lQj^[Pl [l^l ^[Pl [Xl ^P2]X2] j^!'2]X2] _ 15jy[Ml [Pl y^P2] [Xl jy['^lX2] ^l'2]X2] 
_|_X5^[Mi[!^1^[P1 [Xl-^P2]X2]^1^2]p2]^_ 

Along with the second ghost term — Tr log P^, we discern that the effect of the boson and the fermion 
cancels up to the 0{P^^) order. This cancellation is ascribed to the supersymmetry of the JIB matrix 
model. Due to this cancellation, the leading effect of the one-loop interaction is of the order 0{P~^). 
Namely, the interaction complies with the power law ~ (where r is the distance of the two D-objects, 
and the dimensionality is d = 10). In this sense, we regard this one-loop interaction as the gravitational 
interaction. This is an important evidence that the IIB matrix model describes the gravitational inter- 
action, and this is a consequence of the supersymmetry. Especially, the contribution of the block is 
expressed using the notation for the left and right operation as 

\K/(i,3) — Vi. Tr- ( —A( f^^^f f'') _ s/'f(*' J- 9i'f(*' J- /"ft*) ft*) f'*) 1 

— 4(|^(i) _ |^{i))8 \^ t^v Ji'p J px Jxpl °\J t^v J pxJ t^xJ pi' ) ~'~ ^yj pv J nv J pxJ px ) ^ yj t^i^ J PX J t^t^ J PX 1) J 

1 ( _A( fU) _ s( fU) fO) fU) fU)\ 1 of fU) fU) fU) fU)\ , ( fU) fU) fU) f(J)\\\ 

4((i(*) — \ 1^" •' ^ P ■' Px ■> xp I °\J p-v J px J px J P" / ^ P" J pv J px J px ) ~ yj pv J px J P" J px )) J 



1 



■ 4(rf{«) _ dU)r (-^^^'■(■^M.'/ip )^-(/mx /xP ) + 6Tr(/«/«)Tr(/(i'4i')) + - (2.69) 

The tensor structure indicates that the first and second term of the last line of (|2.69|l represent the 
graviton and dilaton exchange respectively. This result augurs very well for the IIB matrix model to be a 
bona fide framework for the gravitation interaction. This argument is limited to the background gauge, 
and it is an interesting future work to extend this argument in a gauge independent way. 
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2.3 Interpretation of the difFeomorphism invariance 

In this subsection, we review the interpretation of the diffeomorphism invariance on the target space 
for the IIB matrix model. In [14], it has been pointed out that the invariance under the permutation 
of the eigenvalues S'at is interpreted as the diffeomorphism invariance of the low-energy action. Since 
the permutation invariance is of course a subgroup of the SU {N) symmetry, this means that both 
the diffeomorphism invariance and the gauge invariance emerge from the SU (N) invariance! This is 
a surprising aspect of the IIB matrix model, in that the unification of these two symmetries has been 
achieved in a natural way. 

In discussing the dynamics of the spacetime, we expand the IIB matrix model around the diagonal 
background 

p^=diagixl,xl,---,x^), x = diag(x\x',---,x'^), (2.70) 

where these satisfy the constraint X^iLi = and X^iLi X* — 0, because they belong to the SU{N) gauge 
group. We delegate the details to the references [11, 14, 16,75], but after we integrate the fluctuation at 
the one-loop level and further the fermion zero mode Xj we finally obtain the following interesting result: 

y'dXdxexp(-5°^^-l°°P[a;,x])= Yl J dXW[X;G]. (2.71) 

G: graph 

Here, G denotes the graphs that connects the eigenvalues X, and W is the Boltzmann weight for the 
graph G and the configuration of the eigenvalues X. An explicit calculation also indicates that the 
dependence on the two different eigenvalues and x^ contributes at the order (a;* — x-')~^'^, when and 
x^ are connected by the graph G. 

The important property of the integral (|2.71|) is that it is invariant under the permutation of the 
eigenvalues Sn- We note that this is of course not the case with each graph G, because the way to connect 
the eigenvalues by the graphs impairs the Sn invariance. Nevertheless, this symmetry is retrieved when 
we take a summation for all the graphs. This phenomenon is somewhat reminiscent of the dynamical 
triangulation approach to the quantum gravity, in which the diffeomorphism invariance is retrieved by 
summing over all the triangulation. 

Here, we see how the permutation invariance of the eigenvalues leads to the diffeomorphism invariance. 
In a sense, this is a very natural interpretation, since we identify the eigenvalues with the spacetime 
coordinate. To elaborate on this viewpoint, we consider the scalar field (j)^ propagating in the distributed 
eigenvalues. Namely, we consider the following effective action as an example; 

^ - E ^-^^^n-' - + E (2-72) 

f{x) is a function damping fast sufficiently at infinity. We introduce the density function p{x) = 
^j5'^^°^(x — X*). Along with the continuous function 0(a;) that satisfies 0* = 0(x*), the action (|2.72l) is 
rewritten as 

S^ I dxdy{p{x)p{y))^-^^^^^^^f{x-y)+m f dx{p{x))^{xf . (2.73) 



Here, we take an average with respect to the eigenvalue configuration X and the graphs G. In order to 
see clearly the correspondence between the gravitational background, we further rewrite the action as 



S ^ ]- I dx{p{x)) 



dy{p{y)){x - y)^,{x - y)„f{x ~ y){l + c{x, y)) 



d^'(t>{x)d''4>{x) 



2 „ 

+ rn j dx{p{x))(jp- {x) + ■ ■ ■ . (2.74) 

Here, we normalize the density correlation as {p{x)p{y)) = {p{x)) {p{y)) [1 + c{x,y)). This urges us to 
identify the eigenvalue density with the gravitational background as 

V^e-*(^) ~ (p(x)), [ dy{piy))ix~y)^ix~y),f{x-y)il + c{x,y)). (2.75) 
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Nextly, we see how the diffeomorphism invariance is reahzed. Since the IIB matrix model itself is 
SU{N) invariant, the model itself is of course invariant under its subgroup Sn- Under the permutation 
— > x'^'^*^ for a e Sn, the fields transform to 0*^^*). We extend the permutation of the eigenvalue to 
the transformation of the continuous transformation x S,{x), where ^(a;) is a continuous function such 
that ^(x*) = x"^^'-*. It is trivial that (j){x) is subject to the general coordinate transformation as a scalar 
field. In addition, the background metric given in H2.75f) also receives a transformation as a rank-2 tensor, 
when the function f{x) decreases rapidly and has a support only near the origin x = 0. In this sense, we 
can fairly interpret the Sn invariance as the general coordinate invariance of the low-energy limit. 

While the IIB matrix model has only a fiat noncommutative background and depends heavily on 
the fiat metric, we have an interpretation for the curved space background. A nontrivial background is 
induced dynamically through the condensation of the graviton. 

2.4 Alternative totally reduced models 

In this subsection, we review the attempts for the generalization of the IIB matrix model. The IIB matrix 
model augurs well for the unification of the gravity. The Af — 2 supersymmetry is a fundamental object 
for the theory of gravity. In the preceding two sections, we have reviewed some of the evidences for the 
gravitational interaction: the graviton-dilaton exchange in the one-loop calculation and the interpretation 
of the diffeomorphism invariance. However, the IIB matrix model suffers a serious drawback in that it 
has only the flat noncommutative background. However, if a matrix model is to be a bona fide framework 
to describe the gravitational interaction, it should accommodate a curved space background in a natural 
way. In this section, we introduce several generalizations of the IIB matrix model, so that they should 
incorporate a curved space background. 

2.4.1 The IIB matrix model with the Chern-Simons term 

Firstly, we introduce the matrix model with the Chern-Simons term, defined by the following action [25] : 

^ = (^-\[A^,A.][A'^,A-'] + ^e^^pA^^A-'A^ + ^^a^[A^,,p]^ . (2.76) 

Here, the indices fj,,v, p, - ■ ■ run over 1, 2, 3 and this model is defined on the three-dimensional Euclidean 
space. A^ are the three-dimensional bosonic vectors and ip is the three-dimensional Majorana fermion. 
This is also a totally reduced model like the IIB matrix model, and ^4^ and tjj are promoted to the N x N 
hermitian matrices, cr^ are the Pauli matrices and in the following a^^i, denotes (Jy^^ — ■^[<Jf^, (Ji,] = ie^^jpCTp. 
This model has the 50(3) rotational symmetry and the SU{N) gauge symmetry. 

In [56] , Tomino performed the calculation of the path integral for the N — 2 case as a toy model, and 
elucidated that the path integral converges for the N = 2 case (while the bosonic model 1)5. l|l which we 
analyze in Sec. 5 does diverge). This is in contrast to the IIB matrix model without the Chern-Simons 
term, in which the path integral diverges for d — 3 and arbitrary N. Austing and Wheater in [57] 
corroborated that adding the Chern-Simons term does not affect the convergence as long as the original 
path integral (without the Chern-Simons term) converges absolutely. 

The interesting property of this matrix model is that it incorporates the classical solution of the fuzzy 
sphere. Its classical equation of motion is given by 

[A^, [A^, A^]] -I- iae^p^[Ap, A^] = 0. (2.77) 

The fuzzy-sphere solution is given by the A^-dimensional irreducible representation of the 51/(2) Lie 
algebra L^: 

Ap = aL^, where [Lp,Lt,] = iep^^pLp. (2.78) 

This is regarded as the sphere, firstly because it satisfies the relation 

_ 1 

Al+Al+Al = RHNxN = a^^^lNxN- (2.79) 

The relation H2.79|l determines the radius of the sphere as — a'^ ^ . In addition, the eigenvalues 
of Ap are distributed sphere-like. In this sense, this classical solution is interpreted as the sphere in the 
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three-dimensional spacetime. That this matrix model has the S*^ fuzzy sphere classical solution is a great 
advantage. While the fuzzy sphere is nothing but a simple manifold, it is of significance as a prototype 
of the curved space background, which is a fundamental principle of the general relativity. 

This matrix model also incorporates the J\f = 2 supersymmetry. Like the IIB matrix model, this 
model also has the homogeneous and the inhomogeneous supersymmetry. 

homogeneous: J^^Up = iecr^^V, S^J-''^: ^ ^{[Af,, A^] ~ iaef,^pAp)a'"'e, (2.80) 
inhomogeneous: (5f ^A,, = 0, S^ip = ^. (2.81) 
This model has the similar commutation relation to that of the IIB matrix model: 

0, [Sil\6il^]^ = 0, (2.82) 
0, [4f,4f]V' = 0, (2.83) 
-ze-a^^, [<5(i),4'V = 0. (2.84) 

The proof of these commutation relations goes in the same way as in the IIB matrix model, and we give 
only the proof of ()2.82|l . Using the formula cr^^a^p = Sp_^ap — Sp_p<7^, we obtain the commutation relation 

[5il\sil^]Ap = t[Ap,X]+ep,pe,Ap, (2.85) 

where A = 2i(e2<yp£i)Ap and 0^ = 2i(e2cr^ei). For the fermion, we obtain the commutator by taking the 
difference likewise: 

= V - ([eiCT^V^, A^]CT^^e2 + a(eiCTpV')o-pe2) + (d ^ £2). (2.86) 
We use the following Fierz transformation whose detailed proof we give in Appendix I A . 1 . 5l 

ieiaptp)af,t,e2 = -(eicr^e2)V'' + (rank-0 terms), {eiap'ip)ape2 = -(eiO'pe2)cr^?A + (rank-0 terms). (2.87) 

In (|2.87|) . we omit the rank-0 contribution, since we are interested in the commutator [6el\ Sel^jtp- Then, 
we obtain the commutation relation 

[,5W,4i)]V = *[^,A]-za^^V- (2.88) 

The difference from the IIB matrix model is that the commutation relation H2.88|l is satisfied without 
using the classical equation of motion. The relations H2.85|l and (|2.88|) indicate that the commutators 
vanish up to the S0{3) rotation by 9p, as well as the gauge transformation by A. The appearance of the 
SO{3) symmetry is a novelty of this matrix model. 

We likewise take the linear combination of the supersymmetry as 

^(1) =<5W+<5(2), ^(2) ^z(5(i)_5(2)). (2.89) 

Their commutation relation gives the translation of the bosonic matrices Ap as 

[^("\4^)]V-0, [S^,''\~df^]Ap^-2tieapOd''^. (2.90) 

The fuzzy-sphere classical solution is a BPS saturated state. Plugging the classical solution (|2.78ll and V' = 
into the homogeneous transformation (|2.80|l . we immediately find that the homogeneous supersymmetry 
vanishes. Whereas, the inhomogeneous supersymmetry H2.81II survives. In this sense, the fuzzy-sphere 
solution (|2.78(l preserves half of the supersymmetry. 

We next discuss the expansion of this matrix model around the fuzzy sphere solution (|2.78|) . In the case 
of the IIB matrix model, we derive the noncommutative Yang-Mills theory by the expansion around the 
flat noncommutative background [17]. This derivation of the noncommutative Yang-Mills theory is quite 
natural, because the matrices, per se, are noncommutative objects. The studies of the noncommutative 
field theory have become explosively popular since Seiberg and Witten [18] elucidated the relation to the 
superstring theory. In the noncommutative spacetime, we impose the noncommutativity on the spacetime 
as [xp, Xi\ = icpy. This relationship is reminiscent of the canonical commutation relation of the space and 



15 



the momentum in the quantum mechanics. Just as the spacetime and momentum have the uncertainty 
principle, the noncommutative spacetime incorporates the uncertainty of the spacetime. In the context 
of the superstring theory, the noncommutativity of the spacetime is introduced by turning on the B 
field [18], which gives the noncommutativity 

Cau = 2Tra'( ^- I ^-(2na')( ^ B ^ ) . (2.91) 

This leads to the uncertainty of the spacetime at the string scale O(Va') ~ 0(ls)- 

When we expand the matrix model (|2.76() around the fuzzy sphere background (|2.78() . we obtain a 
noncommutative Yang-Mills theory on the S'^ sphere. The expansion around the sphere is given by 

= aL^ + ai?a^, (2.92) 

where is the fluctuation. The fluctuation is expanded by the noncommutative spherical harmonics as 

= ^ ^ af_LlmYlm, = X! X! ^IrnYlmj (2.93) 
m— — l /— rn— — / 

where the noncommutative spherical harmonics Y/^ is given by 

Yi„, - J2 C-a.i'^'i''' ■■■i"'- (2.94) 

a 

The indices run over 1,2,3, and /q™..^^ is a symmetric and traceless tensor. We can understand 
that this is a {21 + l)-dimensional tensor as follows. Firstly, when we allocate 1,2,3 to ai, - ■ ■ ,ai in a 
symmetric way, we have 1+2C2 = ways. On the other hand, the tracelessness condition implies 

that /o™ai02 - ai_2 vanishes. This condition stems from the constraint (a;")^ = Inxn- This eliminates 
1C2 = ^^'2 ways. Therefore, the dimensionality of this tensor is _ '^'^ = 2/ + 1. Its 

normalization is defined by 

^Tr{YirnYil^,)^6w6mm'- (2.95) 

This of course corresponds to the expansion in terms of the c-number spherical harmonics 



(f7) = ^ ^ af,i„iYi„i{fl), = X! i^irnYim{fl), whcrc (2.96) 

l—Q rn— — l /— m— — l 

-T-Yi„.Y{:^, = — / smedOYi^Y;^, = SwS. 

, , «/ >/ 



(2.97) 

Now, we can build the mapping rule between the matrices and the functions on the 5*^ sphere. The rule 
is listed below: 

AT-l I oo I 

Op = ^ ^ a^LlrnYhn ^ a^j(f^) = X! X! ".f^lmYlmi^), (2.98) 
l—O in— — I rn — — l 

ah^ a{VL)i.h{n). (2.100) 

We need to look more closely at the star product H2.1()0|l . Firstly, the correspondence H2.98|l is rewritten, 
using the orthogonality condition 12.95|l . as 

N~l I 

^ a,{n) = ^ E E Tria^Y^jYi^in). (2.101) 

1=0 m=-l 
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Here, we require that the maximum value of Hs A'' — 1. This leads us to define the star product as follows. 
Here, we require that / + /' should not exceed A'' — 1. 



= ^EE E / dn'dn'%r^,{n')Yii^^4n)a{n')b{n'')Yi„,{n)T^^^^ 



Im I'm' l"m" 



(2.102) 

Then, the action (|2.76() is mapped as follows: 



a 



a'^(N^-N) a'^R^N f dQ ( 2i fl 1 i 



aRN f dn fl-^^^^^^^-^^^^^^\ ^2.103) 



2^2 J 47r \R 

where we define the field strength Fp^ as 

1 i 
Ffj^i^ = 7;([^A"«!^] ^ [Li^,ap]) ~ —ep„pap + [ap,a^]. (2.104) 
K K 

The subscript ^ indicates that the product is defined by the noncommutative star product. Therefore, 

the commutator survives even for the U{1) abelian gauge group of the Yang-Mills theory. 

This gives the relationship between the coupling constant of the Yang-Mills theory gyM and the 

artifact of the coupling constant of the Yang-Mill theory before being reduced as 

We go on to the symmetry of this Yang-Mills theory. Firstly, the SU (N) symmetry is interpreted as the 
gauge symmetry. For the element of the SU{N) Lie group U, it is expanded as ?7 = exp(iA) ^ l + iX+ ■ ■ ■, 
where A is an element of the su{N) Lie algebra. The gauge transformation for Ap is given by 

6Ap = [Ap,iX\ = [aLp + aRdp,iX]. (2.106) 

This transforms the fluctuation as 

6dp = ^a[Lp, A] -f i[ap, A]. (2.107) 
R 

This is immediately translated in terms of the Yang-Mills theory as 

5ap^ ^Lp\ + i[ap,\]^. (2.108) 
R 

Secondly, we note that this noncommutative Yang-Mills theory preserves the following M = \ supersym- 
metry. 

5sap = ^eap^, Ss^ = '-^^Fp^a'^'^e. (2.109) 
ari Z 

We have hitherto discussed the perturbation around the irreducible representation of the fuzzy sphere 
(|2.78|) . In this case, we obtain a U{1) abelian noncommutative Yang-Mills theory. The matrix model 
(|2.76|) also accommodates the fuzzy sphere solution of the reducible representation 

/ (^Lp \ 
aL,, 



Ap = a{Lp ^ 1 



kxk) 



(2.110) 



V ^Lp J 

When Lp is the nxn representation of the SU{2) Lie algebra, the size of the matrices N satisfies N = nk. 
The perturbation around this reducible representation gives the U{k) nonabelian theory. 
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2.4.2 Other attempts for higher-dimensional fuzzy sphere solution 

We next review the extension of the above discussion to the higher-dimensional fuzzy sphere classical 
solution. To this end, there are two strategies. First is the addition of the higher-dimensional Chern- 
Simons term, and second is the addition of the tachyonic mass term. 
The first attempt is achieved by defining the following action [40, 50] : 



±[lTr[A„A^]- 



(2.111) 



Here, we define this matrix model in the d = {2k + l)-dimensional Euclidean space. The indices ^, ^, • • • 
run over 1, 2, • • • , (2fc-|-l). The supersymmetric Yang-Mills theory can be defined only for the d = 3,4, 6, 10 
case, and the supersymmetric matrix model cannot be defined except for the three-dimensional case as 
we have reviewed before. Hence, we focus on the bosonic matrix model in the following. Except for 
the fc = 1 (three-dimensional) case, the convergence of the path integral is not yet manifestly proven. 
Nevertheless, we surmise that the fuzzy-sphere classical solution might be a metastable state, in that it 
is barriered by a high potential. We attribute this metastability to the similar reasoning to the (j)^ matrix 
model. As we explain in Appendix |21 the matrix model retains the metastability near the origin 
despite its unbounded-below potential. 

The illuminating feature of this higher-dimensional matrix model is that it incorporates the classical 
solution of the 5^'^ fuzzy sphere. The equation of motion is given by 



(2.112) 



It has the following classical solution of the S fuzzy sphere [10,31,43] 



= f where 



^(2fe) 



L2'=x2'= Jsym " 



(I2 



sym- 



(2.113) 



We define the symmetric tensor product more explicitly, for A; = 1,2 case (namely, for the S'^ and S"^ 
fuzzy sphere). We denote the matrix element using the orthonormal state \i) and \j) as {A)ij — {i\A\j). 
The usual tensor product (two- fold product, for brevity) is defined by 



{ii,i2\Ai<» A2\ji,j2) = (ii\Ai\ji)(i2\A2\j2), where |ji,j2) = |ji)|j2). 
On the other hand, the symmetric tensor product is defined as 

sym(ii,«2|-4 «) B|jij2)sym, where |ji,j2)sym = -(|ji)|j2) + |j2)|ii)) 



(2.114) 



(2.115) 



The coefficient n is determined, so that the symmetrized state may retain the orthonormality. Here, we 
give a more explicit definition for the two- fold tensor product of the Paulian matrices (fc = 1 case). For 
the two-fold tensor product, the symmetrized orthonormal state is given by 



|1,1) = |1)|1), |2,2) = |2)|2), |l,2) = |2,l) = i=(|l>|2> + |2)|l)). 



(2.116) 



This difference comes from the fact that the states |1, 1) and |2, 2) are symmetrized ab ovo and hence that 
we do not need symmetrization or normalization. We elaborate on these definitions in terms of explicit 
representation in terms of matrices. The usual tensor product (12 .11411 is rewritten as 



A<RiB 





|J1 = 1,J2-1) 


jjl =2,J2 = 1) 


bi = 1,J2 = 2) 




2) \ 


in 


= l,i2 


= 1| 


aii5ii 


ai2&ii 


011612 


012612 




(^1 


= 2,^2 


= 1| 


021^11 


022611 


021612 


022612 




(n 


= 1,«2 


= 2| 


011621 


012621 


011622 


012622 






= 2,^2 


-2| 


021^21 


022621 


021622 


022622 


/ 



(2.117) 
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On the other hand, its symmetrization is rewritten as 





/ 


1, l)sym 


1, 2)sym 


2, 2)sym 


{A (E) S)sym = 


sym(l, 1 
sym(l, 2 
^ sym(2, 2| 


aii&ii 

(a2ifcii+aiib2l) 


(ai2fcll+ailfcl2) 
V2 

(a22bll+a2l6l2+ai2fc2l+ailft22) 


ai2&i2 

(022612+012622) 


a2ifc2i 


. (V2)' , 
(022621+121622) 
V2 


022^22 



The coefhcient A is related to the radius parameter of the fuzzy sphere a as 



(2.118) 



(2.119) 



where the coefficient nik is given later in the self-dual relation relation H2.12(i|l . 

The S"^*^ fuzzy sphere is mathematically more involved than the S*^ fuzzy sphere. We introduce the 

(2k) 

properties of the higher-dimensional fuzzy sphere solution. Firstly, the gamma matrices is defined 

by the Wick rotation of those for the Minkowski spacetime given in Appendix^ The definition is given 

by 



rf =ap, (p= 1,2,3), Ff)=Ff 



(2fe) 
2k ~ 

(2k) 
2k+l 



l2fe-ix2'=-i 
- l2'=-i x2'=-i 



0-2 



1 0-1 = 




(2.120) 



{2k) 



where p runs over p = 1, 2, • • • , 2fc — 1. Namely, we identify «Fo in the Minkowski space with T^f. 
Euclidean spacetime. 

The size of the representation (|2.113l) is given in [31] as 

(n + l)(n + 2)(n + 3) (n + l)(n + 2)(n + 3)2(n + 4)(,i + 5) 

Ni = {n + 1), A'2 = , A'3 — 



the 



6 

(n + l)(n + 2){n + 3)^{n + Af{n + 5)^{n + 6)(n + 7) 
302400 ■ 



360 



(2.121) 



We see a grave difference between the k — 1 (three-dimensional) case and otherwise. For fc = 1, the 
representation (|2.113|l can be defined for an arbitrary integer-size matrices (of course, except for 1). 
However, this never holds true for k > 2 case. For example, the representation (|2.113|l is possible only 
for a limited integer size; namely N2 — 4, 10, 20, • • •. 

The important property is that G^^*"' does not generically close with respect to the commutator for 
k > 2. This is also a tremendous difference between the three-dimensional case and the higher-dimensional 
case. Here, we denote G^i*^^ = [G^^''^ G^^'^^]. These satisfy the following relations: 

Gf )Gf ) ^n(n + 2fc)l^,xiv,, (2.122) 

- -8kn{n + 2fc)l^,xiv., (2.123) 

[Gll'^', Gf )] = 4(-J,,pG(2'=) + <5.pGf )), (2.124) 

[Gil'lGil'^] = Ai5.,Gif + 5,,G(^/) - J,,Gl^/) - 6.,GI'^^ (2.125) 



In addition, G 



{2k) 



71X ' "l-'-X^vp 

satisfies the following self-dual relation. 



.G(f G(f)...G(2,^)=™,G„ where 



mi = 2j, 7712 = 8(n + 2), mk = —2ki{n + 2fc — 2)mk- 



(2.126) 
(2.127) 



Before verifying these properties, we start with convincing ourselves that the notion of the higher- 
dimensional fuzzy sphere solution is a generalization of the SU (2) Lie algebra. To this end, we focus on 
the relation between the k = 1 (three-dimensional) case and the SU (2) Lie algebra. 
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Firstly, the n-fold symmetric tensor product for A; = 1 is identical to the (n+l)-dimensional irreducible 
representation of the SU{2) Lie algebra. We explicitly give the matrix element of the n = 2 case as 



(2) _ 



[((Ti (g) 12x2) + (12X2 «) fJi)] 



sym 



/ 


1 






/ 


1 




1 
















1 








1 


+ 


1 








V 


1 


0^ 


sym 


V 


1 


/ 



1 
= V2 I 1 1 
1 



sym 



G. 



(2) 



G 



(2) 





512x2)4 


-(12 




( ^ 


~i 







i 
















—i 




I 


i 







512x2)4 


-(12 




/ 1 













-1 










1 







I 





-1 



' 0'2)Jsym 
/ 



V 



sym 
2 ® cr3)]sym 



/ 




sym 



1 



sym 



V 









-1 







= 2 



-1 / 







-1 



(2.128) 



sym 



These are clearly identical to the 3x3 representation of the SU{2) Lie algebra. In addition, the 2-fold 
tensor product obeys the following important relation: 



3 

E 



(cr^ (81 cr^)sym = 
















/ 











-1 

















\ 








1 















1 












-1 













1 








+ 







1 










+ 








-1 







1 








) 


sym 


I 


-1 











) 


sym 


I 








1 






L3x3- 



sym 
(2.129) 



Utilizing this relation, the radius of the fuzzy sphere is derived as follows: 

)sym - 



[G^l^yf ^ {al®l®---® l)sym + + (iDsym + (ct^ ' 



1) 



sym 



n terms 

r2 



n(n — 1) terms 



= (3n + n{n - l))lArixAfi = (^1 - i)iNixNi, 

where we used (|2. 12911 in the second equality and recall that A^i 



(2.130) 



for fc = 1. This is identical to the Casimir of the J^, which satisfies = — ^ 



n + 1. This proves the relation (|2.122|l 
^ Here, the relation 



between G/, and Jf^ is given by iC^i — 



The self-dual relation (|2.12t)|l clearly reduces to the commutation relation of the SU{2) Lie algebra: 

(2.131) 



^(2)(^(2) _ 1 r(^(2) (^(2)1 _ 1 (2if G'-^h - 27(7(2) 



We next go into the proof of the relations (|2. 122(1 ~ (I2.126|l . We immediately discern that 12.1241) and 
(|2.125|l are verified from the commutation relations of the gamma matrices. (|2.122f) can be substantiated 
by repeating the same argument as for (|2.130|) in the fc = 1 case. In promoting to the higher-dimensional 
case, we note that the similar relation to H2.129|l holds true of the general case; namely 



2fe+l 

E 



sym 



(2.132) 



For fc = 2. this relation can be verified as 



E(rf^ 



'r|f'=))sym 
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E 



((ct^ «) <J3) (ct^ «) cr3))sym + ((1 <^ era) (g) (1 cr2))sym + ((1 «) ai) <E) {1(g) cri))syni 



by using (|2.129|l twice. We can verify the higher- fc case recursively in the same way. Once we prove 
(|2.132|l . the proof of H2.122|) goes in the same way as for (|2.13U|) . 
The relation (|2.123|l is verified by the following calculation: 

G^l'^Glfj^^ = pr'^'' ^ 1 • • • ® l)sym + • • f 

= (ir^T'^" ® 1 ® • • • ® l)sym + 4((r^r'' - Sf"") ® (PT" - 5''") ® l ® • • • ® l)sym + • • • 



n terms rt(n - 1) terms 

4(-n2fc(2fc+ 1) l)2fc)l = -8fcn(n-|-2fc)l. (2.133) 

Lastly, we prove the relation H2. 126(1 . The brutal- force proof up to fc = 4 is given in [46]. However, the 
Gordian knot is cut in one stroke in [50] by noting the following reduction law. We consider the algebra 

at the north pole of the fuzzy sphere; Gg^^j^ ^ n. Then, G^^*^' are the generator of the SO{2k) rotation 
around the north pole, where the indices a,b, ■ ■ ■ run over 1, 2, • • • , 2fc. For p, g, • • • = 1, 2, • • • 2fc — 1, the 
elements for the 5^*"' fuzzy sphere can be identified with those of the S'^^^'^ fuzzy sphere; 

iG^S - Gi'"-'^ , G(f ) ^ G(^/-^) . (2.134) 
This gives the recursive relation of the coefficient as 

'^*='-'2fe+i - - eMi---/'2fc(2fe+i)G^, ' ■■■ G)^J - 2fce^i...^,^^_,(2fe)(2fe+i)G;,^^'^ • • • G^^^_^2fe 

= -2fep,...p,,_,G(f-2) ^ „2fczmfc_i(Gg^:f))2 - -2fczmfc_in(n + 2fc - 2). 

(2.135) 

This shows that the coefficient mk obeys the recursion relation 1(2. 127(1 . 

This completes the relations ((2.12211 ^ ((2.126() . We can see that ((2.113() represents the sphere simply 
due to the relation ((2.122(1 . This gives the radius of the S*^*^ fuzzy sphere as 

^Mn + 2k)^ ^2.136) 

The relation l(2.119() has a serious consequence for the radius. For the fc = 1 case, the radius is proportional 
to the coefficient A. However, for fc > 2 case, the radius becomes smaller for the larger coefficient A. 

Nextly, we point out that the S'^^~'^ fuzzy sphere is attached to the S"^^ fuzzy sphere at each point. 
This property is understood by considering the algebra on the north pole G2},_li ^ n. As we mentioned 

before, g')^^'^ are the generators of the SOl2k) rotation. Since identification ((2.13411 gives the algebra 
of the 5*^*^"^ sphere, we see that the 5^'^"^ sphere is attached on the north pole. While we see only 
the north pole, it is clear that this holds true of all the points of the S"^^ sphere due to the rotational 
symmetry. Due to the identification ((2.134(1 . the radius of the S'^^~'^ fuzzy sphere attached to the S"^^ 
sphere is clearly given by 

, a'^n(n + 2k-2) 
R%,-. = (2.137) 

This property can be understood from the viewpoint of the stabilizer group. The stabilizer group of the 
S"^^ fuzzy sphere is S'0(2fc + l)/?7(fc). This is identical to 

SO{2k + l)/U{k) - (S'0(2fc + l)/50(2fc)) X (S'0(2fc)/C/(fc)) 

- (S'0(2fc + l)/50(2fc)) X (50(2fc- l)/C/(fc- 1)). (2.138) 

This solidifies the fact that the S'^^~'^ sphere is attached on the S"^^ sphere. 

We see another grave difference between the S'^ and the higher-dimensional S"^^ fuzzy sphere with 
respect to the classical stability. The classical energy is easily calculated by plugging the fuzzy sphere 
solution 5^^^ into the action as 

^_ kn{n + 2k)Nka^ (\ 1 ^ _ 

^- V U^^fcTlj- ^^-^^^^ 
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In this derivation, we use the relation l|2. 122(1 and H2.12H|I . as well as the self-dual relation (|2.12()|l . The 
classical energy (|2.139f) is negative for the three-dimensional (fc — 1) case, and thus energetically favored 
than the flat background A^^ = diag(2;^, • • • , x^), for which S* = 0. However, the reverse is the case with 
the fc > 2 case. In this case, the fuzzy sphere solution is less energetically favored. 

We have seen the alternative matrix model incorporating the higher-dimensional fuzzy sphere solution 
and the properties of the solution. Another choice is to add the tachyonic mass term to the bosonic part 
of the IIB matrix model as [29] 



S=--TT{-\A^,A,f + \Al\. (2.140) 



1 A, 

This model is defined by the general d-dimensional Euclidean space (here, the indices fj,,!^, - ■ ■ run over 
1,2, ■ ■ ■ ,d). The convergence of the path integral of this matrix model is not clear, while it is balanced by 
the mass term and the quartic term. This comes from the diagonal part of the matrices A^. The diagonal 
part of A^ of course contributes to the tachyonic mass term. Whereas, the quartic term comprises the 
commutator, so that the diagonal part is canceled. This is a serious drawback, compared with the matrix 
model with the Chern-Simons term. 
The equation of motion 

[A,,[A,„A,]] + 2XA^ = (2.141) 

gives a variety of the curved-space classical solution. Due to the commutation relation H2.124|) and H2.125|l . 
(for d > 2k + 1) it accommodates the 5^'"' fuzzy sphere solution 

A, - |g^, (2.142) 

where a and the mass A is now related as 

A = fca^. (2.143) 
The classical energy of the fuzzy sphere solution A^ — ^Gf^ for the action (|2.141|l is 

^=-M!i±p)^. (2.144) 

This matrix model model has a wider class of the non-trivial classical background than the matrix 
model with the Chern-Simons term. While we delegate the explicit definition to the paper [29], this 
matrix model accommodates the classical solution of the fuzzy torus. This property is not shared with 
the matrix model with the Chern-Simons term. In this sense, the matrix model with a tachyonic mass 
term is an interesting attempt for a richer curved space background of the matrix model. 



2.5 Other related studies 

In this brief review, we do not exhaust all the recent developments of the IIB matrix model. Instead, we 
touch on the other important studies of the IIB matrix model to conclude this review. Firstly, in [11], 
the reasoning for the IIB matrix model to induce the four-dimensional spacetime dynamically has been 
found. The novelty is that they related the branched polymer and the eigenvalue distribution of the 
IIB matrix model. In this way, they found that the Hausdorff dimension of the set of the eigenvalues of 
the matrices is four. This is a great discovery, in that the IIB matrix model implies the dynamical 
generation of our real four-dimensional world. This issue is reviewed in detail in [16,75]. 

Another proposal for the dynamical generation of the four-dimensional spacetime is proposed in [36] . 
In this paper, they discussed the Gaussian expansion of the IIB matrix model, by taking the quadratic 
action to be the "test action". They conducted the third-order calculation of the Gaussian expansion, 
and found that the 5*0(10) Lorentz symmetry breaks down to the S'0(4) symmetry. This also serves as 
an important evidence for the generation of the four-dimensional spacetime. In [39,48], they refined the 
calculation by means of the improved Taylor expansion and the notion of the 2PI (two-particle irreducible) 
Feynman diagram. They substantiated the breakdown of the 50(10) Lorentz symmetry to the 50(4) up 
to the seventh order of the Gaussian expansion. 

The IIB matrix model is a successful proposal for the constructive definition of the superstring theory, 
and possesses many exciting properties which solidify the confidence that it is an authentic constructive 
definition. 
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3 Studies of the osp{l\32, R) supermatrix model 

In this section, we discuss the matrix model based on the osp(l|32,i?) super Lie algebra, as a natural 
extension of the IIB matrix model. Firstly, L. Smolin [19] proposed the matrix model based on the 
osp{l\32, R) super Lie algebra as an M-theory matrix model. Smolin discussed the possibility for this 
matrix model to induce both the BFSS model [4] and the IIB matrix model [5]. 

As we have briefly reviewed in the previous section, the IIB matrix model enjoys a lot of successful 
aspects. And, we have also discussed some of the attempts for the extensions of the IIB matrix model 
for a more manifest formulation of the gravitational interaction. Here, we introduce the "supermatrix 
models", which are defined by the "super Lie algebra". There are a lot of illuminating features to the 
supermatrix model. Firstly, the models include both the bosons and the fermions in a unified way, 
since the bosons and the fermions are embedded in the same multiplet. Secondly, as we will see more 
closely in the definition, the supermatrix models accommodate the higher-rank tensor fields, as well as 
the rank-1 bosonic fields A^. Especially, this opens the door to the generalization of the spin connection 
in terms of the large- A'^ reduced models. If we are to describe a gravitational interaction more manifestly, 
we should take this question more seriously. Thirdly, the action is described by the cubic term of the 
supermatrix. This might be a natural aspect for the generalization of the superstring theory, because the 
cubic interaction is the fundamental one in the sense that the four-point interaction is identical to the 
connection of the two cubic interactions due to the conformal invariance. Moreover, this model might 
be solvable in relation to the Chern-Simons theory. The Chern-Simons theory is known to be exactly 
solvable by means of the Jones Polynomial in the knot theory [73] . If we find a correspondence between 
the cubic matrix model and Witten's technique of solving the Chern-Simons theory, we may obtain a 
way to solve the superstring theory exactly in the distant future. This section is based on the author's 
works [26,46]. 

3.1 Definition of the osp{l\32, R) super Lie algebra 

Before entering the investigation of the superstring action, we settle the definition of a super Lie algebra. 
A super Lie algebra is an algebra of supermatrix in which both bosonic matrices and the fermionic 
matrices are embedded in one matrix. Supermatrices possess many properties different from (ordinary) 
matrices, and these properties and the notations are summarized in detail in Appendix. IA.2I 
Let us start with the definition of the osp(l|32, R) super Lie algebra. 

H, If M e osp(l|32, R), then ^MG + GM = for G = 

^ M is traceless with respect to 33 x 33 supermatrix. 
X M is a real supermatrix in that M* = M. 

We confirm that, for the first condition, osp(l|32, R) forms a closed super Lie algebra. Suppose matrices 
Ml and M2 satisfy the condition 

^AfiG + GMi = 0, ^AfaG + GA-h = 0. (3.1) 

If osp(l|32, R) is to be a closed super Lie algebra, we call for the following condition: 

^([Mi, MaDG -t- G[Mi, M2] = 0. (3.2) 

Multiplying G^^ on both H3.1() and H3.2(l from the left, we respectively obtain 

G-^'^ MkG + Mk ^ Q, (3.3) 

G-i^([Mi, M2])G + [Ml, M2] = 0, (3.4) 

where k = 1,2. The proof that osp(l|32,i?) is a closed super Lie algebra is equivalent to deriving H3.4|l 
utihzing H3.3|l : 

n = [G-i'^MaG, G-i^AfiG] + [A^i, M2] ^ [-Afa, -A/i] + [A^i, Afa] = 0. (3.5) 
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This statement, per se, can be satisfied whatever the matrix G may be so long as G has an inverse matrix 

ro 



Here comes one question: Why do we define a metric G as G = ^ This stems from the 

requirement that M should be a real matrix; namely M* = ('"My = M. Let us consider the consistency 
between this reality condition and the very definition of osp(l|32, R) super Lie algebra. Take a hermitian 
conjugate of the definition ^MG + GM — 0. This gives 

G^('^M)^ + M^G^ = 0. (3.6) 

Utilizing the properties introduced in the Appendix. IA.2.31 and the reality condition, this is rewritten as 

= Gt(^M)t + MtGt gHi* + AftGt ^''=^S^n ^t^/* + ^(Af*)Gt gHi + ^MG^. 

(3.7) 

In order for the relationship (|3.7|l to be consistent with the very definition of osp(l|32, i?), we must call 
for a condition G^ = G. And the starting point of osp(l|32,i?) is well-defined if we take a metric as 



i 

The next issue is to investigate the explicit form of this super Lie algebra. This is expressed by 

If M g osp(l|32), then M ^ ^ ^ . (3.8) 

• TO contains only the terms of rank 1,2,5. In other words, m is expressed as 

m = UA,T^' + ^ua.aJ"^'^' + ^^UA,-A,T^'-^'- (3.9) 
2! 5! 

• Here ip denotes ij/^V^ . However, this is equivalent to '^ipT'^ , because we are now considering a real 
super Lie algebra. 



The proof of this restriction goes as follows. Let M be the element of osp(l|32, R) and M — 
• is defined as = ■0'^r° = "^^F" . 



TO tp 

id V 



• Because M is a real supermatrix, to, v and ip are real, while i(j) is pure imaginary (and hence is 
real). 

By substituting this formula into the very definition, we obtain 

^mP" + r"TO ^</. + r"v A „ 

We immediately obtain the relationship between two fermionic fields '0 and (j) from this definition: 

^0 + r> = ^r> + rV = -rV + rV = o. (3.10) 

By multiplying on the both hand sides from the left , we obtain the relationship between ip and (j). 
We immediately note that ^t/iP" — = is equivalent to the equation H3.1()|l . 

The constraint on v is trivial, and v must vanish. On the other hand, the constraint on to is worth a 
careful investigation, m is imposed on the constraint 

^TOP°+r°m = 0. (3.11) 

This is the very definition of the Lie algebra called sp(32), and this statement indicates that the bosonic 
32 X 32 matrix of the osj3(l|32, R) super Lie algebra must belong to the sp(32) Lie algebra. In analyzing 
this supermatrix model, it is more convenient to decompose the bosonic part to G sp(32) in terms of 
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the basis of arbitrary 32 x 32 matrices^. Namely, we decompose the bosonic part m in terms of the 
eleven-dimensional gamma matrices 132x32, T"^, r^^^^^, r^'^^^^" , Y^'-'"^'^ and r'*^ ' "*^, rather than to 
obey the expression in the paper [19]. 

The relationship (|3.11() determines what rank of the 11 dimensional gamma matrices survive. Suppose 
m € sp{32) are expressed in terms of the gamma matrices: 

m^ul + uaT^ + ^UA,A,T^'^' + LuA,A,A,T^'^'''^ + ^UA,...A,r^'-^^ + ^UA,...A,T^'-^\ (3.12) 
The condition H3.11|l is rewritten as 

m = -(r°)-ifTO)r° = r°fm)ro. (3.13) 

Then, performing the following computation for fc = 0,l,---,5, we obtain 

where k is the rank of the gamma matrices. Combining this result with the constraint of the sp(32) Lie 
algebra (|3.13|) . we discern that only the gamma matrices of rank 1, 2 and 5 survive. We are thus finished 
with verifying the explicit form of the osp(l|32, R) super Lie algebra. 



3.2 osp{l\32, R) (nongauged) cubic matrix model 
3.2.1 Action of the cubic supermatrix model 

We first discuss the nongauged action of the osp(l|32, R) super Lie algebra (we postpone the meaning of 
the "(non)gauged action" later). We start with the following action proposed by L. Smolin [19]: 

S - i-2Trf: ((f:M^[M^,M^]]-Mg[M§,M%^]] 

QM=l \ \p=l / / 

= 4 E s^^33x33(M'^M''i^r)Tr(r[^^^^) = -|i|s^r(ArM^A/r). (3.15) 

Here, M is a supermatrix belonging to the osp(l|32, R) super Lie algebra. The indices P,Q,R,--- and 
p,q,r, - ■ ■ respectively run over P,Q,R, - ■ ■ — 1, 2, • • • , 32, 33 and p,q,r, ■ ■ ■ , = 1, 2, • • • , 32. These indices 
are respectively allocated for the whole 33 x 33 osp(l|32,i?) supermatrix, and the 32 x 32 bosonic part. 
Here, we promote each component of the osp(l|32, R) super Lie algebra to the u{N) hermitian matrices. 
In the following, the lowercase trace and the uppercase trace are for the 32 x 32(33 x 33) matrices and the 
N X N matrices, respectively, fabc is a structure constant of the u{N) Lie algebra, and i is necessitated 
in the overall definition of the supermatrix model in order to render the action real. 

This model apparently possesses the following pathological properties [19]. The first is that this theory 
is not bounded from above or below. This pathology stems from the fact that the action is cubic, and can 
be seen by a naive field redefinition M MX. The action is rewritten as S" ^ \^S, and we can set this 
action to be 5" — > ±cx) by setting the parameter to be A — > ±oo. This means that the path integral of the 
theory Z = J does not converge. However, note that this pathology is shared by general relativity. 
This pathological property of general relativity can be seen by Weyl transformation 

Seinste^n = J d^X^R = J d'' XC^i'^^^ ^{R - {d - UJ - fc^M^a^^^A^^). 

Although this pathology of Smolin's proposal is not a happy aspect, this may be regarded as a good news 
because this may indicate that Smolin's proposals may include general relativity by taking a due limit. 
The second problem is that this theory possesses no explicit time coordinate. This is again a property 
shared by the general relativity. We can introduce a time coordinate by expanding the theory around a 
certain background. Once a time coordinate is introduced, we can construct a Hamiltonian of this theory. 

^Actually this is a 32 X 32 = 1024 dimensional basis, because the dimension of this basis is 1 + iiCi + ■ ■ ■ 11C5 = 
1 + 11 + 55 + 165 + 330 + 462 = \{1 + 1)" = 1024. 
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We mention the symmetries of this matrix model. Firstly, this model naturally has the osp{l\32, R) 
symmetry, which is the extension of the SO {9, 1) Lorentz symmetry. We note that this symmetry is 
possible only when we define the theory in terms of the supertrace, which is defined as strM — J2p=i ~ 
M33. The cyclic rule of the supermatrix holds true not for the ordinary trace but for the supertrace. This 
property can be verified by the following argument. We consider the following two arbitrary matrices 

Zi — ( '^^^ ] , Z2 — ( ^Jr ^^'^ 1 1 where An and Bn are 32 x 32 bosonic matrices, aio , 

V "21 0.22 J V 1^21 b22 J 

"21 , /3i2 and /32i are 32-component fermionic vectors, and 022 and &22 are c-numbers. We consider the 
supertrace of .^1.^2 and Z2Z1, which are calculated as str{ZiZ2) = tr{AiiBii+ai2(32i)~ {'^2i(^i2+ci22b22) , 
and str{Z2Zi) = tr{AiiBii + pi2aJi) — {0^iai2 + a22&22)- Since a and (3 are fermionic quantities, the 
sign flips if we change the order. Therefore, we establish that it is not an ordinary trace but a supertrace 
that the cyclic rule str{ZiZ2) = str{Z2Zi) holds true of the supermatrix. 

We next discuss its gauge symmetry. This model is invariant under the following U (N) transformation: 

6M = i[(l33x33 «> u), M], where iu G u{N). (3.16) 

In this symmetry, we do not mix the osp(l|32, R) transformation and the u(N) transformation. We have 
a closer look at the mixed symmetry in an attempt to realize the local Lorentz symmetry, and we call 
such models "gauged models". On the other hand, we call the action H3.15|l "nongauged", because we do 
not mix these two symmetry. This "nongauged symmetry" is close to the symmetry of the IIB matrix 
model, in which the 5*0(9, 1) Lorentz symmetry and the U{N) gauge symmetry are separate. 
In addition, this model incorporates the trivial shift of the supermatrix 

6M = cljvxjv, where c G osp{l\32,R). (3.17) 

The action (|3.15|l is rewritten in terms of the explicit expression of the osp(l|32, R) super Lie algebra 

S = -^{tr{m''m^m'') - 3#'^mV) = \Tr{mp'i[mq' .rUrP]) - ^[m,!^]). (3.18) 
^9 9 

We consider the reduction of the model to the ten-dimensions. To this end, we first rewrite the action 
(I3.15|) in terms of the expansion for the eleven-dimensional gamma matrices as (|3.9f) . In order to see 
the correspondence between the ten-dimensional IIB matrix model, we perform a reduction to the ten 
dimensions. To this end, we specialize the tenth direction x^^ = a;', and we retain the 50(9, 1) rotational 
symmetry. We define the ten-dimensional fields by the following chiral decomposition of the bosonic 
matrix m: 

+ ^ (/;i+^^,r^^-^^(i + r(+)) + /(7.'.^,r^--^=(i + r(-))) . (3.19) 

Here, we also recall the notation of the left(right)-handed fermion: 

V'L = ^ "0, V'i?. = ~~"2~~^' 
The coefficients /^^.-l.^j and J^^-./js self-dual to each other. Namely, they satisfy the relation 

The action is thus expanded as 

Sb = ^TrwxJv(-96[4-[),4;)]0'^i''^ - 96M/[v4(+)^ 4")] + 
9 3 

- AH,,...,A[A'^+Jj<^-^'^---^^-]-[A^-Jj^+^^^^^^ 

Q 

I -ffP (\tM 7-(-)mi-"M51 _ [/-{-) /-(+)mi---P51\ 
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27' 

(3.21) 



3i 
2! 
3i 
4! 

where Sh and 5'/ are the bosonic and fermionic parts of the action S, respectively. The whole action is 
S ~ Sb + Sf. The computation of this action is lengthy, and we delegate the proof to Appendix B.2. 
of [27]. 

3.2.2 Identification of tiie supersymmetry witii tlie IIB matrix model 

We next discuss the identification of the supersymmetry with that of the IIB matrix model. This model 
has two kinds of supersymmetry transformation, corresponding to the osp{l\32, R) rotational symmetry 
and the translation symmetry. 

• The former supersymmetry corresponding to the osp{l\32, R) rotation is called the "homogeneous 
supersymmetry" , in the same sense as for the IIB matrix model. This is defined by the rotation by 

/ Y \ 

the supercharge Q = ( . _ q j • Namely, the transformation of the matrices M is given by 

4"« = IO,M, = l( ° J ) . ( * )l = ( "l^^' -7 ) . (3.23) 

• The latter supersymmetry corresponds to the translation symmetry. This is called the " inhomoge- 

(2) 

neons supersymmetry", and is defined by the simple translation of the fermion 5^ ?/> = ^. 

In this way, the osp(l|32, R) supermatrix model has 32 + 32 = 64 supercharges. This is twice as much as 
the IIB matrix model. This simply leads us to speculate that the double structures of the IIB matrix model 
are embedded in this supermatrix model. We elaborate on this point, by scrutinizing the commutation 
relation. Especially, we elucidate that there is a nice correspondence with the supersymmetry of the IIB 
matrix model, with the chiral decomposition (|3.19|l . Namely, we identify A^/^^ with the vector fields of 
the IIB matrix model, and consider its supersymmetry transformation. 
Since the fields A^^^ are given by'^ 

= T ^trimV^t), (3.24) 

the homogeneous transformation of the vector fields Ajj^-* and ip is computed as 

4'^4+' = ^^K('5(^^"^)r,J + ^tr((<5Wm)r^j) = ^Xi?r^VK, (3.25) 

4'^4"^ = ^iK('5i'^"^)rj-^tr((<5Wm)r^,) = ^XLr^VL, (3.26) 
= -mx- (3.27) 

The transformation of the vector fields AJ^^"* is beautifully reminiscent of the homogeneous transformation 
of the IIB matrix model. In this sense, we see a correspondence between the chirality of the vector fields 
and the fermions as 

4+)^^fl, 4-)<=»Vl. (3.28) 



^We recall that ^trF^F'* = 1 while ^trFAB^^^ = —1- Here, we do not take a summation with respect to the 
duplicate indices A, B. 
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The chirally-decomposed fermions 'ipL,R clearly have 16 components, which coincides the number of 
the supercharges for the IIB matrix model. In this sense, this identification augurs very well for the 
reproduction of the supersymmetry of the IIB matrix model. On the other hand, the homogeneous 
transformation of the fermion is given by H3.26|l . and is devoid of the commutator of the vector fields like 
the IIB matrix model. This is a downside in the identification of the supersymmetry transformation with 
that of the IIB matrix model. 

Nextly, we have a close look at the commutation relation of the supersymmetry transformation. We 
wish to find the similar structure of the commutation relations as that of the IIB matrix model. We 
consider the following commutation relations: 

(2)[<5f (3)[<5W,<5f]4±), (3.29) 

It is trivial that the commutator (2) vanishes, because the inhomogeneous supersymmetry is nothing but 
a translation of the fermion. Here, we scrutinize the rest of the commutators one by one. 



Supersymmetry transformation {3)[Sx\ St'^^]A'"i^^ 

In the IIB matrix model, only the commutator (|2.33|) remains nonvanishing. We attempt to ex- 
tract the similar structure to the commutation relation H2.33|l . The commutator of the supersymmetry 
transformation for each m and ip is given by 

[S['\si'^]m = -^{xe-ex), , J^'^lV = 0, (3.30) 

which can be derived in the same way as (|2.33|) . Namely, we take the difference of the following paths: 



5(2) _ _ 5^1) _ _ 5(2) _ _ 

m m ^ m + z(x'0 ~ V'x): whereas m ^ m + j(xV' ~ V'x) ^ m + ixii^ + £) ~ *(V' + e)Xj 

5(2) 5^1) ^U) ^(2) 

ip ij) + e —f ip + e — TOx, whereas ip ^ — mx ^ ip + e ~ mx- 
The commutator of the supersymmetry transformation for each A^^^ can be extracted as 

[(5W,<5(2)]4+) = lfr-(([<5«,<5(^)]m)r^) + :±tr{{[5^^\5^^^]m)T = '-e^^T ,XR. (3.31) 
[<5W,<5f 14-) = lir(([<5W,<5f ]m)r^) - :±tr{{[5^^\5f ^m)T = '-IlT.xl- (3.32) 
Extracting the specific chirality of the SUSY parameters, we obtain the nonvanishing commutators as 

- r^T^XL. [5^lS?MV = 0. (3.33) 
The commutators with different chirality of the two supersymmetry parameters clearly vanish: 

[s^ieMf' = [s^ieM^^ (3-34) 

These combinations of the supersymmetry transformations clearly resemble the structure of the IIB 
matrix model. The correspondence between the vector fields and the chirality of the fermionic fields 
matches that obtained by the correspondence of the homogeneous transformation itself. 



Supersymmetry transformation (1)[5;^^\ i5e^'']4^'' 

In the case of the IIB matrix model, this supersymmetry transformation vanishes on shell up to the 
gauge transformation, as we have seen in H2.31|l . We want to find the same vanishing in our cubic model, 
if we are to identify the supersymmetry transformation. However, it turns out that this commutator does 
not vanish. In investigating this commutation relation, it is easier to utilize the following identity: 

[jw,,5W]Af = [g^,[g„Af]]- [g„[g;,,M]] = [[q^,q,],m]. (3.35) 

This commutator of the SUSY transformation is thus obtained by 
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This relation reveals that the commutators of the two homogeneous transformation are 

[<5W,,5(i)]TO = i[(xe-ex),"i], (3.37) 
[<5«,(5(i)]V = »(xe-ex)V'. (3.38) 

We have seen in H3.33|l which chiralities of the supersymmetry parameters correspond to the vector 

fields A^i^\ According to this correspondence, we examine the commutation relation with respect to the 

following cases. 

We first investigate the commutation relation [(5^^^\ (5c]:'']. We examine this commutation relation of 
this SUSY transformation: 

= ^tmSi'l5i'^]m)r,) -Ltr{{[S['lsil^m)T,i) = ^(^^[m, Tje,.), (3.39) 

Here, we recall that the commutator [m, F^] belong to the sp(32) Lie algebra (namely, only the rank- 
1,2,5 contributions remain). Unfortunately, this commutation relation does not vanish exactly in the 
ten-dimensional reduction. Now, we have a close look at which components survive or vanish: Since 
the chirality of the supersymmetry parameter is identical, the contribution of the odd-rank fields of m 
(with respect to the ten-dimensional indices) clearly vanishes. On the other hand, the contribution of 
the even-rank fields VF, C^^^^j ; ^pi -m4 survives. For example, the contribution of the rank-0 field W is 
obtained as 

XL[M/F„,F^]ei = XLWT^T^eLi^ 0). (3.40) 

We next investigate the commutation relation of other chirality. [^XR^'^efl]- The corresponding field is 
now aIi^-* , and the commutation relation is now 

[4'i'C^]4^^ = ^H{[Si%S^,]^]m)T,) + -ltr-(([5W,5«]m)F,j) = l(;^^[m, Fje«), (3.41) 

by completely computing in the same fashion as in (|3.39() . and we note that the contribution of the 
even-rank fields W, C^i/xaj H^-^...^^ remains nonvanishing. 

Even the worse news comes in the commutator for the different chirality of the supersymmetry pa- 
rameter. The commutator for the rank-1 vector field is given by 

[4'.^C]4+^ = ^ir-((['5(\,\CV)r.) + ^tr(([5W,5W]m)F,,) = -L(^^^r,e^ + e^F^mXi), 

(3.42) 

= ^^K([4'i'€V)r,) - ^tr(([<5W,,5(i)]m)F,„) = -l-{i^mT ,xl + XlT ,men). 

(3.43) 

The disastrous fact is that these commutators include the fields A^^^ and as 

[<5W,<5W]4+) =-^XL4+^r/e^, + ---, [<5«,CH"^ --^XL4-^r/.«^ + .... (3.44) 

These commutation relations reveal that the two-fold supersymmetry is not independent of each other, 
but are connected by not the impurity W, ^'^'^ ^^mi -M4' fields A\^^ . This is an unfavorable 

situation in the analysis of the supersymmetry transformation of this cubic model. 

Summary of the supersymmetry structure 

We recapitulate what we have obtained from the above argument of the supersymmetry. Firstly, 
we have noted that this supermatrix model incorporates the two-fold structure of the supersymmetry 
of the IIB matrix model. This is predictable from the number of the supercharges. This model has 
the homogeneous supersymmetry and the inhomogeneous supersymmetry coming from the osp(l|32,i?) 
rotation and the translation symmetry, respectively. Therefore, this model has 64 supercharges, which is 
twice as many as the IIB matrix model. The IIB-like structure of the supersymmetry is extracted by the 
chiral decomposition of the rank-1 field and the supercharges, and we have seen the correspondence (|3.28|) . 
While we have reproduced the commutation relation (|2.31() and H2.32|l from the commutator [5x'^ , 5^'^] and 
[(5^'^\ Je^''], the commutator [J^^"* , <5e^^] does not trivially vanish, unlike the IIB matrix model. Moreover, 
the commutation relation (|3.44|l unravels that these two supersymmetries are not independent of each 
other. Its physical interpretation will be reported elsewhere. 



29 



3.2.3 Induction of the IIB matrix model 



We have seen the correspondence of the supersymmetry with the IIB matrix model, and we have elucidated 
the correspondence of the chirality between the rank-1 vector fields and the supersymmetry parameter as 
This leads us to the idea that we may be able to induce the IIB matrix model if we integrate out 



certain fields. To this end, we retain the field A^^'' and -0^;, and integrate out the other fields. However, 
the solid discussion for the extraction of the IIB matrix model is an onerous question, and we contend 
ourselves with a hand- waving argument. This supermatrix model possesses no time spacetime derivative 
from the outset, as is true of the IIB matrix model. Its classical equation of motion Oa g(gjfx) ~ ^ ~ ^ 
clearly incorporates the following noncommutative background 

(4^^)b.g. =Pa.' (other fields) = 0, (3.45) 

where satisfies the canonical commutation relation [p^j,,Pi,\ — ic^y. When we denote the fiuctuation 
as yl^^' = j5p + a|i^\ the commutator is mapped as 

Then, the bosonic and the fermionic terms of the action are expanded as 



Sr, = 



4TrArxjv(-96(a^,4;))C^i^^ + 96(a^W^)A(-)^ + 4(a^,iJ^,...^J/(-)''--^^) (3.46) 



9' 



8 
3^ 

1 

27' 

(3.47) 

Sf = ^Tr^^xJv(-3^^/;L^^a^0L) (3.48) 
+ ^TrNy.N{-i^{~i'L[W,^B]+i'B[W,^L])-Hi'LT^[a^+\^L]+4^ 



3i 
2! 

- |(-V;Lr^^^^^^'^^[i?,.,M2P3M4,'/'i?] +^flr''^'^^'"'"[^MlP2;.3M4,^L]) 

- ^_{24,Lr^^''''''^^'-'''^[I^+^,,,,,,^^^^^^ (3.49) 

The cubic action, per se, does not include a quartic term of vector fields in the action. However, we can 
interpret that the bosonic term is induced by the fermionic term of the the IIB matrix model. The idea 
that the theory consisting of fermionic fields and a Dirac operator induces the Einstein gravity and the 
Yang-Mills theory has long been suggested. The proposal of Connes and Chamseddine is one of these 
suggestions of the induced gravity [3]. Based on these suggestions, we find it a natural idea that the 
bosonic term of the IIB matrix model is induced from its fermionic term. And we hypothesize that the 
whole IIB model should be induced only by the fermionic field. 

Our goal is thus to find the fermionic terms in this cubic matrix model to be identified with that of 
the IIB matrix model. We have seen in the previous section the correspondence between the vector fields 
and the fermionic fields according to the identification of Af = 2 supersymmetry with that of the IIB 
matrix model. And the fermionic terms to be identified with that if the IIB matrix model is 

i^RT'^Al+^^jR<^i,Lr^A[-^^jL. (3.50) 

If either of these terms exists in the action, this can be identified with the fermionic term of the IIB matrix 
model. However, the kinetic term of the fermionic fields (|3.48|) . per se, does not include such terms as 
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(|H.5()(I , and does not serve to induce the IIB matrix model due to the discrepancy of the correspondence of 
the vector fields and the chirality of the fermions. In order to remedy this situation, we consider inducing 
the terms (|3.50|) by the multi-loop effect. The action (|3.4t)|) ^ (|3.49|l just tells us that there is no terms to 
be identified with the IIB matrix model at a tree level. Since the idea of 'inducing the IIB matrix model' 
is to construct the bosonic term by the one-loop effect of the fermionic Lagrangian, the idea of 'induced 
theory' is a notion based on the multi-loop effect of the perturbative theory. In this sense, there is no 
problem if the fermionic term is induced by the multi-loop effect. 

We can read off the Feynman rule from the above action. The typical propagators and the vertices are 
given below: Our goal is to build a fermion vertex (|3.50l) by means of the multi-loop effect. In order to 

Propagators 

Figure 3: Propagators of this cubic supermatrix model. 
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Figure 4: Vertices of this cubic supermatrix model. 

build such a term, the above Feynman diagrams do not suffice. Especially we are lacking the propagators 
of the fields. It is thus necessary to build an induced propagators by means of the multi-loop effect of 
the existing Feynman rule. These induced propagators are constructed by the following loop effect. 

Miscellaneous induced Propagators 
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Figure 5: The induced propagators of the bosonic fields 



• {all al ) : It is easy to construct this induced propagator. All we have to do is to connect the 
fermions using the existing fermion propagators (V'iV'i)- 

• (WW) and {A^i Ai ') : The construction of these propagators is a difficult problem. Unfortunately, 
it is impossible to construct these propagators perturbatively, and we delegate the disproof of their 
existence in Appendix B.4. of [27]. However, this is not the end of the story. Even if we fail to 
induce these propagators perturbatively, we have a choice to induce these propagators by means of 
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the nonperturbative effect. These propagators may be induced by the following recursive structure. 
This structure is reminiscent of the self-consistency condition of Nambu-Jona-Lasino model. And 
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Figure 6: The induced propagators of the bosonic fields 

when we consider the nonperturbative^ effect, there is no particular conservation law which prohibits 
the existence of the propagators (VFVF) or A^"''). Throughout our discussion, we assume the 
existence of these propagators. 

Now we are ready to construct the induced propagators of this cubic model. In constructing the vertices 
of the fermionic terms (|3.50|l . The answer is now easy, and the IIB-like vertex is constructed by the 
following procedure. 
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Figure 7: The induced IKKT-like vertex ^jT'^Ay'ip 



• To construct the induced vertex corresponding to (|3.5()|l . we must first construct two objects. One 
is the propagator (iPr^^r). This is easily constructed once we admit the existence of the propagator 
(WW). 

• Another object is the vertex V^^jF^^a^^Vi?- It is also easy to construct this vertex utilizing the 
induced propagator (W^VF). 

• These two objects serve to induce our desired term ?/'flF''yl^^''?/'/j. Note that this induced propagator 
indicates that the kinetic term ipptV^df^^R. This is diagrammatically regarded as the vertex of V'fl, 
il^R and ^ P^,, where p^, around which we have expanded the theory. Therefore, the sum of 
these two objects is regarded as 

(Mr) + ^BT^^al+^R = -i^PRd^^R + V^flF^a^^+Vfl, = V5i?.F^4+) Vfl- (3.51) 

We thus build the vertex for the bosonic field a[^'' and the corresponding fermionic field ipR- This 
fermionic operator serves to derive the bosonic quartic terms of the IIB matrix model. In this way, we 
expect that the IIB matrix model is induced from the osp(l|32, R) supermatrix model. 

*Here, we mean the word 'nonperturbative' by the effect not stemming from the multi-loop effect of Feynman diagram. 
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3.3 gl{l\32,R) <S)gl{N,R) gauged cubic matrix model 

Wc next investigate the gauged version of the supermatrix modeL OriginaUy, L. Smohn [19] touched 
on the suggestion to enhance the symmetry of the matrix model, by regarding the tensor product as 
not for the Lie group but for the Lie algebra. This attempt is fascinating, in that it broadens the 
symmetry greatly. Smolin's proposal turns out to be essential to build a matrix model with local Lorentz 
invariance. To elaborate on this point, let us think about the extension of the IIB matrix model to the 
local Lorentz invariant model. In the IIB matrix model, we regard the eigenvalues of the matrices as the 
spacctime coordinate. Therefore, in order to render the Lorentz symmetry "local", it is imperative that 
the parameter of the Lorentz transformation should be dependent on the u{N) element in a nontrivial 
way [26,38]. Then, we introduce the so-called "gauged matrix model", in order to build such a matrix 
model. 

In the IIB matrix model, the SO{9, 1) Lorentz symmetry and the SU {N) symmetry are totally 
decoupled. Now, let ^ and u be the generator of the 50(9, 1) and the U{N) symmetry, respectively (i.e. 
£, e so(9, 1) and u G u{N)). The symmetry is a tensor product of the group; namely they are decoupled 
in terms of the Lie algebra as 

exp(^ (g) 1 + 1 (g) u) = «) e". (3.52) 

This holds of the nongauged version of the supermatrix model. 

On the other hand, in the gauged model, we enhance the symmetry drastically by mixing these two 
symmetry. Namely, the element of the transformation group is not exp(^ 1 + 1 (g) u) but 

exp(^(g)«). (3.53) 

In this way, we build a Lorentz symmetry dependent on the u{N) symmetry. We define a generator with 

the tensor product of the two Lie algebra. However, we need a caution in defining a tensor product of 
the Lie algebra. Suppose A,B are two different Lie algebra whose bases are {oj} and {bi}, respectively. 
Generally, the space B, which is spanned by the basis a, (g) bj, does not close with respect to the 
commutator 

[ai (g bi, 02 (g 62] = ^ ([cii, 02] {^1, ''2} + {ai, 02} <Si [&i, ^2]) , for ai, 02 G A, 61, 62 e B. (3.54) 

Rather, we define the tensor product by A^B, which is the smallest Lie algebra that includes A^ B as 

a subset'''. 

In the large-A^ reduced models, there are several ways to interpret the spacetime coordinate. In the 
twisted reduced models [65,66], the matrices represent the covariant derivative in the spacetime. On 

the other hand, in the IIB matrix model, the eigenvalues of the matrices arc regarded as the spacetime 
coordinate. Both relations are linked by the expansion around the flat noncommutative background 

Af,=Pf, + aij,, where [Pf,,Pu\ = ic^u- (3.55) 

The IIB matrix model reduces to the non-commutative Yang-Mills theory after this reduction. The 
fermionic term of the IIB matrix model ^Tr^r^[A^, ^] reduced to the action of the fermion in the flat 
space 

I d'^x^{x)i r^id^^{x) + [af,{x), ^Pix)]). (3.56) 

^To illustrate this idea, we introduce the following simple case. The Lie algebra su(6) is known as the tensor product of 
the Lie algebras su(3) and su{2): 

s«(6) = s«(3)(8)s«(2). 

This fact is discerned as follows. Let A" and ct' be the basis of the Lie algebra su(3) and su{2), respectively. su{3) and 
sm(2) arc 8 and 3 dimensional Lie algebras respectively, and the indices run a = 1, • • ■ , 8 and i = 1,2, 3. The tensor product 

su{3)^su{2) consists of the following elements. 

• A" ig) cr*: These are the elements of the tensor product as a set of matrices: sti(3) su(2). This set, per se, does not 
constitute a closed Lie algebra. 

• A" (g) 1 and 1 igi cr': These are the generators of the group SU{3) X SU{2). 

These elements are known to constitute the algebra of sii(6). This example shows in a pedagogical way how the notion of 
'gauged theory' enhances the gauge symmetry. While the Lie algebra of the gauge group 5(7(3) X 5(7(2) is a 8 + 3 = 11 
dimensional algebra, the tensor product S'u(3)cg)s«(2) is a 8 + 3 + 8 X 3 = 35 dimensional Lie algebra. This is the structure 
of the enhancement of the gauge symmetry. Note that the similar enhancement of the gauge symmetry is seen in Smolin's 
proposal [19,23]. 
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This reminds us of the correspondence between the differential operators and the matrices. 

With this idea in mind, we attempt to formulate a matrix model with the local Lorentz invariance, in 
which the matrices look like differential operators. The action should reduce in the classical low-energy 
limit to the fermionic action on the curved spacetime 

Sf= J d''xe{x)i;ix)irf'e^'{x) (^d,^{x) + + jr''Pu;^,p{x)^P{x)^ . (3.57) 

Here, the indices ^, z^, p, • • • and i, j, fc, • • • both run over 0, 1, • ■ • , d — 1. The former and the latter denote 
the indices of the Minkowskian and the curved spacetime, respectively. In considering the correspondence 
with the matrix model, we need to absorb e{x) into the definition of the fermionic field, since it is not 
/ d'^xe{x) but / d'^x that corresponds to the trace of the large N matrices. And we regard the "spinor 
root density" 

^{x) = e^x)^{x) (3.58) 
as the fundamental quantity. Then, the action is rewritten as 

i'^x^{x)e^x)ir''ep'{x) \^d,{e-^ {x)'^{x)) + [A,{x),e-^x)'^{x)] 



+ ^r''''ij,,p{x)e-i{x)-i>{x) 



d'^xl ■^ix)iT'' 



e^,\x)d, + iep'(a;)w,p^(x) + e^,\x)e^^ {x){die ^ {x)) 



+z'l'(x)r'^epX^)[^^(^),*(a:)]+^^'(x)P'^'^^^^e[^/(x)c.,^,^3](a;)vI/(a;) 

d''xi^^{x)iT^^e^\x){d,^{x) + [Mx),^{x)]) + ^'I'(x)r^i^^^^e[^/(x)c.,^,^3](x)vE'(x) 

(3.59) 

where we have utilized in the last equality the fact that the fermionic field '^{x) is Majorana, which leads 
to the cancellation ^{x)T^^'^{x) — 0. The corresponding matrix model is formulated as 

S'f = -\Tri:V^^[A^,i,] - '-Tri>T^^^>'-^-{A^,^,^,,^l^}. (3.60) 

In promoting the action (|3.59|) to the matrix model, we have identified the covariant derivative with the 
commutator with the rank-1 matrix A^. The rank-3 term is a naive product of the spin connection and 
the fermion, and it is natural to promote the product to the anticommutator of the large N hermitian 
matrices^. In the rank-3 term, the pure imaginary number i is necessary so that the action should be 
hermitian. The action (jH.fiOjl is an addition of the rank-3 term to the fermionic term of the action of IIB 
matrix model. Especially when ?/) is a Majorana fermion, l|3.60|l is equivalent to the following action: 

S'^ = -rr^(r^Ap + zr^i^=^M^,p,^3)V. (3.61) 

The equivalence between (|3.6UII and H3.61|l is verified as 



S'f = --Tri>T^^[A^,i^]--Tri,V^^-^-^^-{Ap,p,^,,4,} 



°For the readers' convenience, we summarize the commutation relations of the hermitian and anti-hermitian operators. 
Let H and A be the set of the hermitian and anti-hermitian operators, respectively. When h,h\,h2 G H and a, ai , a2 G A, 
their commutation relations are as follows: 

[/ii,h2]GA, [ft,,a]eH, [ai,a2]GA, {/ii,/i2}6H, {/i,a}eA, {01,02} GH. 

We prove only the first relation, because the others are easily derived in the same way: 

[hi, ^2]^ = {h\h2 — /i2hi)^ = h\h\ — h\h\ = h2h\ — h\h2 = — [hi, /12]. 
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Let us next consider the local Lorentz transformation of this matrix model. Originally the local 
Lorentz transformation of the fermionic field 5ip{x) = ir^^^^e^j^^ (2^)^(2^) should be promoted to the 
anticommutator of the hermitian matrices as 

5^= ir^^^HeMiM2>}, (3.62) 

in order to retain the hermiticity. However, it is a very onerous problem to find an action invariant 
under this transformation. The conundrum stems from the noncommutativity of the matrices. The local 
Lorentz transformation of the fermion H3.62|l with respect to the action (|3.fil|l entails the terms 

= -^Tr^rf'A^r''^'''^Pe,,,, + ■■■. (3.63) 

It is extremely difficult to absorb this transformation via the local Lorentz transformation of the bosonic 
fields Afj,. 

This obstacle forces us to abstain the hermiticity of the matrices. We instead promote this transfor- 
mation to the matrix version as 

S^P = ip^'iA^^e^^^^TA, (3.64) 

and we take the action (|3.61() . Since the fermion is no longer hermitian, the actions (|3.60|l and H3.61|l 
are no longer equivalent. There are two prices we must pay for this alteration. One is that the product 
does not directly correspond to the covariant derivative. The other is that the matrix model H3.61|l 
is no longer invariant under the translation SA^^ = c^l, which is a spacetime translation in the original 
interpretation of the IIB matrix model. 

The local Lorentz transformation of the action H3.61|l is 

6S'^ = ^Tr^ [F^A^ + iT^^'^^'^'A^,^,,,, , F'^^^^e,,,,] V- (3.65) 

Since we are now considering the local Lorentz transformation and their space-time dependent param- 
eters are promoted to u{N) matrices, the infinitesimal parameters are u(N) matrices. Then, the 
commutator 

L*-L ^^1^2^13 I ~ nV^ 1^ J l^MlM2A'3 ! ^1^11^2 J r, 1- ' J L^Ml M2M3 7 ^ I'l 1^2 J 

^ ^ ' ^ ^ ' 

rank- 3 rank-1, 5 

(3.66) 

must include the rank-5 gamma matrices, and the action (|3.61|) is not invariant under the local Lorentz 
transformation. 

The similar thing holds true of the generator of the local Lorentz transformation H3.64|l . In order for 
the algebra to close, only the rank-2 terms are not sufficient, because the commutator 

r^'^^e^,,„F''^''^4^,j = l [rMiA.2^r-^--] {£^,^„4,j + l {r^iM2^r-^--} [£^,,„4,j (3.67) 

rank-2 rank-0, 4 

includes the rank-4 gamma matrices. 

Let us recapitulate what we have obtained through this observation. 

• Since the information of the spacetime is encoded in the bosonic matrices , the parameter of the 
local Lorentz transformation must depend on the u{N^ element in a nontrivial manner. This leads 
us to consider the "gauged" action. 
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• The correspondence of the matrices Ap, and the covariant derivative motivates us to construct a 
matrix model in which the matrices looks like a differential operators. 

• It is difficult to retain the hermiticity of the fermionic field, when we build a local Lorentz invariant 
action. This forces us to abstain the hermiticity of the fields, despite the following drawbacks. 

T^r The term ipT^A^ip no longer corresponds to the covariant derivative ip{x){dfiijj{x) + [ap(x), 'ip{x)]). 

7k- The matrix model H3.61|l is no longer invariant under the translation SA^^ — Cfj,l, which is a 
spacetime translation in the original interpretation of the IIB matrix model 

• The "gauged" symmetry necessitates the higher-rank fields and the local Lorentz transformation 
parameters. 

This observation leads us to regard the " gauged" version of the supermatrix model as a mesmerizing 
candidate for the local Lorentz invariant extension of the IIB matrix model. To this end, we consider the 
gl{l\32,R) super Lie algebra, which is the analytic continuation of the u(l|16, 16), the complexification 
of the osp{l\32, R). In addition, we mix the gl{l\32, R) symmetry and the gl{N) gauge symmetry. 

This observation is also elaborated in another angle without the supermatrix models in Section^ 



3.3.1 Definition of the u(l|16,16) and gl{l\32,R) super Lie algebra 

The gl{l\32,R) super Lie algebra is a cousin of the m(1|16, 16) super Lie algebra. More accurately, 
gl{l\32, R) is obtained by the analytic continuation of the imaginary part of the u(l|16, 16). Therefore, 
we start with introducing the m(1|16, 16) super Lie algebra. 
The very definition of this super Lie algebra is that 

If MGu(l|16,16), thenM^G' + GM = Ofor G= ( ^ "j. (3.68) 

That this is a complexification of osp{l\32, R) can be seen from the following aspect. Unlike the 
osp{l\32, R) super Lie algebra, we do not restrict M to be a real supermatrix, where the reality of 
the supermatrix is defined as M — M* — Therefore, we must replace the transpose by the 

hermitian conjugate. Note that in the real version the complex conjugate is equivalent to the transpose 
according to the property in Appendix. IA.2I 

We can confirm that the super Lie algebra u(l|16, 16) actually closes in totally the same fashion as 
in osp{l\32, R). The legitimacy of the metric G also stems from the same logic as in osp{l\32, R). We 
specify the element of u(l|16, 16) group according to the above definition. The result is that 

If M e u(l|16, 16), then M ^ ^ ^ . (3.69) 

• t; is restricted to be a pure imaginary number. 

• Myij andu^j^...^g are real numbers, while m, UA1A2A3 SinduAi - Ai are pure imaginary numbers. 

This result is derived from the very definition of the super Lie algebra u(l|16, 16). The complex conjugate 
of supermatrices is defined in Appendix. IA.2I 



JVPG + CM = 



mt (j) \ ( \ / pO \ / m t/. 



0. 



(3.70) 



Since v + v'^ — H3.7()(l immediately constrains t; to be a pure imaginary number. 



• We first investigate the constraint of the bosonic matrix m. These are decomposed as (|3.12ll . like 
the osp(l|32, R) super Lie algebra. We use the relationship of the gamma matrices rO'^rA'i -MfcpO — 
(-l) "°+"2"°"" r^i-w. Therefore, (uai-aJ* = (-1) tt^i-A, . This is real for k = 1,2,5 

and pure imaginary for k — 0,3,4. 
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• We investigate the relationship between two fermions ip and (f)^ utilizing the result icf) + T^ip = 0: 

.^t = (jrV)t = (-i)v^f r°) = (-i)v^(-r°) = i^p. (3.71) 

We can verify that this is consistent with the condition tp^T^ + i(f)^ = 0. 

We are thus finished with the determination of the elements of w(l|16, 16) super Lie algebra. 

The important property of m(1|16, 16) super Lie algebra is that these can be uniquely decomposed into 
the direct sum of two different representations of osp(l|32, R). We introduce two different representations 
of osp(l|32,i?) super Lie algebra 

UAi,UA-i_A2,UAi-A^,i>h&T^, (3.72) 

4k M = {M=(^ ^f; ) |m. = « + ^UA^A.A^T^^^^^^ + Iua,...aJ^^--^\ 

u,UAiA2AsiUAi-Ai,iipa,iv G (purc imaginary)}. (3.73) 

And let H and A' be the element of H and A' respectively. Because these are real (pure imaginary), 
these elements respectively satisfy W = ^H' and A'^ = —"^A'. Therefore, these elements satisfy the 
following property 

'^HG + GH = 0, for H eH, '^a'g - GA' = 0, for A' eA'. (3.74) 

The set H is, by definition, the osp(l|32, R) super Lie algebra itself. We investigate an important property 

of the subset of these two subalgebras. The commutation and anti-commutation relations are properties 
of grave importance in getting acquainted with their algebras. 

(l)[ffi,^2] eH, i2)[H,A']eA', mA[,A'^]eH, 

{4){H,,H2}€A', {5){H,A'}€n, {Q){A'„ A'^} € A' . (3.75) 
where ff, ffi, i?2 e W and A', A[, A'^ S A'. 

Here, we only give the proof of the first property, because the others can be easily verified likewise. 

^ [Hi, H2]G = ^ HiG — ^ Hi^ H2G = ^ H2(—GHi) — ^ Hi[—GH2) = GH2H1 — GH1H2 

= -G[Hi,H2]. (3.76) 

Utilizing these relations, we can discern that A', as well as H' is the representations of osp{l\32, R) 
and also that the algebra A' is a representation of osp(l|32,_R) super Lie algebra by the commutation 
relation {2)[H,A'] e A' for H G H and A' G A'. This commutation relation states that A' remain in 
the super Lie algebra A' after the infinitesimal translation by the elements H gH. In this sense, we can 
understand that A' is another representation of osp(l|32, R). 

The introduction of these two representations of osp(l|32, R) teaches us the relationship of osp(l|32, R) 
and w(l|16, 16) super Lie algebras. H(= osp{l\32,Rj) is a real part of u(l|16, 16) Lie algebra, while A' 
is its imaginary part. It is clear that the elements of u(l|16, 16) can be uniquely decomposed into the 
direct sum of H and A'. 

m(1|16,16) = We^', (3.77) 

where ® denotes the direct sum of two sets. 

Nextly, we introduce the gl{l\32, R) super Lie algebra, especially paying attention to the relation with 
the m(1|16, 16). The definition of gl{l\32,R) super Lie algebra is, per se, simple: 

4k If M e glil\32,R), thenM=(^'^^ (3.78) 

• m is an element of the Lie algebra ,9? (32, i?), id est, m is allowed to be an arbitrary 32x32 bosonic ma- 
trix. Decomposing this by the gamma matrices, this can be expressed by TO = '^].^q jiUAi---Ak^^^"''^'° , 
where the coefficients UAi---Ak are all real. 
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• ip and cf> are independent fermionic vectors. Each of them possesses 32 components, and the com- 
ponents are fermionic real number. 

• t; is also a real number. 

The definition of gl{l\32,R) states nothing. This definition just states that an arbitrary real 33 x 33 
supermatrix is an eligible member of the super Lie algebra gl{l\32, R). Although this definition does 
not give any restriction to the elements, the correspondence with the complex group u(l|16,16) is an 
interesting aspect of gl{l\32, R) super Lie algebra. Since ■0 and (p are independent fermionic vectors, 
these can be rewritten as 

■0 = + "02, 4) = 1pl-Tp2- (3.79) 

And the bosonic 32 x 32 matrices are separated by to = mi + 7712, where toi is rank-1,2,5 and TO2 is 
rank-0,3,4. Then, we define the sets A as follows. 

U,UAiA2Aa,UAi---Ai,i'2,V G U}. 

The super Lie algebra gl{l\32,R) is clearly the direct sum of these two super Lie algebras 

gl{l\32,R) = n®A. (3.80) 

These two subalgebras are also the two different representations of osp(l|32, R) super Lie algebra. H 
is osp(l|32, R) itself, and the same super Lie algebra as in introduced in u(l|16, 16). On the other hand, 
the subalgebra A is A' — iAJ And the elements of these subalgebras readily satisfy 

^HG + GH = 0, for H eH, ^ AG -GA = 0, for A e A. (3.81) 

And it is clear that these two subalgebras obey totally the same commutation relations as those of H 
and A': 

il)[Hi,H2]en, {2)[H,A]eA, (3)[^i, A2] e 

(4){i/i, i?2} e A, (5){H, A} e H, (6){Ai, ^2} e A, (3.82) 

where H, Hi, H2 E H and A, Ai,A2 E A. The proof is completely the same as that of (|3.75|) . and we do 
not repeat it. The commutation relation [7Y,.4] G A indicates that ^ is a representation of osp(l|32._R) 
super Lie algebra. 

Now, the relationship of the three super Lie algebras osp{l\32, R), u(l|16, 16) and gl{l\32,R) is clear. 
We have seen that both u(l|16, 16) and gl{l\32,R) are represented by the direct sum of two different 
representations of osp(l|32, R): 

u{l\16A6)^n®A\ gl{l\32,R) = H®A. (3.83) 

The relationship between A and A' is 

A' = iA^ IfAeA, then iA e A'. (3.84) 

In this sense, we can regard gl{l\32, R) super Lie algebra as the analytic continuation of u(l|16,16). 
Although we adopt a matrix theory with the gauge symmetry gl{l\32, R) unlike Smolin's original proposal 
[23], we note that gl{l\32,R) is a cousin of the original m(1|16, 16). 

Nextly, we explain why we have to introduce the complexification of the osp(l|32, R) from the outset. 
It turns out that the simple tensor product osp{l\32, R) (g) u{N) does not close with respect to the com- 
mutator in mixing these two symmetries. Using the relation H3.75(l , we obtain the following commutation 
relation: 

[{n ® H), {n ® H)] = {{n, n} ® [h, h]) + {[n, n] ® {h, h}) = {a' ® a) -f (h ® h). (3.85) 

''This means that, ii A G A, then iA e A' . 
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Therefore, the tensor product osp(l|32, R) (g) u{N) = (g) H is not a closed Lie algebra. In order to build 
a closed Lie algebra, we are urged to include A' (E) A. This satisfies the commutation relation 

[{A' (g A), {A' (g A)] = {{A', A'} (g [A, A]) + {[A', A'] {A, A}) = {A' ® A) + (H «) H). (3.86) 

Therefore, we establish the smallest closed Lie algebra that include osp(l|32, R) (g) u{N) as a subset as 

osp(l|32, R)®u{N) = (Ti; ® H) + {A' ® A). (3.87) 

This is analogous to the situation in which we had to introduce the higher-odd-rank fields in the action 
and the higher-even-rank fields as a parameter of the local Lorentz transformation once we mix the 
symmetries, as we have seen in (|3.66|l and (|3.67|) . 

Here, instead of promoting each element to a hermitian matrix, we can make a closed Lie algebra 
by restricting them to real matrices. It is clear that {H + A) ig) gl{N, R) = gl{l\32, R) (g) gl{N, R) forms 
another closed algebra. In this case, we embed the spacetime into real matrices, instead of hermitian 
matrices. 

3.3.2 Action of gl{l\32, R) gl{N, R) supermatrix model 

The next job is to investigate the action of this theory. The basic idea is similar to osp(f |32, R) non-gauged 
cubic matrix model, and we proceed rather quickly. The action is 

^ 33 32 ^ 

S - -TrATxjv ^ ((^Mp'^AfQ«Af;^^')-M33'5MQ«Affl33) = -rrjvxw(5tr33x33M3) 

^ Q,R=1 p=l ^ 

= — Str{M''M''M'')Tr{T''T^T''). (3.88) 

^ a,b,c=l 

M is now a multiplet of (7Z(1|32, i?) super Lie algebra, with each component promoted to the element of 
gl[N, R) Lie algebra. The indices P,Q,R, - ■ ■ runs P, Q, i?, • • • = 1, • • • , 33, while p, g, r, • • • = 1, • • • , 32. 

We have promoted the 33 x 33 matrix M to a large 33A^ x 33A^ matrices. However, the structure of the 
promotion is completely different from osp(l|32,i?) model. The gauge transformation is with respect to 
not the separate 5/(l|32, i?) and gl{N, R), but the tensor product of the Lie algebra gl{l\32, R)iggl{N, R). 
This drastically enhances the gauge symmetry of the theory. The gauge transformation is thus M ^ 
M + [u, M] for an arbitrary element of m G gZ(l|32, R) (g) gl{N, R). 

We can express this supermatrix model in terms of the basis of the gl{N, R). Here, Mp is expanded 
as = Then, the trace is written as 

TriT-T'T-) = i^r(^'^[T^ T^]) + i^r(^°{^^ T^}) = ^{fabc + dabc). (3.89) 

Here, the structure constant fabc and dabc are real. Using this decomposition, we can rewrite the action 
as 

S = -^{fabc + dabc)Str{M''M'>M-). (3.90) 

-ig 

This supermatrix model can be expressed using its explicit expression H3.78|l as 

S = ^Tr(tr(m^) - 3i4)mip - Si^ipv - v^). (3.91) 

g 

3.3.3 Wigner-Inonii contraction and supersymmetry 

In considering the relation to the IIB matrix model, we need to reduce the model to the ten dimensions. 
However, this matrix model has a grave difference from the osp(l|32,i?) nongauged supermatrix. Since 
the action H3.88|l is no longer based on the commutator for the gl{N\R) matrices, it is not invariant under 
the inhomogeneous supersymmetry, namely the translation of the fermion. Nevertheless, the matrix 
model (|3.88|l accommodates 32 -f 32 = 64 supercharges. This leads us to speculate that this model may 
also have the two-fold supersymmetry structure of the IIB-like supersymmetry. 
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Here, we reduce the model to the ten dimensions by the Wigner-Inonii contraction, in order to extract 
the translation symmetry. Just as we perceive the earth as a flat two-dimensional space because the earth 
is much bigger than we are, we consider the physics in apparently 'ten-dimensional' space because of the 
large radius of the hyperboloid. To this end, we add the linear term of Af to the action H3.88|l and start 
from 

S = ^TrN^NStr{Mf) - R^Tr^^NStrMt. (3.92) 
By adding the linear term, the matrix model (|3.92l) incorporates the classical solution®. 

The expansion around this classical solution amounts to the Wigner-Inonii contraction. We separate the 
original matrix between the classical solution and the fluctuation and the classical solution as 

Then, the action is 

/ = iir((m + RT^f) ~ i{^mi} + ^iipv + R${1 + r«)V') - - R^{tr{m + RT^) - v). (3.95) 

We ignore the terms of 0{R'^), because this is just a constant. The terms of 0{R^) vanish, because this 
is a linear term with respect to the fluctuation. Then, the action is expressed as follows: 

5 = Ritrim^V^) ~ - i0(l + V^)^) + \r{m^) - — - i{^m^ + v^^). (3.96) 

3 3 

We divide the bosonic fluctuation m into m = TOg -I- TOo, where me and rrio consist of the even-rank and 
the odd-rank components in terms of the ten dimensions. Namely, they should be explicitly expressed as 

me = Zl + WV^ + \{C^,,,,Y^'^- + i?^,^,r^^^^«) -f l(G^^...^^r^i-^^ + i/^^...^,r^--'^^«),(3.97) 

mo = i(4+)r^(i + r«) + 4-)r^(i-r«)) 

+ i(4|.)..^,r^--^^(i + r«) + 47.)..^,r''--^^(i-r«)). (3.98) 

Later, we further decompose irio as irio — mi^'' + mi ^ according to the (±), namely the plus or minus 
chirality, in the above decomposition. The fcrmions are divided according to its chirality. The action is 
then written as follows. 

5 ^ R{tr{mlr^)-v^ ~2i^Ril;L)+tr{^ml + meml) 

- i[<i)R{me + v)iiL + i'L[me + v)'4)R + ci)Lmoi)L + (t>Rm^ (3.99) 

We integrate out the fields of order 0{R) and consider the effective theory. In order to cope with the 
cubic term tr(m^), we consider the following rescaling: 

mt = RT"^ RT^ + R-^m'^ + R'^m'^, vt^ R + v = R + R^K', 

i^ = ^pL + yjR = R-H'L + R~'^'R, ^^^L + ^R^ R^4>'l + R'^-^r- (3.100) 



*This classical solution is impossible in the original u(l|16, 16) gauged theory, because the (33,33) component v is 
restricted to be an anti-hermitian matrix. If we are to consider the Wigner-Inonii contraction, one way is to consider the 

quintic action 5„(i|i6,i6) = \Str{Nr^) - R^StrMt. Then, the classical solution (M> = ®^NxN ^ ^ 

is possible because iR ^ IjvxAf is now an anti-hermitian matrix. Another caution is that the gauge group must be not 
m(1|16, 16) (g) su{N,C) but u(l|16, 16) ® u{N,C). If the gauge group is u(l|16, 16) ® su{N,C), the linear term in 
vanishes because the generators are Tr{T°') = 0, and the Wigner Inonii contraction is impossible from the beginning. 
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Following this rcscaling, tr{mf) and (p'^im'^ + v')iI)l are excluded, because this is rescaled as 0{R 2). 
This theory is thus rescaled to be 

S = (trK^r") - v'^ + tr{m',m'^)) - i{2^'n^'^ + ^'^(m', + t^OV'k + V^W^l + (3-101) 

We integrate out the fields m'^, ■^^ and ^'^ by Gaussian integration. Completing this action square, we 
obtain 

S = tr{{m', + ^K^r" + i(V;k^'jr«)}2r«) - \tr{{m'„' + i(Vk'^'J}^r«) 

- («' + '^(Mr)? - li^'Li'Rf - H^'r + I^'l^oMl + \m'M + {{^Wo^'rI (3.102) 

where we have utilized the fact that (F")^ — 132x32, and thus F" possesses an inverse matrix. The effective 
action is given by the following path integral^. However, this effective action turns out to vanish: 



^ W= -hr{{m',^ + ^{l^/n4>L)r^^) + {{^W^i^R) (3-103) 

= -i^rK^F^) + \{tr{-m'^{^'MY^) + ~4>'l<^'r) + \{tr{{^'RAni) - {^'M') = 0- 

That the first term vanishes is discerned from the anti-commutativity of the matrices m'g and F*'; namely 
rngF' = — F^m^. Noting this fact, we rewrite this term as 

_ itr{rn'„^T^) ^^^■'^°>=° hrim'^'T^m'^) ^^='' ifr(m;4F«) = 0. (3.104) 

While the effective action in the ten dimensions turns out to vanish, we investigate its symmetry. 

This model is invariant imdcr the gl{l\'i2,R) rotation 5M = M] where A = ( _ 5 )■ ^'^'^ ^^'^ 

osp(l|32, R) nongauged model, the fermionic part of this rotation serves as the supercharge. In this sense, 
the fifZ(l|32, R) supermatrix model also has 32 + 32 = 64 supercharges. We consider the following rescaling 
of this parameter: 

A-( < + R~~^a'o ^l + R~~^x'r\ /o^n.^ 



For the 3/(1132,7?) rotation 



RT^ + i{x(l) - ^e) -{m + RTi)x + a^-bilj- 



5M = [A, M] = { rprtr ' -liT ' Y-) '1 ' n T ' ) ' (3-106) 

^ ■' y ie{m + itF") — i{(pa + ve) + tb(p i{eip — cpx) + [b,v\ ) ^ ' 

each component is transformed as follows: 

= [ae,m'J + [ao,m'J+z(xL4 + Xfi4)-*(V'L4 + x'/i4), (3-107) 

5m; = K£«] + K,m;] + i(x'i0i-V'fl4), (3-108) 

5Vl = -(mUfl + <Xi) + (aoVfl + ae^L)-^'VL + ^^Yi:> (3-109) 

H'r = 2x^+(a'eV'fl-&V'fl-(m;x'J), (3.110) 

= ~2e^ + ((eWo)+b'4>'L-4>WJ- (3-111) 

= (4m^ + 4m;) + 6% - (0^a:, + 0>^) - t;'e'^, (3.112) 

Sv' = i(e'j,^'^ + e'Lm',)-t{$'j,x'L+^'LXR) + [b',v']. (3.113) 



^The analogy of this path integral is a following toy model. (i)For the action S = ax + bx + c, the path integral is 
= J dx exp(— (ai^ + bx + c)) = J dx exp(— a(x + ^)^) exp(— c + oc exp(— c + f^)- And thus the effective action 
is = 4a^°^ - (ii)The second example is the following action, S = axy + bx + cy, with x, y being auxiliary fields. This 
integration is performed by e"'^ = J dxdy exp(a{x + ^)(2/ + f ) + a^^d thus the effective action is W - 
even if x, y are Grassmann odd quantities. 
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These results possess two major significances. First is that this reconfirms that the effective action ()3.1()3|l 
vanish. The effective action is invariant under the transformation of the fields m'^, tp'j^ and (/)^. Therefore, 
the effective action makes no difference if we translate these fermions arbitrarily. Taking the fields 771^, 
il^'pj and 0^ to be all zero, we find that the effective action is = 0. Therefore, the effective action 
remains W = Q even if we translate these fields into non-zero values. 

The second significance is that this clarifies the structure of the supersymmctry of the effective the- 
ory. This gauged theory originally possesses no trivial translation like the osp(l|32, R) non-gauged model, 
because this theory is deprived of the symmetry of the commutator. However, the Wigner-Inonii con- 
traction serves to retrieve the translation in the ten dimensions. Actually, the underlined transformation 
of H3. 108(1 . H3.110|l and (|3. 111(1 represents a translation independent of the matter field M. We thus find 
it natural to extract the structure of the IIB matrix model from the fields TOo^\ ip'j^ and (pL- Their 
supersymmetry transformation is read off from l(3.108() . 1(3.110(1 and ((3.111(1 as 

= -iV-i^^, Sm't^^M^J'^, S^'j, = 2x'n-{mWL), S^'^ = -2e'L + {e'j.m'J. (3.114) 

This immediately gives the correspondence of the chirality among the bosonic fields and the fermionic 
fields: {mi~^ and {mi'^^ ^^'r)- This pairing also resembles the case of the osp{l\32, R) nongauged 
model. 

We first scrutinize the former correspondence. In this case, the homogeneous and the inhomogeneous 
supersymmetry is identified as 

4i^-[Qx.,-]=[(o \'),'], c'=[Qe.,-]=[(^^^ ;;),.]• (3.115) 

Then, the rank-1 fields are transformed by this homogeneous transformation as 

Si'^^-^ = ^H^XL4''Lr^)-^tr{txL4>'L^^t) = -^4>'Lr,XL. (3.ii6) 

The inhomogeneous translation of course do not affect the transformation of the vector fields. 

We next explore the commutation relations. Firstly, we discern that the commutator [Sxlj^I>l]i 
well as [Si^ ,Sly^], manifestly vanish. 

[«.«'] = [[( 2 "o'-)-(o "o )'••' = "• 

[€'.«] = [[(.?, S),(4 S)l.-1 = 0, (3.118) 

This is a great advantage compared with the osp(l|32, R) case, in which we are beset by the nonvanishing 
commutator of the homogeneous supersymmetry. The commutation relation [(5^^^ i '^ct''] is obtained by 

m^eMV-'^-^L^.XL. (3.119) 
The correspondence of the latter case goes in totally the same way. The supercharges are defined by 

Then, the supersymmetry transformation of the rank-1 fields is obtained by 

Si'^A'^^^ - ^tr{'^i^'R^Br^) + ^tr{-zij'j,enr^t) = ^enr^^'j,. (3.121) 
The commutators [6^]^ , Sr^]] and [S^r, S^^r] vanish likewise, and the commutation relation [Si]^ , S^^] gives 

[C^'5S]4*^^-f^-H^.X«. (3.122) 

We see a better correspondence of the supersymmetry with the IIB matrix model than the osp(l|32, R) 
model. However, the supersymmetry transformation of the fermion is not balanced by the commutator 
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[Afj., A^]. This obstacles are shared with the osp{l\32, R) case. The decisive drawback of this argument is 
that the ten-dimensional effective action totally vanishes. This comes from too large a symmetry of the 
model. However, it may be related to the notion of the topological matrix model [8], in which the IIB 
matrix model is induced ex nihilo. The exploration of such a gauged matrix model are interesting from 
various points of view, especially the existence of the local Lorentz invariance, and it is worth further 
investigations. 

3.4 Curved-space classical solutions of the massive osp{l\32,R) nongauged 
supermatrix model 

In this section, we review the work [46] , in which we elaborated on the relation between the supermatrix 
model and the curved-space background. In Sect ion we have introduced several generalizations of the 
IIB matrix model to accommodate the curved-space background. The work [46] focuses on the similarity 
between the massive IIB matrix model (|2. 14011 and the osp(l[32, R) supermatrix model with the quadratic 
term. This has led us to expect that the osp{l\32, R) supermatrix model could also have interesting 
noncommutative static solutions. We delegate the detailed properties of the higher-dimensional fuzzy 
sphere, which plays a main role in this analysis, to Section Ti. 4. 21 



3.4.1 Action 

Here, we go back to the osp(l[32, R) nongauged supermatrix model instead of the gauged matrix model, 
since our aim is to unravel the similarity with the alteration of the IIB matrix model described in 1)2.140(1 . 
We start with the following massive model with the cubic interaction^*^: 

33 / / 32 \ \ 33 / / 32 \ \ 

S = 3^,Tr E E - MgM^Q' ] - ^Tr E K' 1 " [M§, Mf ] 

(5=1 \ \p=i / / QM=l \ \p=l / / 

= S^iTrimlmP - 2ii>'4j) - iTr{ml[ml,mP] - 3ii^p[ml, Vg])- (3.123) 

Here, P,Q,R,--- and p,q,r, - ■ ■ likewise run over 1, 2, • • • , 33 and 1, 2, • • • , 32. Here, we focus on the 
nongauged action, and this action has the separated osp{l\32, R) rotational symmetry and the U{N) 
gauge symmetry. We reduce the above matrix model to ten dimensions, by specializing the tenth space 
direction x^" = x". In Section f3.2l we focused on the identification of the supersymmetry with that of the 
IIB matrix model, and we resorted to the chiral decomposition. However, we now focus on the nontrivial 
curved-space background, which leads us to reduce the bosonic 32 x 32 matrix m to the ten dimensions 
without the chiral decomposition. Namely, m is reduced as 

m^WTi + A^r^ + B^r^« + Ic^^^^r^i^^ + li/^^...^^r^i-f^» + lz^^...^,r^i-^\ (3.124) 

Then, the relevant action ((3.123|l is rewritten as 



32iTr ^ - 3C^,^, {A^- ,A^-\^ 3C^,^, {B^- , B^^] + 6W'[^^, BP\ + C^,^, {C^\, , C^^^-\ 



^"^We warn the readers that the notation of this section is different from that of [46]. In this paper, we define the path 
integral of the action S as 

where Se is defined in the ten-dimensional Euclidean space by the Wick rotation of the components Aq, ■ ■ • and the gamma 
matrices F", Therefore, the action is different from that of [46] by the overall sign. 



43 



200^ ..^^^j l-l-in/T \H P H 1 

+ 6-ff'''*^^;^p2 [■^A'3A'4A'5PX' ■^MB'-'Piq]) ^ 

+i.^r^i-M4S[^^^ ^^^^] + i.^r'^i-M5[z^^...^^,^]^ . (3.125) 

In the purely cubic supermatrix model (without mass term, which has been studied in [19,23,26,27]), the 
rank-2 field C^i^a possesses a cubic interaction term but has no quadratic term. This has been a severe 
obstacle to the appearance of a Yang-Mills-like structure in the supermatrix model, because it has been 
impossible to identify C^^^a '^ith the commutators of the rank-1 fields [B^^.B^^] (or [A^^^, A^^]). In the 
eleven-dimensional case, this difficulty has been overcome in [37] through the addition of a mass term, and 
we thus expect this model to contain the massive IIB matrix model, the bosonic part of which has been 
studied in [29] to investigate perturbation theory around noncommutative curved-space backgrounds. 

3.4.2 Resolution of the equations of motion 

We proceed to search for possible curved-space classical configurations solving the equations of motion 
that follow from the action l|3.123|) . To get a clearer picture of the problem, we now set the fermions and 
the positive mass-squared bosonic fields to zero: 

^ = W = A^^ H^,...^, = Z^,...^, = 0. (3.126) 

Since their masses are positive (at least in the spatial directions, while the time-like direction of quantum 
fields is generally unphysical), ()3.126|l is a stable classical solution. Furthermore, we choose to identify 
the tachyonic ten-dimensional vector field S^, rather than the well-defined A^^ with the bosonic fields of 
the massive IIB matrix model, in order to obtain a possibly stable curved-space classical solution. The 
classical equations of motion for the remaining tachyonic fields B^ and C^i^ following from (|3. 123(1 are 

Bf, = -i^i~'[B-',C^,,l (3.127) 
C^, = -tfi-\[B^,B,] + [C/,C,p]). (3.128) 

Although it is difficult to solve these equations in full generality, the equation of motion for C^^ suggests 
to take C^i/ oc [Bf^ ,B^] for B^'s satisfying a fairly simple commutator algebra. If wc look for objects 
having a clear geometrical interpretation, it is tempting to look for solutions building fuzzy spheres. 



X X classical solution 

The simplest tentative solution is the product of three fuzzy 2-spheres with the symmetry 5*0(3) x 
S0{3) X 5*0(3). Such a system is described by x iV hermitian matrices building a representation of 
the so(3)(~ su{2)) Lie algebra in the following way: 

- 1 

[B„B,] = ifire,,kBk, Bf + + BI = fi^r^^—lN^N for (ij, k ^ 1,2, 3) (3.129) 

with similar relations for fc = 4, 5, 6 and k = 7,8, 9, trivial commutators for indices that do not 
belong to the same group of 3, and Bq — {eijk is defined as usually). 

This set of fields (|3.129|l describes a space formed by the Cartesian product of three fuzzy spheres 
located in the directions {xi, X2, x^), {xa, x^, xq) and (xy, x^, xg), whose radii are all fir\/ N-^ — 1/2. (A^ — 
l)/4 is the quadratic Casimir operator of the so(3) Lie algebra. Note that any positive-integer value of 
A^ is possible here, since A indexes the dimensions of irreducible representations. For 5*0(3), the irreps 
have dimensions N = 2j + 1, for all integer values of the spin j. However, we can also use spinorial 
representations with half-integer spins in this case. We have to consider this classical solution instead of 
the single 5^ fuzzy sphere 

[B„ Bj] ^ ifire,,kBk (for i, j, fc = 1, 2, 3), = (for ^ = 0, 4, 5, • • • , 9), (3.130) 
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because the solution i34 = • • • = Bg = is unstable in the directions 4 to 9 due to the negative mass- 
squared^^ of the rank-1 fields B^. Without restricting the generality, we can focus on the first sphere 
located in the direction (xi, 2:2, X3), since the three fuzzy spheres all share the same equations of motion. 

In the framework of fuzzy 2-spheres, we can solve the equations of motion H3.127() and (|3.128|l with 
the following ansatz for the rank-2 field C^ : 

C^3 = f{r)e,jkBk, (3.131) 

where f(r) is a function depending on the radius parameter r. Indeed, the equation of motion (|3.128|l 
reduces then to: 

e^jkBki-fir) + r + r,f{r)) = 0. (3.132) 



(|3.132|l has two solutions: /±(r) = 2r ■ When we plug this result in the equation of motion for B 
(|3.127|l . this leads to 

B,(l-2r/±(r)) = 0. (3.133) 
This gives the same condition on the radius parameter r for both f+{r) and f-{r), namely: 



y/l ~ 4r2 = 0. (3.134) 

Therefore, when we assume the ansatz (I3.131|l . we obtain the classical solution l|3.129|l with the radius 
parameter set to r = i, which is fortunately real. Indeed, < would indicate that the fuzzy sphere 
solution is unstable. For example, in the IIB massive matrix model described by (|2.14U|I . the sign of the 
squared radius of the fuzzy 2-sphcre is linked to the sign of the mass term in the action and it would 
become negative for a correct-sign mass term, which is to be expected, since in that case, the trivial 
commutative solution becomes the stable vacuum of the theory. 

We next want to discuss the stability of the x x S'^ classical solution in more qualitative 
terms [25]. To this end, we compare the energy of the trivial commutative solution = with that of 
the fuzzy-sphere solution. The classical energy for = is obviously i?_B^^=o = Sb^=q = 0.^^ 

In the X S'^ X fuzzy-sphere background, the 2-form field Cij is 



Cij — €ijkBk. (3.135) 



Therefore, the total energy is 



9 3 
Es2 = 552 = -64:fi^Tr{B^B'') = -3 x 64:fi^Tr{B,B,) 

fj,— l i—1 

= -l2^i^N{N ^1){N + 1). (3.136) 

This result shows that the S'^ x S"^ x fuzzy-sphere classical solution has a lower energy compared to 
the trivial commutative solution and hence a higher probability. 

Other curved-space solutions and the fuzzy 8-sphere 

So far, we have considered the simplest curved-space solution x x S^. Here, we consider the 
other curved-space solutions. We focus on the fuzzy 8-sphere solution, which has the iS'0(9) rotational 
symmetry and is spanned by the nine-dimcnsional space. We assume the following fuzzy sphere solution 

i?(«)=0, s(«) = ^Gf , (3.137) 



where p runs over the space component 1, 2, • • • , 9, and Gp is already defined by H2. 113(1 . We likewise 
assume the ansatz for the rank-2 field C^^ as 

gg=0, C ^l'> = -^^^-'g{r)[B(,'\B!^% (3.138) 

^^The classical solution with Bq = has no problem, because it has a positive mass unlike the other directions of the 
field B. 

Since we now consider a classical solution with Bq = (thus no need of Wick rotation), the energy is simply equal to 
the classical action in which we substitute the solution. 
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Then, the equation of motion (|3.1'28|l implies 



— [Bf ](-5(r) + 1 + lr^g\r)) = 0. (3.139) 



We again have two choices for the function /(r): 



9±{r)^^j^, . (3.140) 

The equation of motion (|3.127|l for the rank-1 field b'^^ gives 

4'Hl-8rV(r))=0. (3.141) 

Now, imlike the case of the S*^ x S*^ x fuzzy sphere, 1 — 8r^g- (r) = does not have any real positive 

1 

8' 

1 ± - 4(2/c - l)r2 



solution for r. However, there is exactly one such solution for 1 — 8r'^g-i-(r) = 0, which is r — ^ 
More generally, for an S'^'' fuzzy sphere Sp^''-* = {^)G^''\ the same ansatz would give 



9±ir) 



2(2fc- l)r2 



1 — 2kr^g±{r) = , solvable only for g+{r) at r = — . 

We discuss the stability of the fuzzy-sphere classical solution by computing its classical energy. At 
the classical level, we obtain 

5 

Egs = --fi^n{n + 8)7V4. (3.142) 
8 

is given by H2.121|l . and we recaU that this is given by A^4 = {n + l){n + 2){n + 3)'^{n + A)'^{n + 
5)2 (n + 6)(n + 7)/302400. In contrast with the x x case, N can take here only certain precise 
values. For example, the smallest non-trivial representation (n = 1) has dimension 16, the following one 
{n — 2) 126, then 672, etc... The classical energy for the 5*2 x S'^ x fuzzy-sphere solution is of the 
order 0{—^^n^) = 0{—iM^N^) while that of the fuzzy 8-sphere solution is of the order 0{—fi^n^'^) = 
0{— N's). Therefore, at large iV, the 5^ x 5*2 x 52 triple fuzzy-sphere solution is energetically favored 
compared to the solution at the equal size N of the matrices. The presence of a spherical solution for 
all N in the 5*2 x S'2 x 5*2 case may indeed be a stabilizing factor. On the other hand, at equal value of 
n, whose physical meaning is less clear, the fuzzy 8-sphere solution has lower energy. 

We have seen that the fuzzy 8-sphere is energetically favored. This is in contrast to the matrix 
model with the Chern-Simons term (|2.111l) . in which the classical energy is given by l|2.139l) and thus the 
fuzzy 8-sphere is positive. In this case, the fuzzy 8-sphcre is less energetically favored than the trivial 
commutative background. In this sense, this situation has more similarity to the IIB matrix model with 
the tachyonic mass term (I2.140f) . in which the higher-dimensional fuzzy sphere, as well as the 5*2 fuzzy 
sphere, is more energetically favored, as we have seen in (|2.144f) . 

The single fuzzy g-spheres for q < 7 do not constitute a stable classical solution of our model. When 
the S** sphere occupies the direction xi,X2, - ■ ■ ,Xq+i, the solution -6^+2 = = • ■ ■ = Sg*-* = is 

trivially unstable because of the negative mass-squared. Whereas, the Cartesian product of several fuzzy 
spheres, such as S'2 x 5^, is a possible candidate for a stable classical solution. 



3.4.3 Nucleation process of spherical branes 

Starting from a vacuous spacetime, it is interesting to try to guess how spherical brane configurations 
could be successively produced through a sequence of decays into energetically more favorable meta-stable 
brane systems. The reader may have noticed that we have so far limited ourselves to the study of curved 
branes building irreducible representations of their symmetry groups. This could seem at first to be an 
unjustified prejudice, but it turns out that such configurations are energetically favored at equal values 
of N . For example, for 5*0(3), an irreducible representation 7?. at of dimension N contributes as 

En^ = -4^i^iN^ - N) (3.143) 
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per fuzzy 2-sphere, while a reducible representation TZni®- ■ -(BTZn^ of equal dimension A^i + . . ■+Nm = N 
would contribute as 

■m 

En.,e...en.^ = -V^(iVf - N,). (3.144) 

i=l 

This is obviously a less negative number, especially for big values of N. A similar conclusion was reached 
in [25] for the case of a Euclidean three-dimensional IIB matrix model with a Chern-Simons term and 
it seems to be a fairly general feature of matrix models admitting non-trivial classical solutions. This 
property is particularly clear for low-dimensional branes, since the classical energy is of order 0{—^^N^) 
for the 5*^ fuzzy sphere, but it remains true for any fuzzy 2fc-sphere solution, whose energy is of order 
0(— /x^A^^+''/('^('^+^^^), which also shows that low-dimensional configurations are favored. As hinted for 
in the preceding subsection, this latter fact can be physically understood by remarking that there are 
more irreps available for low-dimensional fuzzy spheres, which makes it easier for them to grow in radius 
through energetically favorable configurations. A third obvious fact is that configurations described by 
representations of high dimensionality are preferred. 

Put together, these comparisons give us a possible picture for the brane nucleation process in this and 
similar matrix models. As they appear, configurations of all spacetime dimensions described by small 
representations will be progressively absorbed by bigger representations to form irreducible ones, that 
will slowly grow in this way to bigger values of N. Parallel to that, branes of higher dimensionalities 
will tend to decay into a bunch of branes of smaller dimensionalities, finally leaving only 2-spheres and 
noncommutative tori of growing radii. If the size of the hermitian matrices is left open, as is usually 
the case in completely reduced models, where the path integration contains a sum on that size, no 
configuration will be truly stable, since the size of the irreps will grow continuously. 

Of course, this is a relatively qualitative study, which could only be proven correct by a full quantum 
statistical study of the model. However, it seems to be an interesting proposal for the possible physics of 
such theories. 



3.4.4 Supersymmetry 

We next comment on the structure of the supersymmetry. The biggest difference with the purely cubic 
supermatrix model, due to the addition of the mass term, is that this model is not invariant under the 
inhomogeneous supersymmetry 

(5i„homogc„cousm = 0, (^inhomogcncouaV' = 1 (3.145) 

which is a translation of the fermionic field. However, this model has 2 homogeneous supersymmetries in 
ten dimensions, which are part of the osp{l\32, R) symmetry: 



e \ f m ip 
ie ) '[ iip 



i^eip — ipe) —me 
iem 



(3.146) 



which transforms the bosonic and fermionic fields as 

SeTn — i{e'ip — i^t), Sei^ — —me. (3.147) 

In the IIB matrix model, the supersymmetry has to balance between a quartic term Tr([Afj_, A^])'^ 
and a trilinear contribution Tr^r^[A^, in the action H2.ll) . which implies that the supersymmetry 
transformation of the fermionic field has to be bilinear in the bosonic field. On the other hand, the 
homogeneous supersymmetries are all linear in the fields in the purely cubic supermatrix model [26,27]. 
By incorporating the mass term, we are allowed to integrate out the rank-2 field C^i^ja by solving the 
classical equation of motion iteratively as in [37]^'^. 

^^The Yang-Mills-like structure of the homogeneous supersymmetry transformation on the fermion comes from 
— iCfn^rf"^ t, which is a part of S^tp = —me. The exphcit form of the iterative solution of the equations of motion 
i!j.lij71 and 11^81 is 

Cf,^ = -ifi-^[B^,B^] + ifi-^[[Bi„Bp],[B^,BP]] 

-if,~'[[B^, Bp],[[B,,B^], [B", B>^]]] + i,,-^[[B,,Bp],[[B^,B^],[BP,B^]]] + 
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Thanks to this procedure, the homogeneous SUSY transformation for the fcrmionic field becomes 



(3.148) 



while the transformation of the field i?^ is 




(3.149) 



In that sense, the mass term is essential to realize the Yang-Mills-like structure for the homogeneous 
supersymmetries. On the other hand, if we want to preserve the homogeneous supersymmetries, we cannot 
just put a mass term for C^i/^a by hand. The osp{l\32, R) symmetry forces all fields to share the same 
mass, since they all lie in the same multiplet. In particular, we are forced to introduce a mass term for 
the fermions as well, which breaks the inhomogeneous supersymmetries. In other words, it seems difficult 
to have super- Yang-Mills-type structure for both homogeneous and inhomogeneous supersymmetries in 
the context of supermatrix models. 

Indeed, in contrast with the purely cubic supermatrix model [26,27], which has twice as many su- 
persymmetry parameters, the massive supermatrix model has only M — 2 supersymmetries in ten di- 
mensions, because it lacks the inhomogeneous supersymmetries. In consequence, we cannot realize the 
translation of the vector field as a commutator of two linear combinations of the homogeneous and 
inhomogeneous supersymmetries (|3.147|l as in the IIB matrix model, where it leads to the interpretation 
of the eigenvalues of as spacetime coordinates. On the contrary. 



vanishes up to an sp{32, R) rotation. This problem is a serious obstacle for the identification of the 
supersymmetry of this model with that of the IIB matrix model. More analysis will be reported elsewhere. 

3.5 Other related studies 

We have so far gone over the author's works [26,46] concerning the supermatrix models. We conclude 
this section by touching on other related studies. In [37], Bagnoud, Carlevaro and Bilal discussed the 
osp{l\32, R) supermatrix model in the twelve-dimensional and the eleven-dimensional contexts. Espe- 
cially, they indicated how to perform the infinite momentum frame limit and used the matrix version of 
the T-duality to obtain a supersymmetric matrix quantum mechanics. They introduced the quadratic 
terms as in Sec. 13.41 and resorted to the perturbative technique that led to an infinite tower of the 
higher-order interaction among the physical fields. On the other hand, the lower-order term reproduced 
the BFSS model with an additional mass term together with interaction terms involving the five-brane 
degrees of freedom. 

Other attempts for a background-independent matrix models have been achieved by the cubic action 
of the exceptional Jordan algebra. Firstly, L. Smolin proposed a matrix model based on the exceptional 
Jordan algebra whose automorphism group is F4 (in this sense, we refer to this model as "the F4 model" 
here) . The exceptional Jordan algebra is related to the octonion, and gives a theory that has a compact- 
ification which reduces to the matrix string theory [6] . Its complex extension has been proposed in [35] . 
The Eq exceptional Lie algebra is known to be the automorphism of the complexification. Therefore, the 
Eq matrix model has twice as many degrees of freedom as the F4 model, and derived the effective action 
of the matrix string type. 

4 Matrix model with manifest general coordinate invariance 

In the pervious section, we have investigated various aspects of the supermatrix models based on the 
osp(l|32, R) super Lie algebra. This extension of the IIB matrix model poses a lot of interesting questions. 
Especially, in Sec. 13.31 we have investigated the gauged gl{l\32, R) (g) gl{R) model as a candidate for the 
extension of the IIB matrix model equipped with the local Lorentz invariance. The key idea is to render 
the parameter of the 5*0(9, 1) Lorentz transformation dependent on the u{N) matrices in a nontrivial 
manner. This is an essential alteration from the IIB matrix model, since the spacetime is interpreted 
as being immersed in the eigenvalues of the bosonic matrices A^. In addition, the u(l|16, 16) (and its 



[S^,S^]m = i[{ex - Xe), m], [S^,S^]'iIj = i{ex - xe)V' 



(3.150) 
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analytic connection gl{l\32, R)) supermatrix model incorporates the higher-rank tensor field. Especially, 
the rank-3 field has a possibility to be identified with the spin connection in the quantum field theory on 
the curved spacetime. 

In this section, we review the author's work [38] that inherits the above idea. The most serious obstacle 
of the gl{l\32,R) supermatrix model is that the effective action reduced to the ten dimensions by the 
Wigner-Inonii contraction finally vanishes. This is ascribed to the fact that the element of gl{l\32,R) of 
all ranks are allocated for both the matter fields and the parameter of the local Lorentz transformation. 
These two elements had to be mixed in the eleven-dimensional theory, due to the self-duality condition 
of the gamma matrices (jA.llll . Namely, even when the rank-6 fields come in the local Lorentz parameter 
due to the rule (|3.67|l . this is identified with the rank-5 fields due to the self-duality rule (|A.lip . By the 
same token, the rank-7 matter fields coming from (|3.66(l are tantamount to the rank-4 fields. 

In [38], we address the above approach in a different angle without the supermatrix model. We 
start from the ten-dimensional action from the outset, without the reduction from the eleven-dimensional 
model. Since we start from the ten-dimensional field theory, we can separate the matter fields and the 
local Lorentz parameter. We allocate the odd-rank fields for the matter fields, and the even-rank fields 
go for the local Lorentz parameter. In addition, we explicitly give an expansion of the matrices with 
the differential operators, to embody the idea of the identification of the matrices and the differential 
operators. In these senses, the work [38] gives a clearly different angle from our previous work in Sec. 



4.1 Bosonic part of the matrix model 

This time, we start from the following action, defined in the ten-dimensional spacetime from the outset: 

S = TrNy,N[tr32x32V{m'')]. (4.1) 

The uppercase Tr and the lowercase tr denote the trace for the N x N and 32 x 32 matrices, respectively. 
The indices fj,,!/, - ■ ■ and i,j, ■ ■ ■ both run over 0, 1, • • • , 9. The former is allocated for the flat Minkowskian 
spacetime, and the latter is for the curved spacetime. V{x) is some function determined later, and m 
consists of the odd-rank fields: 

7 1 7 1 

Here, ^^^...^2^-1 ('^ = lj2,3,4, 5) are all hermitian matrices, and this satisfies F^to^^F" — m. In this 
sense, our action is shown to be hermitian as 

^ Tr[trV{{rri'f)] ^ Tr[trV{{T^mT"f)] = Tr[trVim^)] = S. (4.3) 

Here, we exclude the odd-power of rn in the action, because 771^'^+^ all consists of the odd-rank gamma 
matrices and thus couple the fermions of the different chirality. 



4.1.1 Identification of large N matrices with differential operators 

We identify the space of the large- matrices with that of the differential operators. By this identification, 
we can describe the differential operators on an arbitrary spin bundle over an arbitrary manifold in the 
continuum limit simultaneously, because they are embedded in the space of large N matrices, as is 
depicted in Fig. |H1 

To illustrate this idea, we visit several simple examples on the differential operators of the scalar fields 
on two different bundle over the Si circle (see Fig. EJ. We first consider the trivial bundle with the 
periodic condition /(I) = /(O). We discritize the region < cc < 1 into small slices of spacing e = j^. 
Then, the differential operator is approximated by the finite difference as 

1 \ / 1. _ 1 \ \ 

(4.4) 
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Figure 8: The way the spaces of the differential operators are embedded in the space of large N matrices. 
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(1) trivial bundle with the periodic (2) Z2-twisted bundle with the antiperiodic 
condition f(1 )=f(0). condition f(1 )=-f(0). 



Figure 9: Differential operators on (f)the trivial bundle and (2)the Z2-twisted bundle over 5*1. 



Due to the periodic condition, this finite difference is expressed by the large N matrix as 



A = 



N 



/ 
-1 



V 1 



1 

1 
-1 



-1 \ 







(4.5) 



We next consider the similar problem with respect to the Z2-twisted bundle, in which the antiperiodic 

condition /(I) = — /(O) is imposed. Paying attention to this antiperiodicity, we understand that the 
finite difference in the discritized space is expressed as 



2 



/ 
-1 



V -1 



f 

1 
-1 



1 \ 



-1 oj 



(4.6) 



In the following, from the big space of large N matrices, we pick up a subspace consisting of the 
differential operators over one manifold. We regard m as the differential operators over this manifold. 
Then, wc analyze the effective theory of the fields appearing in the expansion of the differential operators 
(the explicit form is given later). 

The space of the differential operator is infinite-dimensional. The trace Tr for this space is generally 
divergent. However, wc can render this trace finite by choosing the function V{'m.'^) damping rapidly. Wc 
would like to identify this matrix m with the Dirac operator in the curved spacetime. Clearly, the Dirac 
operator has a dimensionality [(length)"^], while m is a dimensionless quantity. Then, we introduce a 
constant r with the dimensionality [(length)^], and express m as 



m = T^D, where 



5!^" 



La 

7! 



7 ^ 



3 



(4.7) 
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Here, D and ^pi- -^2ri-i ^-re matrices of dimensionality [(length)"^], r is similar to the Regge slope a' 
in string theory, and is introduced as a reference scale. This parameter r is not a cut-off parameter. 
y(rn^) is an exponentially decreasing function, and the damping factor is supplied by the action itself. 
Since V{m'^) is a function of the dimensionless quantity m, r represents a damping scale. When we 
approximate the differential operators by finite N matrices, an A'^-dependent ultraviolet cut-off naturally 
appears. When we take N to infinity, this cut-off becomes infinitely small. But the scale r is completely 
independent of this ultraviolet cut-off, and takes a constant value even in the large N limit. Thus, we 
can fairly take N to infinity and identify the large N matrices with the differential operators, at least 
when we investigate the effective theory at tree level. In this sense, our model differs from the induced 
gravity. 

The matrices A^j...^2^_j^ are expanded by the number of the derivatives. The hermiticity of A^^...^^^_^ 
leads us to expand it by the anti-commutator of the derivatives as 

°° jfc 

k=l 

Here, the indices ii,---,ik are symmetric, while the indices //i, • • • , /i2fc-i are anti-symmetric. Since 
we identify D with the extension of the Dirac operator in the curved spacetime, we find it natural to 
identify the coefRcients aj^ (x) with the vielbein of the background metric. A simple dimensional analysis 
immediately indicates that the coefficients al]ll'''li'll_-^{x) has a dimensionality [(length)~^+'^]. Then, D is 
identified with the Dirac operator in the curved space as 

e~^{x) + (higher-rank and higher-derivative terms). (4.9) 

4.1.2 Local Lorentz invariance 

We next clarify that our matrix model is invariant under the local Lorentz transformation. Here, the 
local Lorentz transformation is given by 

Sm — [m,e], where 

1 7 1 

i 1 

£ is the parameter of the local Lorentz transformation, and it satisfies F'^e^^F" = e. The coefhcients 
£07 ■ ■ ■ J ^^^l■■■^llo ai'e a-H hcrmitian matrices. The discussion in H3.67() immediately indicates that the closure 
of the algebra of the local Lorentz transformation necessitates all the even-rank fields. The invariance of 
the action is verified as 

5S = Tr[tr{2V' {m^)m[m, e])] = 0. (4.11) 

Here, V'{x) denotes V'{x) = The vanishing of PTT|l can be verified by noting the cyclic rule of 

the trace Tr. Since we are now discussing the space of the infinite-dimensional space, it is not necessarily 
trivial whether the trace satisfies the cyclic rule. However, the cyclic rule holds when the coefficients 
q(«i -'a;)^^ ^^^^^ (^^"^ damps rapidly. The only nontrivial part for the cyclic rule is the commutator between 
the derivative and the fields 

Tr([a„a(j--^'=)^,...^,„_, (.:)])= J d'*x(a;|(a.a(^--^'=)^,..^,„_, (x))|x) = J d'^xid^a'^^^-^"^ ^,...^,„_,{x)){x\x). 

However, this surface term does not contribute when the coefficients a^^^'"^''^ ^-^...^2^_-^{x) damp rapidly. 

4.1.3 Determination of ¥{171^) 

We next determine the function V{in'^). We require that the action H4.1|l should have the ffat metric 
Too = ir^V^dfi as a classical background. However, this turns out to be a draconian constraint on 



D = e^{x) 
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V{x). The bosonic part is expressed by the heat-kernel expansion. We delegate the details of the heat- 
kernel expansion to Appendix ^ We find it more convenient to express V{u) in terms of its Laplace 
transformation 



V{u) = / dsgis)e-'\ (4.12) 





in order to see a more transparent correspondence with the heat kernel. Then, the action is now expressed 
as 

/■oo /• id ( ^ / f'OO \ \ 

TT\tTV{m')\ = dsg{s)Tr\tre-^^^^\ = J ^ ^ ^ ^ dsg{s)s-i+'' j r'^ Ak{x) j , (4.13) 



where Ak{x) are the Seeley-de-Witt coefficients of the heat kernel Tr[tr{e )], namely 



-f^ ^ t'A,{x). (4.14) 



Here, k runs all the integers. We can show that the half-integer k can be excluded in the follow- 
ing way. The coefficients Ak{x) has a dimensionality [(length) ^^'^J. On the other hand the terms 
rife^iC^ii ■ ■ ■ ^i^^o.^^^'"^'''^ iii-ii2r^^-iix)) l^as a dimensionality V = -n + Y.7=i(~Pi + ^0- Therefore, this 
belongs to the coefficients A_2.{x). On the other hand, this term has —n + X]^=l(2'7^^ +Pi + U) indices, 
which is equal to T) modulo 2. Therefore, when T) is odd, these terms cannot contract the indices, and 
thus cannot survive in the action. If such terms existed, they would no longer be a scalar and thus would 
violate the Lorentz invariance of the action. In this way, we eliminate the contribution of the half-integer 
k. 

Then, let us think about the condition for mp to be a classical solution. To this end, the linear 
terms of all the fluctuations must vanish in the action. Now, only the scalar fields can contribute to the 
action, we have only to consider the linear terms a^^'^'i'^'^'^ ' *'*'' (x) (for ^ = 0,1,2,- • •). Since a'^'^\{x) 
is identified with the vielbein, its cancellation means the cancellation of the cosmological constant. This 
amounts to the cancellation of the Seeley-de-Witt coefficient Aq{x). We disregard the derivative terms 
di-^ • • • a(^^ "^^j^...^2^_j^ (x), because these can be eliminated by the partial integration. Since the fields 



a 



(A'*in»2»2 - 2i2i) (3;) Jias a dimensionality [(length)^'], it belongs to the Seeley-de-Witt coefficient A-i{x). 
Therefore, we generally demand that the coefficients A-i{x) for Z = 0, 1, 2, • • • should all cancel. 

However, it turns out that this imposes a demanding condition on the function V{m?). The above 
condition is translated as 



dsg(s)s" = 0, for n = 0, 1, 2, • • • . (4.15) 
This condition can further be rewritten as 

/•oo /"OO /"OO 7 /"OO 

/ duV{u)u"-^+i = du dsg{s)e-'''u''+i-^ =T[n+^) dsg{s)s-i-" = 0, (4.16) 
Jo Jq Jq 2 Jo 

for n = 0, 1, 2, • • •. One example for such a function is the following: 

We can show that the function (|4.17|) satisfies the condition H4.16fl as follows [71]. Firstly, we note the 
following integral 

dyy^'e-'^y = m\a-"'-\ where a = exp(^) = (4.18) 

4 V 2 

This is a real number when to — 3 is a multiple of 4. Taking the imaginary part of the both hand sides, 
we obtain 

/•oo 

dyy^"+' sin(^) exp(--^) = 0, forn = 0, 1, 2, • • • . (4.19) 
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We make a substitution u = ^ and perform the partial integration to obtain the solution H4.17(l "'^^. 

We then derive the Einstein gravity in the low energy limit, from this bosonic part. The linear term of 
the vielbein a^*^^(x) is eliminated, and thus the graviton is now massless since the mass term vanishes due 
to the general coordinate invariance. When we retain the curved-space Dirac operator in the expansion 
(113, the Seeley-de-Witt expansion is now obtained by 

TrMe--')] = / d'^J^ei.) (r^ + • • •) ■ (4.20) 

This term should be clearly allocated for the coefficient Ai (x) . 

We next give a qualitative argument for the mass and the kinetic terms of the matter fields. While 
the explicit heat-kernel expansion of l|4.9|l entails an involved calculation, it is easy to estimate which 
coefficient Ak{x) the mass term and the kinetic term belong to via the dimensional analysis. Since the co- 
efficients a('^'"*'=)^i...^2„-i (^) ^-iid '^''"'^ '''Vi - M2,i-i(^) respectively have the dimensionality [(length)~^+'''] 
and [(length)~^+'], we have 

mass terms:a^*i'"*''Vi---M2„-i(a^)«*^'"'^'Vi---A'2„-i(a;) G A^^k+i{x), (4.2f) 
kinetic terms:(afe,a(*i-'^V,...^,„_,(x))(9fe.a(^'^--"Vi-P2„-i(a;)) e A^_t+i{x). (4.22) 

Here, k+l must be an even number. Firstly, the mass term and the kinetic terms of the fields a^^...^^^_-^ (x) 
belong to Ai{x) and A2{x), respectively. The explicit heat-kernel calculation reveals that the relevant 
mass terms are given by 

rf'^^^^^T i^-a,,...,,{x)a'^^-'^^{x) + -a,,...,,{x)a^--'^^^{x) 

+ ^a^,...^,ix)a'^^-'^^{x) + |a^,...^,(a:)a^--^n^)) + ' ' ' • (4-23) 

We give the derivation of this result in Appendix. ^ Thus, the fields a/2i - /i2n-i (2^) fo'' — 2,3,4,5 
are clearly massive. While (|4.23f) cancels the mass term a^(x)a^(a;), the field a^(x) is also presumably 
massive because there is no reason to prohibit the cross terms a^(a;)a*^*i ' *''^^(a;). 

Next, we consider the fields a'*^ai---a2„-i (2^)- Especially the fields a'^^\ai---a2„i^) (with n — 1,2,3,4) 
are identified with the anti-symmetric tensor fields in the type IIB supergravity. Their mass terms are 
included in ^0(2;), while the kinetic terms reside in Ai{x). These fields are therefore regarded as massless. 
This observation augurs well for the matrix model which reduces to the type IIB supergravity. 

However, it is not clear whether the higher-spin fields a'*^ "*'°-'ai - a2„-i (2;) (fc = 2,3, • • •) are massive, 
since the mass terms and the kinetic terms belong to the coefficients Ai-k{x) and A2-k{x) respectively 
and both of them vanish in the action. 



4.2 Supersymmetric action 

In the previous section, we have discerned that the bosonic part of the matrix model reduces to the 
Einstein gravity in the low energy limit. Our goal is to establish the extension of the IIB matrix model 
that reduces to the type IIB supergravity in the low energy limit. This leads us to scrutinize the fermionic 
part. We pay attention of the correspondence of the supersymmetry. 
We start from the following action: 

Ss = Tr[tr{V{m'^))] + Tri)mip. (4.24) 

Here, m is defined in the same way as in the bosonic model. ip is a, Weyl fermion. However, the crucial 
difference from the IIB matrix model is that we abstain its Majorana property. Namely, we abstain its 

^^This is the famous counter-example against Hausdorff's moment problem when we replace the finite integral with the 
infinite integral. When the function V(u) satisfies 

duu^Viu) = 0, forn = 0, 1, 2, • • ■ , 

The function V{u) must be zero constantly. This problem is called "Hausdorff's moment problem" [71]. However, this is 
not the case with 14.161 because we are considering the infinite integral. 
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hermiticity. It is extremely difficult to build a matrix model that closes with respect to the local Lorentz 
transformation 6ip = jT^'^ {Sf^,^, That is why we abandon the hermiticity and define the local Lorentz 
transformation as 6ip = \T^^'^ e^^ij}. The fermionic field is also expanded by the number of the derivatives 
as 

^= + (e-^^^'r- (4.25) 

Clearly, the fields x^*^ ^'H^^) -^i^^ the dimensionality [(length)']. We put the damping factor (e^'^^^)", 
so that the action {prml> should be finite. We choose this power in accordance with the function (|4.17|) . 
which leads us to set a = |. 

4.2.1 Local Lorentz invariance 

We discuss its local Lorentz invariance, which is defined as 

5m =^ [m,e\, Sip — eip, (4.26) 
where e is defined in H4.10|l . ?/' is transformed as 

= <5(V'^(r°)) = -?/'^(rO)2£tro = -^e. (4.27) 

The action is transformed by this local Lorentz transformation as 

5Ss = 2Tr[tr{Vs{m^))[m,e]]+Tr[tr{^[m,e]ij)] = 0. (4.28) 

This holds true when the coefficients cS^^"'^''\^...^^^_^{x) and x^^^ ^''\^) damp rapidly in the infinite 
distance. 

4.2.2 Determination of Vsim?) 

We next discuss the determination of the function ^5(771^). Here, we do not give a specific form of Vsirn?) 
and we just list up the conditions for Vs{m?)- 

Firstly, the fiat background mo = ir^V^d^ must be a classical solution. As we have seen in the bosonic 
case, this imposes a strict condition on the function Vg{m^); namely this must satisfy the criterion (|4.1t)|) . 

Secondly, this model must retain the even-rank anti-symmetric tensor fields a'^'^i^j^...^^^ (x) (n = 
1,2,3,4). The action Vsim^) must be determined so that these fields should be massless and thus be 
bequeathed in the low-energy limit. 

It is not clear whether the function H4.17|l satisfies the latter condition, since we have yet to elucidate 
the mass of the higher-spin fields. 



4.2.3 Supersymmetry 

We next discuss its supersymmetry. We define the supersymmetric transformation as 

S,m = eip + i:e, 6^^ ^ 2Vsim'^)e. (4.29) 
Here, F^(x) = Now, we assume that the function Vs{u) can be expanded around li = as 

oo 

^«(") = E5"'- (4-30) 

k=l 

Therefore, the following argument does not apply to the function H4.17|l . because it is not possible to 
expand this around the origin u = 0. The supersymmetry transformation of "0 is derived as 

oo 

oo 

= -6tr«(^(-i)^'=-ia2(fc+i)(rVtr«)2'=)r" = 2eV^im'). 

(4.31) 
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Then, the supersymmetry invariance of the action can be easily verified as 

SeSs = Tr[tr{2Vs{m'^)m{eij + tpe) + 2tpmVs{m'^)e + 2emVs{m'^)i))] = 0. (4.32) 
We next analyse the commutation relation of the supersymmetry transformation. This is given by 

[S„6^]m = 2[{^e-e^),V^{m% (4.33) 

/ yi(m2)-yi(o)^ ,n{m^)-v^io) 



with the help of the equation of motion 

^ = 2mil) = 0, ^ = 2^m = 0. (4.35) 

0%) dtp 

This can be verified by taking the difference of the two transformations, as we did for the IIB matrix 
model or the supermatrix model. For the boson, we compare the following two transformation: 

m ^ m + Ci/' + V'C^™+(e + 6V' + '0(e + O + 2eeV"5(m2) + 2y^(TO2)ef, 
m m + ei/^ + Ve^m+(e + 0^ + V'(e + + 2e^^s(m2) + 2T/^(m2)5e. 

For the fermion, wc compare the following two paths: 

^ % V + 21^s(to^)C 

OO 

^ V + 2V^^(m2)(e + 0+1^ a2k[{ei> + i>e)m''^-^ + m{etP + i,e)m''^-^ + • • • m''^-^{ei) + V'e)]^, 

fe=2 

V' ^ il) + 2V^{m^)e 

OO 

h 

fc=2 



In order to see the structure of the supersymmetry transformation of the Af = 2 supersymmetry. 
Since we now abstain the hermiticity of the fermion, we separate the fermions into the real part and the 
imaginary part as 

e = ei+ie2, ^ = ^i+i^2, (4.36) 

where ei, £2, and ^2 arc Majorana-Weyl fermions. Now, we assume that the supersymmetry parameters 
are proportional to the unit matrix for brevity, namely 

ei,2 = jjliVxAT, ^1,2 = C^i,2'^NxN, (4.37) 

where 0^12 and c^^ ^ real Grassmann-odd c-numbers. Another assumption is that the background 
metric is fiat. 

Firstly, the transformation of the bosonic field is calculated as 



k=2 

^ 00 
k=2 

= ^(^-[^''^^-^^Je - e[r'^^r-%r^]OA., a., + • • • 

= ^(frV - eVOiAi, A^] + --- = ^(aPei + 6r'e2)[^i, A^] + ---, (4.38) 

We focus on the supersymmetry transformation of the vector field a^{x). The commutator 
represents the translation and the gauge transformation: 

[Ai,Aij\ = [idi + ai{x),idij, + aij,{x)]-\ = i{diafj,{x)) -i{di_,ai{x)) + [ai{x),a^{x)]-\ . (4.39) 

^ V ' ^ V ' 

translation gauge transformation 



55 



Therefore, the coramutator of the bosonic fields represents the right translation of the vector fields. 
However, this does not hold true of the fermionic fields. The commutator is computed as 

n 
k=2 

= -2a4(Cir^'ei + 6r^e2)V^^j + • • • • 
We have a closer look at the term il)Aj: 

^A, = ii,d, + • • • = ^{x)d, +f2x^''■■■"H^)^^, ■ • •a.a,^ [e-^^'r + ■■■, (4.40) 

where • • • denotes the omission of the non-linear terms of the fields. Therefore, each fermionic field is 
subject to the following transformation: 

[V-,, i^M^) = + • • • , [V-,, i^i]x'^''-"+'Hx) = -2a4(eir^ei + 6r^'e2)x^'"-"' {x)6''^'^' +■■■. (4.41) 

Therefore, the fermionic fields are not subject to the translation by the commutator of the supersymmetry. 

However, the argument in this section is only applicable to the function Vg{u) that can be expanded 
around u = 0, and does not hold true of the function H4.17|l . We speculate that we may find a right 
supersymmetry for such functions. 

We surmise that the fermionic fields and x*-*-* (x), which are respectively identified with the dilatino 
and the gravitino, would be massless due to the supersymmetry, when the even-rank fields a^'^i^ij...^j2,^ {x) 
are massless. It is an intriguing future problem to seek such an action more rigorously. 



5 Monte Carlo simulation of the fuzzy sphere solutions 

We have seen several generalizations of the HB matrix model to accommodate the curved-space classical 
background in Sec. 12.41 We have reviewed the author's work [46] in Sec. 13.41 about the curved-space 
background in the context of the supermatrix model. 

In this section, we address the stability of the curved-space background in the quantum sense, based 
on the author's work [60]. It has been a conundrum to discuss the stability of the fuzzy sphere in the 
quantum sense. There have been hitherto several attempts for the quantum stability. In [25], the quantum 
stability has been discussed in terms of the one-loop effective action. In [51,55,59], they performed the 
two-loop perturbative calculation of the Feynman diagram, to unravel the stability of the fuzzy sphere and 
other backgrounds. Another attempt is to investigate the quantum stability via the Gaussian expansion, 
and the first-order calculation is given in [52]. 

Here, we address the quantum stability of the fuzzy sphere via the heat bath algorithm of the Monte 
Carlo simulation, whose detail we delegate to Appendix lO Our approach is totally different from these 
foregoing works [25,51,52,55,59], in the sense that we discuss the stability nonperturbatively via the 
niunerical simulation. The Monte Carlo simulation has played a pivotal role in elucidating the nonper- 
turbative aspects of the lattice gauge theory. By the same token, the application of the Monte Carlo 
simulation to the HB matrix model has given us many interesting nonperturbativc insights. In [13], the 
bosonic IIB matrix model has been scrutinized nonperturbatively via the Monte Carlo simulation. This 
analysis [13] was followed by the supersymmetric extension [20,22] via the hybrid Monte Carlo simulation. 

Here, we scrutinize the simplest case for the matrix model with the curved-space classical solution. 
Namely, we concentrate on the three-dimensional bosonic IIB matrix model with the Chern-Simons term 
which incorporates the fuzzy sphere solution: 

S = NTr i-\ t^^^'^'']' + ^ E ) • (5.1) 

The notation is the same as that in Sec. 12.4.11 except we set p- = iV in this section. We delegate the 
miscellaneous properties of the fuzzy sphere classical solution — aL^ to Sec. 12.4.11 
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5.1 Nonperturbative studies of the fuzzy sphere initial condition 

We start with the Monte Carlo simulation for the initial condition of the irreducible representation of the 
N X N fuzzy sphere classical solution: 

4") = aL^. (5.2) 

Starting from this initial condition, we see whether or not the fuzzy sphere decays due to the quantum 
effect of the matrix model. Especially, the histogram of the eigenvalues of the Casimir operator 

Q = Al + Al + Al (5.3) 

describes the fuzzy sphere's stability visually. In the following, x denotes the eigenvalues of Q, and 
f{x) denotes the eigenvalue distribution function normalized as Jp°° dxf{x) — 1. At the initial state, the 

histogram is peaked at = ^ ^^ a^, because the Casimir is proportional to the unit matrix for the 
fuzzy sphere. If the fuzzy sphere is stable in the quantum sense, the histogram does not decay greatly 
from the peaked distribution. On the other hand, if the fuzzy sphere is unstable, the histogram undergoes 
a drastic deformation. 

5.1.1 Nonperturbative stability of the fuzzy sphere 

We first investigate the TV = 16 and a = 0.5, 1.0,2.0 case starting from the initial condition H5.2() . We 
see a clear difference between the a = 0.5 and the a = 1.0, 2.0 case. 





100 150 200 

X 



Figure 10: The instability a,t a — 0.5(left) and the stability at a = 1.0, 2.0(right) of the fuzzy sphere for 
N = 16. 

We see from Fig. ^]the instability of the fuzzy sphere at a = 0.5. The vertical lines in the histogram 
denote the radius-square of the fuzzy sphere classical solution i?^ = ^ . Now, this value is 15.9375 
for a = 0.5, and the eigenvalues of Q were originally peaked there. This fuzzy sphere decays due to the 
quantum effect and the eigenvalues are concentrated far away from the original fuzzy sphere configuration. 

However, this situation changes for a = 1.0,2.0. Fig. EH indicates that the eigenvalue distribution 
constitutes only one lump in the vicinity of the classical value of the radius-square (which is 63.75 and 
255 for a = 1.0 and a = 2.0 respectively), and that the eigenvalues are not dissipated otherwise at all. 
This implies that the eigenvalues remain constituting a thin sphere shell near its original fuzzy sphere 
configuration. In this sense, we regard the fuzzy sphere classical solution as stable. This stability is also 
reconfirmed in Section 15.3.11 by demonstrating the dynamical evolution of the single fuzzy sphere state 
from the initial condition A^j^'' — 0. 

5.1.2 First-order phase transition 

In the previous section, we have discerned that the fuzzy sphere solution. In this section, we deepen the 
above observation and demonstrate the first-order transition of the matrix model H5.1|l . 
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Firstly, we ascribe the above stability of the fuzzy sphere to the meager quantum effect of the model. 
The effective action W of this matrix model is given by 



W = -\og\^J dAe"" j . (5.4) 

The effect of the classical fuzzy sphere solution is of the order 0{a*N'^), as one can easily verify by 
plugging the fuzzy sphere solution l|2.78|l into the action (jOJ- On the other hand, the effect of the path 
integral measure is of the order 0{N'^). This comparison leads to the observation that the matrix model 
receives only a small quantum effect at 

a>C'(^). (5.5) 



To elaborate on this observation, we next calculate the following quantities, for N = 8, 16, 32 with 
the fuzzy sphere initial condition 1)5. 2|l . 

1. The action (S): This quantity gives the energy of the numerical configuration. 

2. The spacetime extent {j^TrA^): This provides us with the information of the eigenvalue distribu- 
tion, and hence the stability of the fuzzy sphere. 

3. The Yang-Mills term: {j^TrF'^^,): This quantity, as well as the spacetime extent, have been already 
investigated in the bosonic IIB matrix model without the Chern-Simons term (this amounts to the 
a = case) in [13]. For the bosonic IIB matrix model, both are shown to behave at the order 0(1). 
It is interesting to compare the behavior of this quantity with that of the bosonic IIB matrix model. 

4. We derive the following exact result from the Schwinger-Dyson equation of the matrix model: 



d 



0= / dAj^{Tr{eA^)e-% (5.6) 



dA 



Here, (a = 1, 2, • • • , N — 1) are the basis of the SU{N) Lie algebra, and the matrices are 
expanded as A^ — X^^T^ ■^lit"'- The Schwinger-Dyson equation (|5.6|l gives 



{K) = l(rr(-[A^, + 2tae^,pA^A,Ap)) = 3(1 - ^) = Ko, (5.7) 



as we prove explicitly in Appendix lC.4l While K is given analytically, the numerical computation 
of its vacuum expectation (K) is significant for the legitimacy of the algorithm. 



5. The Chern-Simons term: 



(M) = {^Tre^.pA,A,A,). (5.8) 



We have obtained the analytical result K, which amounts to 



K - ^TrF^, + 3aM. (5.9) 



While (K) itself obviously behaves as 0(1), it is interesting to discern the behavior of each of the 
separate terms M, as well as j^TrF^^. 

We postpone the plot of (K) later, and plot the other four quantities. The above power-counting obser- 
vation for the quantum effect motivates us to plot these quantities against a = ay/N, instead of a. For 
the fuzzy sphere classical solution, they respectively scale as 

l^TrF^^=0{a*N^)^0{a^), M = 0{a^ N^) ^ 0{&^Vn). (5.10) 
This leads us to plot the quantities j^(S'), j^{jjTrA^^), {jjTrF"^^) and --j={M) against a in Fig. ^ 
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Figure 11: -^{S) (upper left), ^(^TrA^^ (upper right), {j^TrF^ (lower left) and ^ 
right) against a, for = 8, 16, 32 with the fuzzy sphere initial condition (|5.2|l . 
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This plotting poses a serious consequence. This matrix model has the first-order phase transition with 
the change of the parameter a. In the above four cases, we see a discontinuity at the critical point 

(5 = a^t^-2.1. (5.11) 

We also find that this discontinuity is steeper for the larger iV. This indicates that there is a first-order 
phase transition between the two phases. 

We call the phase at a < air the "Yang-Mills phase". In this phase, the quantum effect of the 
matrix model becomes strong, and thus the fuzzy sphere classical solution is no longer stable. Rather, 
the behavior for this phase has a similarity to that of the bosonic IIB matrix model without the Chern- 
Simons term. Fig. ^2 clearly indicates that this model undergoes a smooth transition from a = 
(the bosonic IIB matrix model itself) to the small a. The vacuum expectation values of {j^TrA^^ and 
(j^TrF^^) have been analyzed in [13], where they are shown to behave at 0(1). This behavior is inherited 
in the Yang-Mills phase. We elaborate on this phase more closely in Section 15.2. II 

We call the phase at ci > air the " fuzzy sphere phase" , where the quantum effect is so small that the 
classical fuzzy sphere is stable. Now, we discern that the fuzzy sphere is stable at iV = 16, a = 1.0, 2.0, 
because the system is in the fuzzy sphere phase in this case. They are scaled as (|5.10f) . and this behavior 
is totally different from that of the Yang-Mills phase. 

This first-order phase translation is a product of the the Monte Carlo simulation, which could not be 
seen by the perturbative approaches. In this sense, we accentuate that this first-order phase transition is 
a nonperturbative phenomenon. 

We next exhibit the corresponding plot of the quantity (K) in Fig. El In the Yang-Mills phase, 
the observables {j^TrF"^^) and (Af) behave as 0(1), similarly to the bosonic IIB matrix model. The 
fluctuation of K is thus small in the Yang-Mills phase. 

On the other hand, (jjTrF^^^) and 3a{AI) both behave as 0{a^) in the fuzzy sphere phase. Nev- 
ertheless, the quantity {K) = {{j^TrF'^j^ -t- iaM)) is of the order 0(1), as is analytically verified. This 
is ascribed to the huge cancellation between these two terms in the fuzzy sphere phase. Therefore, the 
fluctuation of {K) is naturally bigger in the fuzzy sphere phase than in the bosonic Yang-Mills phase. 
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Figure 12: The plot of the quantity for the fuzzy sphere initial condition (|5.2|l . 



5.1.3 One- loop dominance in the fuzzy sphere phase 

The behavior of these four quantities in the fuzzy sphere phase gives an even more striking consequence. 
The numerical value of these four quantities coincides their one- loop results. The one- loop vacuum 
expectation values are calculated in the large- A'^ limit as 




(5.12) 



We delegate the derivation of these results to the appendices of [60]. The first term comes from the 
classical effect, while the one-loop correction comes in the second term of each quantity. The dotted 
lines in the graphs denote these one-loop result of these quantities'^ . The coincidence of the numerical 
values of these observables with the one-loop result suggests the higher-loop effects are suppressed in the 
large- iV limit and hence that the one- loop effect is dominant in the matrix model H5.1|l . 



5.1.4 Width of the eigenvalue distribution 

We next discuss the dependence of the width of the lump in the histogram of the Casimir Q on the 
parameter a and N. We define the width a as 

Jo \Jo / 

(5.13) 

Here, Ai, A2, • • • , Aat are the eigenvalues of Q. The above three expressions are all tautological. For 
the convenience of comparing the histograms for different N and a, we plot the eigenvalue distribution 
against the (x — xq) in Fig. ^1 instead of the usual argument x, where xo{N,a) are the mean of the 
observed eigenvalues for a particular N and a and are functions of A'^ and a. 

The overlapping of the different a histograms for a particular TV shows that the width of the eigenvalue 
distribution does not depend on a. We also observe that the width of the distributions for different N 
has a linear dependence on log N in Fig. 1141 The width a turns out to behave as 

- 2 log iV- 0.84. (5.14) 

This is in agreement with the one-loop calculation = 21og-/V, which is also derived in the appendices 
of [60]. 

-•^^For {jjTtF^jj), we omit the dotted line, because it is not subject to the one-loop correction. 
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Figure 13: Eigenvalue distribution of Q around its mean, for = 8 (upper left), N = IQ (upper right), 
iV = 32 (lower left) and TV = 64 (lower right), respectively. 
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Figure 14: The square of the width of the distribution a"^ against logiV. 
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5.2 Analysis of the Yang-Mills phase 



Wc have so far investigated the numerical simulation starting from the fuzzy sphere initial condition 
(|5.2I) . We have found that there is a phase transition between the Yang-Mills phase and the fuzzy sphere 
phase at air ~ 2.1 (namely, air ^ "^)- ^^^^ section, we elaborate on the behavior of the Yang-Mills 
phase, and we launch the Monte Carlo simulation from the initial condition 



= 0. 



(5.15) 



5.2.1 Eigenvalue distribution in the Yang-Mills phase 

We start with discussing the eigenvalue distribution of the Casimir Q, for N = 8, 16, 32 with a = 0.0 to 
indicate the N dependence and a = 0.0,0.2,0.4,0.6 with = 32 to indicate the a dependence. Here, 
we plot the eigenvalue density p{r) against r = (eigenvalues) in Fig. ^1 This r trivially represents the 




0.06 



0.04 



0.02 



a=0.0 
a=0.2 
a=OA 
a=0.6 



Figure 15: The plots of the eigenvalue density against the radius r, for = 8, 16,32 with a — 0.0 to 
indicate the N dependence(left), and a = 0.0, 0.2, 0.4, 0.6 and = 32 to indicate the a dependence(right). 

distance from the origin. p(r) is defined as 



(distribution of -v/ (eigenvalues) of Q) 
= ^2 ' (5-16) 

so that this may represent the eigenvalue density. p{r) is normalized so that drAirr^ p{r) — 1. 

The histogram for a — 0.0 is devoid of the eigenvalue in the vicinity of the origin, and the histogram 
soars around the radius 0.8. In this sense, there is an ultraviolet cutoff in the eigenvalue distribution. 
This phenomenon is interpreted in the following way. If a certain state is an eigenstate of the Casimir 
Q with almost zero eigenvalues, this implies that the same state is almost an eigenstate of Ai, A2 and 
^3 individually with almost zero eigenvalues. However, this violates the uncertainty principle, because 
a generic configuration of do not commute with each other. In the bosonic IIB matrix model [13], 
the vacuum expectation value (jfTrF^^) is of the order 0(1). Therefore, the scale of the uncertainty 
[A^, Ajy] is also of the order 0(1). The above histogram is consistent with this consideration. 

The histogram for a = 0.0,0.2,0.4,0.6 with iV = 32 demonstrates that the eigenvalue distribution 
undergoes a gradual change as we vary a within the Yang-Mills phase. 

5.2.2 Phase transition and the hysteresis 

We next investigate the phase transition of the simulation initiated from H5.15|l . We find that the critical 
point differs in the initial condition (|5.15|) from that for the fuzzy sphere initial condition (|5.2|l . When 
we initiate the simulation from = 0, the critical point is found at 

4;;) ~ 0.66. (5.17) 

We plot (j^TrAl) for N = 8, 16, 24 with the two initial condition and in Fig. UHl The 

(u) 

important difference in the Yang-Mills phase is that the critical point acr is independent of N. Since the 
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Figure 16: The plot of (j^TrAl) against a, for iV = 8, 16, 24. 



critical point acr is larger than that for the fuzzy sphere initial condition air ^ for the sufficiently 
large N, we call the critical point for the Yang-Mills (fuzzy sphere) phase 'the upper (lower) critical 
point' respectively. The above hysteresis cycle is realized for the sufficiently large N. For the N — 8 
case, the lower critical point is approximately ac) ^ 7i ^ 0.742 • • •, which is larger than air^ . In this 
case, the hysteresis disappears, as is discerned from Fig. ^] On the other hand, this is not the case with 
N = 16, 24 in which the hysteresis cycle is realized. 



5.3 Metastability of the multi-sphere solutions 

In this section, we investigate the multi- fuzzy-sphere solution of the matrix model H5.1|l . It accommodates 
the following multi-fuzzy-sphere solution 

A^^ax diag(L|:^),4"^),...,4"''^). (5.18) 

Here, is the rta-dimensional irreducible representation of the SU{2) Lie algebra. The whole size of 

the matrices is = X]a=i "-a- For this multi- fuzzy-sphere solution, the eigenvalue distribution of Q 
is peaked at 

rl = ^{nl-l), (a = l,2,...,fc). (5.19) 
The classical value of the action is easily calculated as 

a=l 

The value of the action H5.20|l for the multi-fuzzy-sphere condition is bigger than that for the single fuzzy 

sphere state (namely, S — — " ^ )■ Therefore, the multi- fuzzy-sphere state can be realized as a 

metastable state. In the following, we report the results of the simulations to unravel this metastability. 
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5.3.1 Evolution of the multi-fuzzy-sphere solution 



Firstly, we report the simulation with the initial condition (|5.15l) , and observe how the fuzzy sphere state 
evolves. Here, we focus on the large-a case in the fuzzy sphere phase = 16 and a = 2.0. We plot the 
history of the vacuum expectation value of the action (S) , and the eigenvalues of Q against the sweeping 
time in Fig. 1171 This indicates the process in which the multi-fuzzy-sphere state is generated dynamically. 
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Figure 17: The history of the vacuum expectation value of the action (5) (left), and the eigenvalues of 
Q (right) against the sweeping time, for = 16, a = 2.0. 

Firstly, we start to see a tiny plateau in the graph of (S) corresponding to the multi-fuzzy-sphere state 





(10) 



When we further undergo the sweeps, this multi-fuzzy-sphere state falls into 

the lower-energy state = a ^ ^-^^^ j . In this way, the state falls off to the lower-energy state 

\ J 



A,, = a 



^(4.3,24) 



\ 412,13,14,15) J- 

Finally the system arrives at the one-fuzzy-sphere irreducible representation ~ aL^}^'' . Namely, 
the one-fuzzy-sphere state is dynamically generated in the course of the thermalization. When we further 
undergo the sweeps, the system is stuck in this one-fuzzy-sphere representation, and never falls off to the 
other state any more. This observation reinforces the stability of the one-fuzzy-sphere state in the fuzzy 
sphere phase. 

5.3.2 Metastability of the fuzzy sphere solution 

In this section, we obtain an insight into the metastability of the fuzzy sphere state, by focusing on the 
dependence on fc, a and N. We adopt the initial condition of the n — ni = n2 = ■ ■ ■ Uk = ^ case in 
(|5.18|) for brevity. Namely, we initiate the simulation from the multi-fuzzy-sphere state 

4"' -"^"^^Ifcxfc. (5.21) 

We have observed in the previous section that the multi-fuzzy-sphere state falls into the lower-energy 
state in the course of the thermalization of the system. Firstly, we observe how the multi-fuzzy-sphere 
classical solution (|5.2HI decays, for = 16, a = 10.0 and A: = 2,4, 8 in Fig. El This indicates that the 
multi-fuzzy-sphere state (|5.21|l decays faster for larger k (namely, for the smaller representation of the 
size n = This is a natural result, and we elaborate on this dependence on k later. 

Nextly, we show the first-order phase transition and the one-loop dominance of the k = 2 reducible 
representation. The observables ■^{jj^TrA'j^) , (jjTrF^^,) and -^(Af) are calculated in the appendices 
of [60] as 

(S) 1 , 1 „ ,9, 1 , 1 „ , <5^ 1 , , 1 , 

¥—24^ + ^' N^N^^''l^ = W^-a^' <]v^^(^-)^> - 2P' 7N^''^--6k^ + a- (^-22) 
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Figure 18: The eigenvalue distribution of the Casimir Q, as the multi-fuzzy-sphere solution ()5.21|l starts 
to decay, for iV = 16, a = 10.0 and A: = 2 (left), fc = 4 (middle). A: = 8 (right). 

for the one-loop effect around the multi fuzzy sphere (|5.18|) with n = ni = n2 = ■ ■ ■ = Uk = 
launch the simulation from the initial condition H5.21|l for fc = 2. We plot these four observables in Fig. 
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Figure 19: ^(S) (upper left), ^(^Tr^) (upper right), (iTrF^ (lower left) and ^(M) (lower 
right) against a, for N = 32, 64 for the multi fuzzy sphere state (|5.21|) . 



We have seen in Fig. 1181 that the multi fuzzy sphere decays in the course of the thermalization. Here 
we undergo a sufficiently short sweep in the fuzzy sphere phase, such that the multi fuzzy sphere H5.21|l 
may not decay. In the fc = 2 multi fuzzy sphere, we find two consequenses similar to the fc = 1 irreducible 
representation. Firstly, we note that there is a first-order phase transition. This time, the critical point 
lies at 



_ 2.8. 



(5.23) 



Secondly, the fc = 2 multi fuzzy sphere has the one-loop dominance in the fuzzy sphere phase in the 
large- iV limit. 
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We can understand the critical point (|5.28(l from the one-loop dominance. In the fuzzy sphere phase, 
the spacetime extent (j^TrAj^) behaves at the one-loop level as (j^TrAj^) = ^ — -i^ in the fuzzy 
sphere phase. This is positive only when a > \/2k. In this sense, we obtain the lower bound on the 
critical point from the one-loop dominance as ficr'" ^ 0{\/k). The critical point (|5.11() for fc = 1 and 
(|5.23|l for k = 2 are consistent with this observation. 

We next discuss the sweeping time r that takes for the initial condition H5.21|l to decay more qual- 
itatively. To this end, we plot log(Tfc'^) against log a for = 8, 16 in Fig. 1201 This indicates that the 
sweeping time for the multi-fuzzy-sphere state to decay depends on a and k as 

r-a^fc-^ (5.24) 

Fig. 1201 indicates that this power law is independent of the size of the matrix N. 
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Figure 20: The plot of log(Tfc^) against log a for = 8, 16. 

This is in contrast to the intuitive observation that the decay probability P obeys P ^ e^^ — 
exp(— TV^a^fc^^); namely the sweeping time for the decay is p ~ exp (iV"'a'*fc^^) . It is interesting to 
ruminate on the reasoning for the deviation from this intuitive observation. 

5.4 Miscellaneous future directions 

We conclude this section by listing the future approaches to this work. There are a lot of interesting 
future works in this direction. Firstly, this toy model is expected to serve for the dynamical generation of 
the gauge group. As we have reviewed in Section f2.4.1l the expansion around the reducible representation 
(|5.18|l for m ~ ■ ■ ■ rik = ^ gives the noncommutative Yang-Mills theory with the U (fc) gauge group. Iso 
and Kawai [14] suggested that the dynamical generation of the gauge group is ascribed to the cluster 
distribution of the eigenvalues. By examining the eigenvalue distribution of this model, we may gain 
insight into the dynamical generation of the gauge group. 

Secondly, it is also exciting to extend our analysis to the supersymmetric case. For the simplest 
case, we start with the four-dimensional supersymmetric model with the three-dimensional Chern-Simons 
term^^. It is intriguing to investigate how the supersymmetry affects the results we have obtained for the 
bosonic case. 

Thirdly, the extension to the higher-dimensional extension is also an intriguing issue. We can extend 
our analysis to the higher-dimensional fuzzy sphere, as we have reviewed in Section 12.4.21 There has 
been hitherto no works which discuss the quantum stability of the higher-dimensional fuzzy sphere S'^'', 
since it is much more mathematically involved than the simplest S*^ fuzzy sphere. On the other hand, it 

'^^Austing and Wheator [57] corroborated that the three-dimensional action 12.761 has a divergent path-integral. There- 
fore, the numerical simulation of the three-dimensional model is not possible. 
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is easy to extend the algorithm for the matrix model with the three-dimensional Chern-Simons term to 
the higher dimension. We expect that the Monte Carlo simulation may open the door to the study of 
the higher-dimensional fuzzy sphere. 

The fourth future direction is to extend our analysis to the general homogeneous space, which is 
proposed by Kitazawa in [44]. By investigating several homogeneous spaces (such as CP-^ = SU (3)/C/(2)), 
we may be able to understand which homogeneous space is the most favored in the quantum sense. We 
expect that this direction may give us a clue to the compactification of the spacetimc. 

In this way, the Monte Carlo simulation we have described here gives us a lot of exciting prospects. 

6 Conclusion 

In this thesis, we have reviewed the author's works [26, 38, 46, 60] about the relation between the grav- 
itational interaction and the matrix model. Now, it is believed that the 'Theory of Everything', which 
unifies all the interactions in the nature, is realized by the siiperstring theory. The large- TV reduced model 
is now regarded as the promising candidate for the nonperturbative formulation of the superstring theory. 
If we build the nonperturbative formulation of the superstring theory, this will be the right framework 
to unify all interactions. Then, we will be able to solve all the riddles of the Standard Model, such 
as the dimensionality of our spacetime, the generation of the quark, the gauge group and the eighteen 
parameters of the Standard Model, from the parameterless theory. 

In order to arrive at this ultimate goal, it is an important step to see the correspondence between the 
matrix model and the gravitational interaction. In the future, it is indispensable to elucidate the more 
manifest correspondence with the gravitational interaction. For example, it is interesting to pursue the 
matrix model that is equipped with the local Lorcntz symmetry and reduces to the supergravity in the 
low-energy limit, by inheriting the idea we described in Section 4. If we find such a matrix model, it will 
be certainly a better extension of the IIB matrix model, in the sense that it has clearer relation to the 
gravitational interaction. In addition, it is also interesting how the matrix model realizes the curved-space 
background. The curved spacetime is also an essential feature of the general relativity. So far, it is known 
that the IIB matrix model, per se, realizes the curved-space background through the condensation of the 
graviton. In addition, several alterations of the IIB matrix model have been advocated, so that they may 
be equipped with the curved-space manifold. On the other hand, the alterations of the IIB matrix model 
is merely limited to the simple manifolds; such as the 5^*^ fuzzy sphere, fuzzy torus or CP^ manifold. It 
is exciting to pursue how the matrix model realizes the general curved spacetime more manifestly. 

We finally mention the relation between the IIB matrix model and its miscellaneous extensions. In 
the quantum field theory, some different models which share the same symmetry are equivalent in the 
continuum limit, known as universality. We expect that the similar mechanism may hold true of the 
matrix models and hence that various matrix models may have the same large N limit. Thus, there is a 
possibility that these new extensions are equivalent to the IIB matrix model. We can expect that both 
the IIB matrix model and the new extensions arc equally authentic constructive definition of superstring 
theory in the large N limit. It is an intriguing future work to elaborate on this conjecture of the 
universality. 

Mankind has yet to grasp what is the true 'Theory of Everything'. We must exert more effort to 
arrive at the answer to this ultimate and most difficult question of the elementary particle physics. 
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A Notation 

A.l Definition and properties of the gamma matrices and the fermions 

Basically, we follow the same notation as for [27]. Throughout this thesis, the metric of the Minkowskian 
spacetime is taken so that the time component should be minus and the space components should be 
plus. Namely, the metric is defined as 

r?'''' = diag(-l,+l,---,+l). (A.l) 

The gamma matrices are defined so that it may comply with the following Clifford algebra: 

{r^r''} = 27?''''l32x32. (A.2) 

Here, the gamma matrices F'* are the 32 x 32 real matrices and the indices fjb,v, - ■■ run over 0, 1, • • • , 9. 
The explicit components are given by 

r» = 1...,. « (-...) = ( ^i^^ ) , = « .3 = ( ^; ; ) . 

Here, are the nine-dimensional (and thus 16 x 16) gamma matrices, and obey the Clifford algebra for 
the nine- dimensional Euclidean spacetime {7^, 7'} = 26P'^liex 16, where p,q, - ■■ run over 1, 2, • • • , 9. 
The transpose of the gamma matrices are given by 

^(F°) = -r°, ^{TP)=+rP (p = l,2,---,9).. (A.3) 

The chirality matrices r^'^(= F") are defined by 

r« = rW...F9 = ii6xi6®ai = f ^ ° I16XI6 y ^A.4) 

V Il6xl6 J 

For the indices A, S, • • • = 0, 1, • • • , 10 in the eleven-dimensional Minkowski spacetime, the gamma ma- 
trices comply with the Clifford algebra {r^,r-^} = 277^^132x32- 

Next, we define the charge conjugation matrix C. C should satisfy the conditions 

CF'* = -^(r^)C, C+^C = 0. (A.5) 

In our case, C is identical to F°. 

The anti-symmetrized gamma matrices are abbreviated as 



k\ 



^ y sgn(CT)F^<'(i) r^-(2) • • • F^-C") . (A.6) 

Here, iS„ is the set of the permutation of the integers 1, 2, • • • , fc. 
A. 1.1 Chirality of the Weyl fermion 

We summarize the notation for the chirality of the fermion. The left-hand (right-hand) fermion is defined 

as 

1 + F* 1 - F" 

V'L = V', V'K = V'- (A.7) 
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The left-handed (right-handed) fermions satisfy the following identities for the ten-dimensional gamma 
matrices: 

XLT^'-'^^'eL = XRr^'-''''eR = 0, xlT^'-^^''+' eR = xrT'''-''''^'^l = 0. (A.8) 
Here, we prove only xl^^^ ^^''^l — 0, so that the rest can be verified likewise. 

^ ^;^i±l!l_l!rOr'^i-M-, = o. (A.9) 

A. 1.2 Epsilon tensors and the duality of the matrices 

Here, we define the anti-symmetric epsilon tensors as 

eoi-fc = 1, (so that e^^- *^ = -1). (A.IO) 

This definition leads the duality relations of the gamma matrices 

fc(fc+i) 

i - (10 -fc)! ' ^A.+i-Ar (A.iij 

A. 1.3 Multiplication law of the gamma matrices 

The anti-symmetrized gamma matrices obey the following formula of the product. 

+ (-l)('"-i)+("-2)„C2„C22!r7l^^"''r;-4^'''r-43-^™l^^-^"l (A.12) 



A. 1.4 Flipping property of the Majorana fermions 

Here, we introduce a useful formula for the flipping of the Majorana fermion. The gamma matrices has 
the following property: 

pO(TpA)pO ^ pA^ (^^^3) 

This leads to the following relation: 

pOj-TpAi-.-Afc-jpo ^ (-i)''^i(r°(^r"*'')r°) • • • (r°(^r^^)r°) = („i)fe-ir^'='"'*i 

= (-l)'=-i(-l)^^r'^i-^^ = (-1)^ — ^ — '-Y^i-Au^ (A.14) 
Then, the Majorana fermions are flipped as 

^pAi...A.g ^ ^(xr'4i-^'=e) = -^e(^r'4i-^'=)(i^r0)x = ~^e(r")2(^r^i-'4^)(r0)x 

= (-l)^eT-4--^'=x. (A.15) 
Namely, the sign is minus for k = 1,2,5, and plus for k = 0, 3, 4. 



A.1.5 Proof of the Fierz identity (|2740|l and (1278 

We give a proof of the formula of the Fierz identity H2.40|l and (|2.87|) . which contributes to the proof of 
the supersymmetry algebra of the matrix model. 

We start with the ten-dimensional formula (|2.40() . The fermions ei and £2 are Majorana- Weyl fermions, 
and effectively act on the 16 x 16 space of the gamma matrices. The left hand side is written as 

e-ir,Vr^''e2 = (ei)"(r.)„%(r^'')/e25 = -(r^'')/(eie2)5"(r.)„''#. (A.16) 



69 



Here, the minus sign emerges because we have flipped the order of the Grassmann odd fermions. The 
matrix (€162) 5" can be decomposed as 

Since the fermions ei and 62 have the same chirality, the terms of the even rank vanish due to the relation 
(EH. We ignore the rank-3 term because we are interested in the difference eiTi.ipT'^'^ e2 — ^2^vi^^^'^ ^i, 
and the rank-3 term does not contribute due to the relation ljA.15|) . This leads us to focus only on the 
rank-1,5 terms. The relevant gamma matrices are computed as 

- 7(rpr^ - 77%) - 977^p = 7r,r^ - 1677%, (A.is) 
V p-^^...p^Ti, = r'* (^Tiypi-.-ps + 577y[pj^rp2...p5]) 

= ^{^"n^ ^ypi---p5 ~^ 'n lyT'^pi-.-ps — Sjy^fp^r^ P2---P5]) + ^(-'^''pi-'-PS ~ '?^[P2-'^P1P3P4P5]) 

= -rp,...p/ + 5r7'^[p^rp,...p,] = rp,...p,r^ (A.19) 

When we substitute ljA.17p . ljA.18|l and (|A.19|I into (|A.16p . the computation of the Fierz transformation 
goes as follows: 

= --L(e-irpe2)(r^''r'T,)V - ^^^(eirp,..p,e2)(r'^"r''--'''^r,)^ + (rank 3 term) 

= e^T^^e2^ - -^eir''e2rpP'V - — ^eir''--''=e2rp,...p,r'V + (rank 3 term). (A.20) 
16 16 X 5! 

We next go to the proof of H2.87|l for the three-dimensional spacetime. This time, the fermions are 
Majorana ones. The proof goes in the same way as for H2.40|l . Here, we note that the 2x2 matrices M are 
decomposed by the Paulian matrices as M = i(TrM)l + ^{TrMa^)a^. Then, the Fierz decomposition 
is verified as follows: 

eicr piper 62 = {ei)°'{(Tp)a'^ipp{ap„)^^{e2)s = -{cr,ii^)',^i^i<^2)s''io-f^)a'^i>f3 ^ -'^{lia^t2){(Jp,^a^ap)il}. 

(A.21) 

Here, we drop the contribution of the rank-0 term (eile2), because we are interested in the difference 
^lO' piper pi,e2 — 620' pipa pi,ei and the rank-0 term cancels. Likewise, the latter formula is rewritten as 

licrp'ipape2 = (ei)"(a'p)a^V'/3(o-p)7''(e2)A- = -(o-p)7''(eie2)5"(a'p)a^V'/3 ~ -^(eia-xe2)(a'^crxCrp)V'- 

(A.22) 

The relevant Paulian matrices are calculated as 

erpv^r^Op — i^ipypCr p(^5y^p ^- it-^prjO'q) — ^e-^i^pOp {^ur^^px ^i^X^P^)(^P'? ^ ^^prjfi^rp^ 

^^X^P^P ^^x^P^P ^^X ^ ^^X ^pri^pri ^^i^X' (A. 23) 

=3 

er pCrx^ p- ^pi^xp ^ ^^xpv^v) ^ ^x ^ ^^xpvi'^^pvp^p ^ ^pv^ ^x i^w^xp ^vx^vp^^ p 

= - 3 + 1) = -ax- (A.24) 

Plugging HA.23|) and HA.24|) into (jA.2ip and (|A.22p . we complete the proof of the Fierz identity (|2.87|l . 
A.1.6 Proof of the Fierz identity ((TTn\i 

We next verify the Fierz identity H2.10|l . This formula plays an important role in the supersymmetry 
of the super Yang-Mills theory or the Green-Schwarz action of the superstring. This holds only for the 
following cases: 
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1. d = 3, and ^1,2,3 are Majorana. 

2. d — 4, and "01,2,3 are Majorana or Weyl. 

3. d — 6, and V'1,2,3 are Weyl. 

A. d— 10, and 4'i,2,?, are Majorana- Weyl. 
To substantiate this relation, we carry out the following decomposition. 

(er^Vi)(V52r^V'3) + (er''V2)(^3r^0i) + (ePV3)(V5ir^V'2) 
+ (e)„(rOr'^)„5(V'3)5(V'i)/3(r°r^)0^(V2K. (A.25) 

Namely, our job reduces to verifying the vanishing of the following quantity: 

(rOr'^)„^(rOr^)^,- + (ror^)„^(r°r^)50 + (rOr'')c.5(r°r^)/3^. (A.26) 

We multiply (xi)7 and (x2)(5 with (|A.26|) to obtain^^ 

{T>^UCxiUT^.X2) + (r''xi)arx2r''r^)/3 - {T^X2).CxiT^T^)f,. (A.27) 

We verify the vanishing of (|A.27|I by decomposing them by the gamma matrices. For example, the 16 x 16 
matrices for the ten-dimensional case are decomposed in the same way as (IA.17|) . This leads us to verify 
the vanishing for the components of each rank. 

For (r^)^^ of the first term, it is trivial that only the rank-1 components contribute. For the sec- 
ond and third terms which entail the fermions Xi,2, we can limit the argument at most to the rank-5 
components. We note that the rank-0,3,4 components do not contribute because of the anti-symmetry 
of TqTp^...p^ (for k = 0,3,4). The even-rank components do not contribute for the Weyl fermion either, 
because of the Weyl projection. We list up which rank we should consider in the following: 

1. For d ~ i Majorana case, only the rank-1 matrices may contribute. 

2. For d = 4 Majorana case, only the rank-1, 2 components may survive. For the Weyl case, only the 
rank-1 may survive. 

3. For the d = 6 Weyl case, we consider the rank-1 components. 

4. For the d ^ 10 Majorana- Weyl case, we consider the rank-1, 5 components. 
For the rank-1 case, we multiply (Tp)i3a, and we find that this quantity vanishes: 

rr(r^rp)f xiror,,x2) - rx2ror^rpr^xi) + rxiror,,rpr^x2) = i-v + 2{d - 2))f xirorpX2) = o. 

V is the size of the gamma matrices, and d is the spacetimc dimensionality. Here, we utilize the following 
formula for the general rank-fc gamma matrices Tp-^...p^,: 

rprpj...p^r^ ~ (r^pi-.-pfc + fc'7^t[pirp2...p^])r'' — ^{^i^i^ pi---Pk ^''i[pj' tiP2---Pk\) ^''Ip.ipS P2---Pk\p. 

= (-l)'=(d-2fc)rp,...p,. (A.28) 

The rank-1 terms cancel for d = 3, 4, 6, 10, where V — 2{d — 2) = 2, 4, 8, 16 respectively. 

Next, we go on to the generic rank-fc (fc = 2, 3, • • •) case. We multiply {Tp^...p^)f:ja to obtain 

- rx2ror^rp,...p,r''xi) + rxirorprp,...p,r^x2) = 2f xiror^rp,...p,r''x2). (A.29) 

Due to the formula of the gamma matrices l)A.28|l . these contributions actually vanish for k — 2,5 for 
d = 4, 10 respectively. 

fc(fc+i) _ 

'The flipping property Xl^pi---PkX2 = (—1) ^ X2^pi--pkXi holds true only for the Majorana fermion, because 

k{k + l) 

X = ^xr" holds only for the Majorana fermions. However, the property ^Xi^^'^pi---PkX2 = (— 1) ^ '^X2^'^^Pi --PkXl 
holds true without assuming the Majorana fermion. 
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A. 2 Supermatrices 

This section is devoted to introducing the definitions of the notion of supermatrices. In treating super- 
matrices, there are many points we should be meticulous about, because what holds true of ordinary 
matrices is not applicable to the supermatrices. 

A. 2.1 Transpose 

We first introduce a notion of the transpose, emphasizing on the difference from the ordinary matrices. In 
considering such objects, it is extremely important to settle the starting point, because the other notions 
are defined so that they are consistent with this starting point. 

Transpose of Vector 

The guiding principle in considering the transpose of the supermatrices is the transpose of the vector. 

• The guiding principle is that the transpose of a vector is defined as 

T 

= {xi,---,Xn). (A.30) 

We denote {x^} as the components of v, and these components mean both bosons and fermions. 

• We define the vector as w = ^ ^ ^ where 77 and b are fermionic and bosonic real fields respectively. 

Transpose of Supermatrices 

The transpose of supermatrices must be defined so that the definition is consistent with the transpose 
of a vector. Therefore the transpose of a supermatrix must satisfy 

'^{Mv) = '^v^M, (A.31) 

where M is a supermatrix and w is a vector. Following this rule, the transpose of a supermatrix is defined 
as follows: 

• a and d are bosonic (i.e. Grassmann even) m x m and n x n matrices, respectively. 

• (3 and 7 are m x n and n x m fermionic (i.e. Grassmann odd) matrices, respectively. 

• r] (&) denote the upper m (lower n) bosonic (fermionic) components of the supervector respectively. 
(Proof) This can be verified using the very definition of the transpose. 

Then the transpose of this vector is by definition 

'^{Mv) = fri'^a + ^b^P, + '^b^d) = '^v^M. {AM) 

The point is that the sign of 777 has changed because these are Grassmann odd. Noting this fact, we can 
read off the result (|02|l . (Q.E.D.) 

We have one caution about the transpose of the supermatrix. The transpose of the transpose does not 
give an original matrix. This 'anomalous' property can be immediately read off from the definition of the 
supermatrix ljA.32|l : 
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Transpose of Transverse Vector 

We have seen an important fact that the transpose of the transpose of a supermatrix does not give 
the original supermatrix. In fact, the same holds true of the transpose of the transpose of a vector. 
Conclusion coming first, the definition is 



h 

We confirm that this is actually a well-defined settlement. This notion is defined so that 

^{yM) = ^M^y, (A.37) 

where y is a. transverse vector y — (r/, b) and M is a supermatrix Af = ^ ^ ^ ) ' compute both the 
L.H.S and the R.H.S and verify that they actually match if we follow the above definition. 



We have used the fact that in the transpose of fF, we must multiply —1 because a fermion jumps 
over another fermion. 



R.H.S. : ' AP y = 



-^7 \ f \ _f ^^a^rj-^^^b 

'^d ) [ ^b ) - [ -^/?^77 + ^7^6 



Thus we have verified that the above definition of the transpose is consistent with the condition "^{yM) = 
M^y. As we have mentioned before, the transpose of the transpose of a vector does not give the original 
vector, because 



A. 2. 2 Hermitian Conjugate 



We introduce a notion of hermitian conjugate of the supermatrix. This notion is much simpler than the 
transpose, and we do not see anomalous properties as emerged in the transpose. The starting point of 
this notion is 

• For a fcrmionic number a, P, the complex conjugate is {ctf3y — (/3)^(a)^. 

• For a vector u = ^ ^ ^ , the complex conjugate is ^ ^ ^ — (rj^b^). 
Under this definition, the hermitian conjugate of a supermatrix is defined as 

F»„=(; •2), M.= (;; 2;). (A.39) 

(Proof) The guiding principle to determine the hermitian conjugate of a supermatrix is the condition 

(MvY = (u^)(Af^). (A.40) 

For Af = ^ ^ and w = ^ ^ ^ , {MvY is computed to be, utilizing the definition for the vector, 

(A//«)t = y = ((ary)t + m\{^rj)^ + {db)^) - {rj^ + b^ fil rj^^ + bU^ ) . (A.41) 

The explicit form of the hermitian conjugate of a supermatrix can be read off from HA.41|I . and this completes 
the proof. (Q.E.D.) 

The definition of the hermitian conjugate of a transverse vector is now straightforward. This is given 

by 

y^ = iv,by=( )• (A.42) 
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It is straightforward to verify that this definition is consistent with the guiding principle 

(yM)^ = M'fy^ {A A3) 



because {l.h.s.) = {r.h.s.) 



A. 2. 3 Complex Conjugate 

We define a notion of complex conjugate for a supermatrix. The guiding principle to define the complex 
conjugate is to require the matrices and the vectors to satisfy the condition 

{Mvy=M*v*. (A.44) 

In order to satisfy this condition, we define the complex conjugate of the vectors and the matrices as 
follows^®. 

V* "if f i;)t, M* (^M)t. (A.45) 
It is straightforward to verify that this definition is consistent with the guiding principle ljA.44|) . 

{MvY = {(^v)(^M))^ = {^M)^{^v)^ = M*v*. (A.46) 

Combining the results obtained in the previous section, the explicit form of the complex conjugate of the 
vectors and the matrices are 

M- = ( ; 2 )■ - ( :;.*)."-( '0' ^ ( ? ) ■ = <"■"• = <-•''■'''■ 

This is clearly consistent with the guiding principle {Mv)* = M*v* because (l.h.s.) — {r.h.s.) — 
a*ri* +I3*b* 
-J*ri* +d*b* 

We have following properties which relates the transpose, hermitian conjugate and the complex con- 
jugate. 

(Prop) (1)^M = (M*)t, (2)Aft = ^^(M*). 

(Proof) These properties can be verified by noting that the hermitian conjugate of the hermitian conjugate 
gives back the original quantity, which can be readily verified by definition. 

1. (Af*)t ^ ((^M )t)t = ^Af . 

2. -^(Af*) = {{M*)*)^ = AI^ . In the last equality, we have utilized the fact that, for a supermatrix, 
(Af *)* = Af , which can be readily verified from the explicit form of the complex conjugate HA.47|I . 

This completes the proof of the above properties. (Q.E.D.) 

Now we are ready to answer the question : what do we mean by 'a supermatrix is real ?'. In consider- 
ing physics, we must take into account the reality condition. We utilize supermatrices in the context of 
expressing the action of superstring theory, the action must be real, and we are required to solidify the 
definition of the reality of a supermatrix. 

def 

(Def) A supermatrix A/ is real <^ M is a mapping from a real vector to a real vector. 

This statement is equivalent to, for the above definition of complex conjugate, 

M* = M. (A.48) 



i^Be careful about the faet that ^(Mt) is different from (^M)t = M*. They are computed to be ^(Mt) = 
(7* d* )' ^^''^^ ~ ( ''^ ) ^"^^ ^(y^) ~ ^'^'^ these are not the complex conjugate of the vectors or 



the matrices. 
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This can be verified by noting the starting guiding principle that M is designed to satisfy [Mv)* ~ M*v* . 
If M* = M is satisfied, (Mv)* is a real vector if v is real, because 

(Mv)* = M*v* " '^=^^^ M*v ^^ = '^' Mv. (A.49) 

The relationship ljA.48p tells us the conditions for the components of M to satisfy. Noting the explicit 
form of complex conjugate (|A.47p . we can derive a* — a, d* = d , (3* — (3 and 7* = —7, id est, 

• a, P,d should be real. 

• 7 should be pure imaginary. 

B Calculation of the Seeley-de-Witt coefficients 

This appendix is devoted to introducing the basic technique of computing the Seeley-de-Witt coefficients 
in the heat kernel expansion. There are a number of ways to compute the coefficients, and here we focus 
on the calculation based on the Campbell-Baker-Hausdorff formula. We consider the general elliptic 
differential operator 

In this section, we consider the d-dimensional spacetime, and the indices ••• run over ••• = 
1, 2, • • • , d. These indices are allocated for the curved spacetime. 

We consider the trace of the large N matrices in terms of the heat kernel. The trace of the operators 
are expressed using the complete system as 

Tre-^^' = J d'^xixle-^^'lx), (B.2) 

where the bracket |a;) and {x\ satisfies J^x — 1- However, it is difficult to consider the trace of a 

general operator, and we regard the operator as the sum of the first term of (jB.ip and the perturbation 
around it. This is a famous procedure, and the perturbation is expressed in terms of the Seeley-de-Witt 
coefficients. 

It is well known that the Green function is computed to be 

(x| exp [rg (,) — — j |,) = ^ exp j . (B.3) 



Therefore, its trace is easily derived as 



T.(exp(rg'^(y)A_l 



{2t:t) 

On the other hand, the heat kernel expansion of the general elliptic operator IjB.ip is expanded as 



Tr(e-"^') = / d'^x{x\e-^^' \x) = / d'^x^^iao + rai + r^aa + •••)• (B.5) 
J J (27rT) 2 

These coefficients ao,ai,a2, • • • are called "the Seeley-de-Witt coefficients". (|B.4|) trivially gives oq = 1. 
However, the other coefficients give nontrivial results. In the following, we derive the first order ai for 
the general elliptic operator (|B.1|) . 

To this end, we divide —tD^ as —tD^ = X + Y, where 

Y = r ({g^^{^x)~g^^{y))^^-^+A\x)^^+B{x)\ . (B.7) 
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From now on, wc compute the exponential exp(X + Y) via the Campbell-Baker-Hausdorff formula, while 
the computation is rather involved: 

e^e^ =e^p(^A + B+ ^[A, B] + ^{[A, [A, B]] + [B, [B, A]]) + --^. (B.8) 
Since wc know that {x\e^\y) = — ^^^t- exp (— ^^(a; — 2/)*(a; — y)-^gy (y)), the quantity in question is 

(27rT) ^ 

computed as 

^x+Y^-x ^ ^^p(y+^[X,Y] + ^{[X + Y,[X + Y,-X]] + [-X,[-X,X + Y]]) + ..^ 
= exp (^Y + Y] + ^{2[X, [X, Y]] - [Y, [Y, X]]) + ■ 
= 1 + y + + [X,Y]] + ^[Y, [X,Y]] + ■■■ 

+ + [x,Y]] + [x,y]] + ■■■)' + 



= i + Y + ^[x,y] + ^[x,[x,y]] + Iy^ + ^ [X, F]2 + ir [X, F] + i [X, y]y + . 



(B.9) 



Before we enter the computation of the quantity (a;|e^+^|t/), we summarize the formula of the differen- 
tiation of : 

flp^ 1 

= -^^{x-yy9,Ay)e'', (B.io) 



e 



(B.ll) 



_^j.i2dxis ^ - y)\gi,i^{y)gi^i{y) + gi^i^{y)gi^i{y) + Qi^i^^^^ 



d^e^ ( 1 

-■^{x- yY'i^ - yyH9iit2{y)9i3hiy)9iii2iy) + 9t2iAy)9iih{y)9ui2{y) ^ 9tiiAy)9t2h{y)9iAi2{y) 
+9iiu {y)9i2h {y)9i3i2 (y) + 9i2i4 {y)9tih {y)9i3i2 (y) + .9*3*4 iy)9nh {y)9i2i2 (y)) 

+ Y^(a; - vYHx - yt{x - yY^ix - yy'9niAy)9i2i2{y)9i3h{y)9uh{y)^ e^- (B.13) 

Computation of Fe^ 

We start with the computation of the easiest case: 

Ye- = r(igHx)-9Hy))^^+A^ix)^+Bix))e- 

= T {{g'\x) - 9'^{y)){-\9iM + ^(x - yt{x - ytguAy)9oiM) 
+ B{x)-'^A\x-yy9i^{y)^ e^. 

(B.14) 

Therefore, the trace is obtained by 

Tr{Ye-) = [ d'^x{x\Ye-\x) = [ d'^x^^!^B{x). (B.15) 
J J {2'kt) 2 
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Computation of F]e^ 

We next go on to a bit more complicated case, and we compute the operator [X, Y] itself: 



dx^^ dx'^^dxi^dx^^ dx^^dx^'^ dx^^dx^^ 



Therefore, the trace is computed to be, with the help of the formulae ljB.13|) . 
rr(i[X,r]e^)= j d'x{x\\[X,Y]e''\x) 



(B.17) 



Computation of [X,Y]]e^ 

We compute the operator [X, [X, Y]] as 



da;*-'! dx^^ dx^^ dxi^ dx^^ 

A ii^2i \ k,k2( w d^g^^J'^jx) d^ 
+ 9 [y)9 ^y'^dx^^dx'-^dx''^' dx'^-'dxi^dxJ^ 

iii2( kik2( d'^gi^^'^{x) d^ d'^A^x) d^ 

+9 ' ^{y)9 ' ' <• dx^^ dx^^dx^^ dx^^ ' dxi^ dxi^ + ^ ' ^^^'^ ' ^ ^'^'^ dx^^dx^^' dx^^-dx'^^dxi 

+4g''*' (y)^''' ^y^^dx^^dx^^dx^^^ dx^^dxi ^^^^dx^^dx^^dx^^dx^^^ d^ 



dx^^dx^^ dx'^^dx'^^ dx'-^ dx^^ dx'^^ dx'^^ 

Therefore, the trace is computed as 



Tr(i[X, [X,r]]e^)- J d'^xixl^iX, [X,Y]]e''\x) 



(27rr)5 I V6 dx^^dx^^ ' 3^ dx^dx^ ^ 

r2 r ^„»l»2/^^„JU2/^^„/clfe2/^^/ d'^g^'^'ix) 



T7;9'''Kx)9''^Hx)9 ' , ■ , - , 1 , ^ y 

^-[^2^ ^ ^ ^ '^dx'^dx'^dx'^^dx''^' 

+ off W( J ) + off ^(3-773-77) + 7r(5 ^ 



3 ' 3^ ^ ' ^ dx"-^ dx''2 ' J 6 dx'^^dx'^dx^^dx^^ ' 

(B.19) 

Computation of iy^e'^ 

The next job is the computation of the term ^F^: 
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^ dx'^dxi^dxi^ ^ ^^'^ dx'^dx^^ ^dx^^dx^^ 

d^ dAUr) r/^ 



dx^^dx^^dx^ da;'i dx'-^dxi 

^ dx' '^(ixJidx-?^ ^ ' dx'^dxi ^ ' ^ dx' ^ dx'' ' 

+ {9''-'{x) - g^^^^iy))B{x)^^^^ + Bix)A\x)^ + B{x)B{x)\ . (B.20) 



The trace is thus 



Tr(V2e^) = J d'^x{x\^Y^e^\x) 
= I d'^J^ (- J^Xx)5,u.(x)(^^^^) - \A\x)A\x)g,,{x) 

+r' + lBix)B(x)^ I . (B.21) 

Computation of l[X,Y]'^e^ 

We next compute the commutator [X, Y]"^. Prom now on, the computation becomes more compUcated 
than before, and we give only the trace: 



Tr(i[X,y]2e^) = J d''x{x\^[X,Y]^e''\x\ 



Computation of iy[X, Y\e^ 

rr(ly[x,y]e^) = y (i''x(x|iy[x,y]e^|x) 

= / d'^x — V <^ r -A'(x)( , / ) 9,1,0 (a;) + - A''{x)g^Ax){ , . M 
3 ^ dx'dxJ 3 ^ da;' 7 
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Computation of l[X,Y]Ye^ 

Tr{^[X,Y]Ye^) = J d'^x{x\^[X,Y]Ye^ \x) 



Summing up these terms altogether, we obtain the general form of the first-order coefficient ai{x) as 

1 . ..dgi^^^{x)..dg'^'^{x). 1,,,,,,,, 1,,,, , ^^dg^^'^ix). 

^ 1 k,k,, ^ . ^ I ,, dg'^'-{x) ^^ dg^-^-{x) ^ 

1 , ^, dg^-^Hx) ,, dg^-^^{x) ^ 1 , ^, dg^nx) ,, dg^^{x) ^ 
Especially, the Seeley-de-Witt coefficient of the Laplace Beltrami operator 



da:* 



(B.26) 



is important on the practical ground. This amounts to A^{x) — (y^^^j^f^^ ~ \9^'' {x) (^-^g'^'' (x)) gki{x) 
and B{x) — 0, in the general form IjB.f p . Plugging this into the general formula of Oi (jB.25|) . we obtain 
the following important result: 

a^{x)^^, (B.27) 
6 

where R{x) is the Ricci scalar of the spacetime metric. This result is used in the context of the superstring 
theory in deriving the Weyl anomaly and the corresponding Liouville action. 

Derivation of (|4.23|) 

With the above idea in mind, let us derive the result (|4.23l) for the simplest example of the calculation 
in our case. We expand the differential operator —tD^ given by H4.7|) around X = rdfi^d^. Namely, we 
express —tD^ as —tD^ = X + Y , and perform a perturbative expansion around X. The heat kernel 
for X is given by {x\e^\y) = J '^'^^ ^ exp | _ fcz^M^:!^ j _ Here, we identify the field a^'-\{x) with the 

vielbein, and we expand this around the flat background metric as a*^*)^(x) = (a^*)^)o + /i*p(a;), where 
(a^*V)o(«^''V)o = Then, D is more explicitly expressed as 



(B.28) 
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Now, we are interested in the mass term of the vector fields a/ii --/i2„-i (^)- the Campbell-Baker- 
Hausdorff formula ljB.9|l . only the terms Y and iF^ contribute to their mass terms. Now, we have a close 



look at these two contribution. 

For Y, only the elements of the rank-0 gamma matrices contribute to the heat kernel, because the 
trace tr for the 32 x 32 gamma matrices projects out the higher-rank gamma matrices. Therefore, the 
contribution to the mass term is clearly 

-la^^...^,ix)a^^-^^ix) - ^a,,...,,{x)a'^^ -^'^ix)^ , (B.29) 

where the constant 32 comes from the trace tr-l32x32- We next compute the effect of ^Y"^. This time, the 
higher-rank terms of Y contribute, because the product of the gamma matrices survives. The relevant 
terms are as follows: 

Y = T (^2m^(x)9^ - lr''^-2^^^^^^(^)5^ „ ^r-''- -^a^,,...,,{x)d^ + |r''--''^a^,,...,, (x)^^ 
2i 



where • • • denotes the omission of the irrelevant terms. This gives the contribution 

(xl^trY^e'^lx) = x32x (^-4a^{a:)a''{x) - ^a%,,,{x)a-'^^P%x) - ^a%,...M^^^^^^ 

-±a%,...,,ix)a'^''--P<^{x) - ^a'^,,...,,{x)a-'P^-'>^{x)^ d.d^e'^lx) 
'^\32r (a,{x)a'^ix) + ^a^,,,,,ix)a^^^-^-^ (x) + ^a,,...^,{x)a'^^-^Hx) 



where we have utilized the formula (|B.11|I for gi-^i^ (y) ~ riiii2 ■ We obtain the mass term H4.23|) in question 
by the sum of ||R29|) and l|R30|) . 



C Basic knowledge for the numerical simulation of the matrix 
model 

This appendix is devoted to introducing the basic knowledge of the Monte Carlo simulation of the matrix 
model. While the following ideas can be readily applied to the QCD, the programming for the matrix 
model is much simpler than for the quantum field theory. 



C.l Review of the Monte Carlo simulation 

The Monte Carlo simulation plays an extremely important role in the numerical simulation. In many 
cases, the analytical computation may be too complicated to handle. Let us think about the four- 
dimensional 10,000-site lattice gauge theory. This system has 40,000 sites. Therefore, even for the 
simplest Z2 gauge theory, we have to compute the sum of the 2"*°'"°° ~ 1.58 x lO^^^'*! configurations. 
The analytical computation of the partition function is almost impossible and this leads us to resort to 
the numerical treatment. 

To put it simply, the fundamental idea of the Monte Carlo simulation is to produce the equilibrium 
ensemble artificially. In the case of the matrix model, the state C is identified with the element of the 
N X N matrices; namely 

C-{(^^),,, •••,}. (C.l) 
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For the quantum field theory, C of course corresponds to the gauge field configuration Wj{x). Since 
the sequential integral is represented as / dA^ = X^c'''' partition function is rewritten as Z = 
Xp exp(— S'(C)). In the Monte Carlo simulation, we derive the configuration Ck that complies with the 
Boltzmann probability 

P^ol^Ck) = exp(-5(C,)). (C.2) 

To this end, we consider the following Markov process. Namely, the probability to transform from the 
configuration Ck^i to the new one Ck depends only on the state Ck~i and Ck- This probability never 
depends on the previous history Ci, C2, • ■ ■ , Ck-2- In this sense, this probability is denoted as 

P{Ck-i-^Ck)=^P{Ck-uCk), (C.3) 

since the probability is the function of only Ck~i and Ck- The following "detailed balance condition" is 
a vital constraint on the Boltzmann distribution: 

g-s(c)p(^^ C") = e^^(^')p(C', C). (C.4) 

For the transformation complying with the detailed balance condition, we have the following two crucial 
properties: 

1. An equilibrium sequence of the state transforms into another equilibrium state. 

2. A nonequilibrium sequence approaches to the equilibrium state. 

It is easy to prove the first statement. While it is trivial, we recall the following properties of the 
probabilities. Firstly, J2c ^i^^ ~ Tic ^(^' = 1- Secondly, when E' is obtained from E through 
the Markov process and P' is the probability for the ensemble E' , we have P'{C) = Tc' ^(^' C')P{C')- 
The first property is verified by taking the sum of the both hand sides of (jC.4ll with respect to C. Now, 
the previous state P\yQ\{C) complies with the Boltzmann distribution P\^q\{C) = Z~^e~^^'-^\ Then, we 
have 

P'iC) = J2 Pho\iC)PiC, C) = P^^^{C')P{C\ C) = Z-^e~^^^'\ (C.5) 
c c 

Thus, the newly emerging state also complies with the Boltzmann distribution. 

The second property is substantiated by defining the distance of the two ensembles E and E' as 



\E-E'\=Y,\P{C)~P'{C)\- (C.6) 

c 

Then, the distance from the Boltzmann distribution is given by 



i^^'-^boii = Ei^'(^)-^boi(c^)i = E 



5]P(C7,C')(P(C')-Pbol(^')) 



^ E E ^(^' C')\P{C') - Pbol(^')l - E - ^bol(^')l = \E- i?boll- (C-7) 

c c c 

This implies that the distance from the Boltzmann state is smaller as we iterate the Markov process 
obeying the detailed balance condition. 

There are two major algorithms that satisfy the detailed balance condition. One is the heat bath 
algorithm. The heat bath algorithm is used when we have a clear correspondence between the uniform 
random number and the configuration of the fields. In this case, the probability to transform to the 
configuration C" is given by 

P(C,C")oce-^(^'). (C.8) 

In this process, the probability P(C, C") depends only on the new configuration, and does not even depend 
on the previous state. It is trivial that this satisfies the detailed balance condition. In analyzing the matrix 
model, we resort to the heat bath algorithm. We postpone the explicit example to the subsequent sections 
for the matrix model simulation. 
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The other is the MetropoHs algorithm. Firstly, we choose the configuration C at random. Then, we 
compute the difference /S.S = S{C') — S{C). We adopt the new state at the probability min(l, e"'^'^). 
Namely, the probability is given by 

^ J I g-AS ^ g-S(c')eS(C) (if s{C) < S{C')). ^^'^^ 

We take a uniform random number r e [0, 1]. When r < e^"^^ we adopt the new configuration C", and 
otherwise we jettison the new configuration C" and keep the status quo. It is easy to verify that the 
Metropolis algorithm also satisfies the detailed balance condition, by noting that the probability for the 
inverse transformation is given by 

P(C' C) -^i '^' = e^(^'^e-^(^) (if S{C) > 3(0')), . . 

'^^ 1 1 (if s{c) < s{c')). 

The advantage of the Metropolis algorithm is that it can be applied to more various systems. When it is 
impossible to build the heat bath algorithm, the Metropolis algorithm is here to stay. A good example is 
the supersymmetric matrix model. When the model has fermions, we often use the hybrid Monte Carlo 
simulation. It is impossible to handle the Grassmann numbers numerically. Therefore, if we are to treat 
the model equipped with the fermion, we must integrate out the fcrmion analytically. Then, we evaluate 
the Pfaffian composed of the bosonic fields with the numerical method. 

C.2 Simulation of the quadratic matrix model 

We first start with the simplest example of the matrix model. While it is a very simple toy model, 
we scrutinize this case in full detail because this is the heart of the heat bath algorithm of the more 
complicated matrix models. We deal with the following quadratic matrix model: 

S=^Tr^\ (C.ll) 

Here, is an iV x hermitian matrix, and this model is invariant under the U (N) unitary transformation. 
Its path integral is given by 

Z = J d^Ve"^ = J d^Vexp(-yrr02). (C.12) 
Analytically, the propagator of this matrix model is evaluated as foUows^^: 

(0ji0fci) = ^ j d^"" 4'<t'i3<t'kieyiY>{~'^Tr(l)'^) ^ ^SiiSjk- (C.13) 

The proof of this propagator goes as follows. We note that, due to the hermiticity of (/>, the trace is 
written as 

N N ^ N ^ 

i,j=l i<i<j<N i=l 



^^The result iC.13i is applicable to the case <fi £ U{N). Here, we touch on the propagator for the SU{N) case. This can 
be derived from 1C.131 by subtracting the trace part. Namely, we replace (pa with <f>'^^ = (pa — X/i— i 't'"' propagator 
{4''ii4''ii) is computed as 

JV JV 
{<l>'rr<l>'u) = {<^>^^M " 2^ ^i'Prr'I'j j) + ^ ^'^n'Pkk) = ^ - ^) . 

Therefore, the Feynman rule is replaced for SU{N) case by 

{<f>i3<t>ki) = — (^SiiSjk - —S^jSki^ . 
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We separate (pij into the real and the imaginary part as 



(C.15) 



'2N 



Here, Xij and Yij are real c-number. The diagonal parts are real numbers, and they are rewritten using 
the real numbers W, as 



(C.16) 



Then, the quadratic term is written as 



N 



^i=l l<i<j<N 

The derivation of the propagator reduces to the simple Gaussian integral: 



(C.18) 



• {4'ii4'ii) = j^{WiWi) survives only for i = L 



N 

• For {ipijipki) {i ^ j), we note the following two results. Firstly, {(f)ij(f)ij) is shown to vanish as 



1 - 1 

cancelled (*) 



= - +2* X,,Y,,)) = ^ = 0. (C.19) 



Here, (*) does not contribute ab initio, since this is a linear term of each Xij and Yij. Secondly, we 
note that 



= ^{iX^,X,, + - 1 (C.20) 



survives (namely when i — l,j = k). 



This completes the proof of the Feynman rule ljC.1311 . 

Using the Feynman rule (|C.13|) . we can easily derive the following vacuum expectation values: 

(lTr0^) = l, (lTr^4^ = 2+i^, ((lTr</>^)2) = i + J_. (C.21) 

We next go into the numerical treatment of this matrix model, by means of the heat bath algorithm. 
The rewriting of the action ljC.17p plays a crucial role in the analysis. The advantage is that all Wi, Xij 
and Yij are decoupled. Therefore, in the heat bath algorithm, we can update each of these real numbers 
independently. Recall that, in the heat bath algorithm (jCSp . the probability for the new configuration 
depends only on the new configuration. Now, the partition function is rewritten as 



N ( \ ^ I 

WdW, n dX,,dY,e^v\-i^Y.^^-2 E 

1=1 l<i<j<N \ i=l l<i<i<JV 




(C.22) 

Thus, for example, the new configuration of Wi complies with the following probability distribution. 

P(M^.) = -^exp(-^), (C.23) 

where P{Wi) is normalized so that ^^'^ dWiP{Wi) = 1. Namely, the updating of the quantities 
Wi, Xij, Yij all reduces to the Gaussian distribution. The process of updating the matrix elements can 
be divided into the following three steps. 
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Generation of the uniform random number 

Firstly, we introduce the numerical method to generate the sequence of the pseudo-random number. 
A simple way is the "congruence method". We generate the sequence of the random numbers {zi} by 
the following recursive formula: 

Zi+i ^ az, + c (mod 2^^ - 1). (C.24) 

The initial term zi should be given by hand as a random seed. When we normalize this sequence as 
{ 231-1 }i obtain a uniform pseudo-random number [0,1]. Especially, it is known that we obtain a 
sensible pseudo-random sequence by taking a = 5^^ and c = 0. The drawback of the congruence method 
is the price of its simplicity. This random sequence has too short a cycle. However, this is not so 
problematic so long as we undergo the simulation of the matrix model, while this is a serious problem for 
the QCD simulation. 



Generation of the Gaussian distribution 

We next introduce the way to generate the Gaussian distribution from the uniform random numbers, 
generated in the previous step. We take two independent uniform random numbers x,y € [0, 1]. Then, 
we introduce the quantity r as 



r = V-a^ logx2. (C.25) 
The probability distribution of r is given by 

P{r)dr = P{x)—dr = ^ exp(-(-)2). (C.26) 
dr a a 

We next introduce the quantities X, Y as 

X = rcos(27r?/), y = r sin(27rj/). (C.27) 
We discern that this complies with the Gaussian distribution as 

P{r)drdy cx exp ^--^(X^ + Y'^)^ dXdY. (C.28) 

Especially when a — 1^ X and Y obey the normal Gaussian distribution. 



Updating of the matrix elements 

Now that we introduce the random number complying the Gaussian distribution, we update each of 
the matrix elements Wi, Xij and Yij by the Gaussian distribution. When we finish updating all of them 
once, this means that we have completed one sweep. Then, we reiterate the sweeps until the system is 
sufficiently thermalized. 



C.3 Simulation of the quartic matrix model 

We next investigate the (j)* matrix model. This model has a lot of lessons in treating the IIB matrix 
model and its extensions. The action is defined by 

S = N (^^Trcj)^ ~ ^Tr<j)^^ . (C.29) 

In order to perform the heat bath algorithm of this matrix model, we introduce the auxiliary field Q as 

S^^Tr<j)^ + ^TrQ^-aNTr{Q(j)^). (C.30) 

Here, a = and the auxiliary field Q is an x A'^ hermitian matrix. It is easy to verify that the action 
■3011 reduces to (IC.29|) by integrating out Q: 

S = Y^r{Q - a4?f + S. 
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The idea to update the quantities Q and (f> is similar to that of the quadratic case, and we proceed rather 
quickly. We rewrite the elements of Q as 



V -iV V ^jV 



(C.31) 



In this way, the action S is written as 



yi=i i<i<j<N 



s. 



In the following, VF^ {i — 1, ■ ■ ■ , N), Xij and (1 < « < i < ^) denote the real c- numbers that obey 
the normal Gaussian distribution. 

We next update the elements of the matrix (f). To this end, we extract the dependence of S on each 
component 0^ . Firstly, the dependence on the diagonal part (jiu is extracted as follows. In the following 
(up to (jC.47p l. we sometimes do not take a summation with respect to the duplicate indices. 



S = 



Na 



6„;)^ - aNQu{<i)iif' - aN(l)ii^{(f>ijQj^ + Qij(f)ji) 



• , where 



a 2 Ci 



(C.32) 



The • • • denotes the terms independent of (jju . We omit them because these do not concern the updating 
of 4)ii Therefore, the diagonal parts (pa are updated as 



JNci Ci 



(C.33) 



We next go on to the updating of 0i j , whose dependence is extracted as 



S 



f^hlt + . . . , where 




;Qkj + ^ Qjk<Pkj 
k^j k^i 



(C.34) 



This leads us to update the nondiagonal parts as 



(C.35) 



This completes the algorithm of updating the matrix elements. However, we must recall that this 
matrix model is unbounded below. The eigenvalues near the origin is stable only for the large- limit, 
and only metastable for the finite N. The function V{x) — N{^x^ — |a;^) has extrema at x = ±-^, 

at which V{±^) = ^. Therefore, the origin is barriered by the potential with the height For the 
large- limit, this barrier prevents the eigenvalues from dissipating outside the potential, and the system 
is kept finite. In [61], the following vacuum expectation value is computed analytically in the large- A'^ 
limit for < (7 < yj- 



lim ( — TrdP') 



a2(4-a2) 



where 



1 + VI- 12g 



(C.36) 
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Figure 21: The metastability of the origin for finite N. The height of the potential barrier is j-. 



Here, the parameter a gives the support of the eigenvalue distribution. Namely, the eigenvalues are 
distributed only in the region [—2a, +2a\. More accurately, the eigenvalue distribution is also derived as 

p(X)^i,/I^^(^AgX^-ga' + iy (C.37) 

where p{X) is normahzed as p{X)dX — 1. We delegate the detailed derivation to [61]. 

However, for the finite N, the height of the potential barrier is also finite, thus we see the divergence 
in the course of updating the eigenvalues. Namely, when the eigenvalues surpass the potential barrier, we 




Figure 22: The reguralization trick utilized in [64]. When any one of the eigenvalue surpasses the potential 
barrier, we take back the whole system several steps. 

see an avalanche of the eigenvalues rolling off and off the potential. This is what happens in the course 
of the divergence of the system. In order to evade the divergence, we need some regularization. Here, we 
utilize the same regularization trick as for the research of the Weingarten model [64] . Namely, when any 
one of the eigenvalues surpasses the potential barrier we take back the whole system several n steps. 

In explicitly writing our code, we store the history of several previous steps for the matrix elements 
of (j). When the eigenvalues happen to cordon the potential barrier, we recall the record of the history, 
and take back several steps (practically, around n = 3,4 steps). In this way, we evade the divergence for 
finite N. Here, we plot the numerical result based on the heat bath algorithm for the quantity {j^Tr4P') 
for N = 32. 

C.4 Simulation of the bosonic IIB matrix model with the Chern-Simons term 

We next introduce the algorithm for the IIB matrix model with the Chern-Simons term. In this section, 
we define the action as 

S = ——Tr[A^,Ai,]'^ — pYA^e^i...^2fc+i^^^Mi^M2 ' ' ' ^tiik+i- (C.38) 

Namely, we take the g in the action (|2.111|) to be p- = N . When we set A = 0, this is of course identical 
to the bosonic IIB matrix model, which has been scrutinized in [13]. There is not so big a difference from 
the bosonic IIB matrix model, because the Chern-Simons term is nothing but a linear term with respect 
to each A^. 
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Figure 23: The plot of {^Tr(\?) for the 0^ matrix model for TV = 32. 



We define this matrix model in the d = (2k + l)-dimensional Euclidean space, and the indices ^, z^, • • • 
run over 1, 2, • • • , d = (2fc + 1). The following argument is applicable to the bosonic IIB matrix model 
without the Chern-Simons term (namely when A = 0) for an arbitrary dimensions d > 3(the even d is 
also acceptable) For the bosonic IIB matrix model without the Chern-Simons term, the readers are to 
ignore the Chern-Simons term in the following argument. 

The action (|C.38() is also analyzed via the heat bath algorithm. Firstly, we note that the quartic 
commutator in IIC.381) is rewritten as 



^ Y.Tr[A,,A,r = -^ ^ Tr[A,,A^r^N ^ [Tr{AlAl) - TriA.A^AM] 



4 ^ L--M'--j 2 

p. 1^=1 l<p,<i^<d l<p,<i^<d 

N 



^ Y: TrGl + 2N Y: TriAlAl), 



(C.39) 



l<fj.<i'<d 



l<fj.<u<d 



where G^^ = {Af^,A^}, and these are hermitian matrices because these are anti-commutators of Ap 
This leads us to introduce the auxiliary fields Q^i/ as 



^ = ^ E (IttQI, - Tr{Q^,G^,) + 2Tr{AlAl) 



l<fj,<v<d 



(C.40) 



Here, Q^i, are hermitian matrices, and satisfy — Q^^. Q^i, is defined only for [i ^ v. Of course, the 
action IIC.40|l is equivalent to ljC.38|l after we integrate out Q^^y. 

N 



s = - - ^A^-)' + ^- 



l<^i<u<d 



^''The path integral of the two-dimensional bosonic IIB matrix model (without the Chern-Simons term) 

2 



-^Tr[AuA2? 



^,1^—1 



is easily shown to diverge [24] for any size of the matrices A''. When we diagonalize Ai as Ai = diag{Ai, A2, ■ • • A]v), the 
diagonal parts of A2 clearly do not contribute to the action. However, in the measure of the path integral 

the integral for the diagonal parts of A2 runs vacuously. This causes the divergence of the path integral. 
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Therefore, the matrices Q^i, are updated as 



(C.41) 



We next update the matrices Ap. To this end, we extract the dependence of the action ljC.40p on each 
component of {Ap)ij. To this end, we rewrite the action HC.40p as 



-N J2 TrQpMpA. + A,Ap) - nJ^ TrQpp{ApAp + ApA^,) 

p<v fj,<p 



2k + 



TrA 



' ■^V2k-^P 



Vi+P 



So = 



-TV J2 TrQ^piA^Ap + ApA^) + 2iV ^ Tr(44) - A7V ^ e.,...,,TrA, 
-NTr(TpAp) + 2NTr{SpAl) + ■■■, where 

^p; ^p — ^(^A'QpP + Qpp^pi) + ^ ei/i - i^2fc^'^^i'i ■ ' ' ^V2k ■ 

p¥=p Mt^p '^i¥'P 



(C.42) 
(C.43) 



The difference between the pure bosonic IIB matrix model and the model with the Chern-Simons term 
only comes in the definition of Tp. Otherwise, both cases go totally in the same way. 

Using the rewriting l|C.42|l . we first update the diagonal part {Ap)ii. (|C.42|I is further rewritten as 

S = -N{TpUAp)u + 2N{SpUAp)l + 2iV^[(5,),,(A,),,(A,),, + {Sp),,{Ap),Mph] + ' ' ' 



2N{SpU{Ap)l ~ ANh, + • • • = 2N{Sp\^ {Ap) 



{Sph 



2N- 



iSp)^ 



where (C.44) 



2^[('5p)ji(^p)y + iSp)ij{Ap)ji] 



{Tp)ii 

Therefore, the diagonal parts {Ap)ii are updated as 



+ 



^ANiSp)u (Sph 

We next update the nondiagonal parts {Ap)ij, whose dependence is now extracted as 



(C.45) 



-7V[(r,),,(Ap)y + (rp)y(Ap),,] + 2N[iSp)u + (^p),,]|(^p)yf 



+ 2N 



''^iSp)jk{Ap)ki{Ap)ij + 'y^^{Sp)ik{Ap)kj{Ap)j 



+ '^^{Sp)kj{Ap)ji{Ap)ik + 'y^^{Sp)ki{Ap)ij{Ap)jk 

k^j k^i 

2Nc,j\{Ap),j\'' - 2N{Ap),jhj, - 2N{Ap)j,h,j + • • 



= 2Nc,, 



(^p)ij - 



-2N^ 



{Sp)ii + {Sp)jj, hij — ^{Tp)ij 



■ ■ , where 
^{Sp)tk{Ap)kj 

.k^i 



(C.46) 



J2iSp)k,iA, 



p }ik 



Therefore, the nondiagonal part {Ap)ij is updated as 



This completes the algorithm for the heat bath algorithm for the bosonic IIB matrix model with the 
Chern-Simons term. In order to corroborate the legitimacy of the code, it is useful to exploit the following 
relation coming from the Schwinger-Dyson equation 

= / '^'^ E E aMiTr{e^,)e-% (GAS) 

J a=l M=l 

Here, i° is the basis of the SU [N] Lie algebra and satisfies the following identities: 

N 



Trim'') = J2 = Sri^jk - -^kAi. (C.49) 



0=1 



The matrices are expanded in terms of t"" as — 'Yl!l=i ^ ^u^"- This leads to the following identities 



Tr{eA)Tr{eB) = ^ Aj,Bik{e),j{e)ki = A,,Bik{5aS,k - jz5^,5ki} 



N 

a— I a— I 

= Tr{AB) - ^TrATrB = TrAB. (C.50) 

At the last equality, we utilize the tracelessness of the SU{N) Lie algebra (namely, TrA = TrB = 0). 
Recalling the fact that 

dS 

— = -NTri^A,, [A„A,]]) - XdNe,,,...,,_,Trit'^A,,A,., ■ • • A,,_J, 
we rewrite the Schwinger-Dyson equation (jC.48|l as 

= d'^AY^ [Tr(tH'')de-^ + NTr{t''Af,)Tr{t''[A^, [A^, A^]])e-^ + 

a=l 

+XdNe,,,,...,,_,Tr{t''A^)Tr{eA,,A,, ■ • • A,,_Je-^] 
d'^A [d{N^ - 1) + NTr {[A^, A,f + \e^,...^^A^,A^, ■ • • A^J] e"^ 

(C.51) 

This implies that the vacuum expectation value of the following quantity is analytically computed as 

- ^(rr[A^, - h.l^e^^,„^^TrA^^A^, ■ ■ ■ A^,) = d{l - (C.52) 

C.5 Jackknife (binning) method 

Finally, we introduce the jackknife method to estimate the errorbar of the quantities obtained by the 
numerical simulation. Now, we have a sample of independent measurements of a primary quantity A. 
The measured samples are 

Ai,^2,---,^jv. (C.53) 

The sample average is simply defined as 

N 

N 



1 

A=-YAs. (C.54) 



Next, we define the jackknife average. We exclude one of the sets {An\ in taking the average: 



A 



{.J)s 



N 



^—YAr. (C.55) 
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The variance of the jackknife estimators can be obtained as 

N-^ ^ 

^l)A = ^T.(^iJ)s-Ar. (C.56) 



Note that the average ^(j)s is trivially equivalent to the average Ag. We can verify (7(^j)a = (^A 
following way. 



Ar_1^/_ 2-- 1 



N-1 



(TV 

Now, cr^ is an unbiased estimator of the variance: 



E 



Next, we explain the jackknife binning. This method is utilized when Ag's have correlations. Espe- 
cially, the configurations generated by the Markov process are not statistically independent. This poses 
the problem of the "autocorrelation". In coping with this problem, we use the jackknife binning. We 
consider the following bins: 

(C.59) 

one with /' = elements »^v- biu one' bin 



M bins 

Suppose one bin contains n elements, and there are M bins. Then, we take the average for each bin as 
follows: 

a^^(Bin), 

where b runs over 1, 2, • • • , M. We take the variance for these ^(j)6: 

M- 1 ^ 



6=1 



Note that = A. This variance depends on the number of the elements in one bin n, and this 

variance increases with n and we finally see the plateau. 

To illustrate this process, let us take one simple example. We have the data 

Ai = l, A2 = 2, ^20 = 20, (C.62) 

and we construct the bin as 

{A,,---, A4, {A, • • • , As}, • • • , {A,r, • • • , A^o}. (C.63) 
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In this case, N = 20, n = 4 and M = 5. Namely, we define the average for each bin as 

(5 + 6 + 7 + 8) + (9 + 10 + 11 + 12) + (13 + 14 + 15 + 16) + (17 + 18 + 19 + 20) 



^(J)b=i 



^(J)b=2 



^(J)&=3 



^(J)b=4 



20-4 

(1 + 2 + 3 + 4) + (9 + 10 + 11 + 12) + (13 + 14 + 15 + 16) + (17 + 18 + 19 + 20) 

20-4 

(1 + 2 + 3 + 4) + (5 + 6 + 7 + 8) + (13 + 14 + 15 + 16) + (17 + 18 + 19 + 20) 



= 12.5, 
= 11.5, 



A 



(J)6=5 



20-4 

(1 + 2 + 3 + 4) + (5 + 6 + 7 + 8) + (9 + 10 + 11 + 12) + (17 + 18 + 19 + 20) 

20-4 

(1 + 2 + 3 + 4) + (5 + 6 + 7 + 8) + (9 + 10 + 11 + 12) + (13 + 14 + 15 + 16) 

20-4 



= 10.5, 
= 9.5, 
= 8.5. 



The average is of course ^(,/)b — 10.5. The variance is computed as 



plateau 




Figure 24: The plateau emerges as we increase the bin size. 



7L,„) = ^—^[(12.5 - 10.5)^ + (11.5 - 10.5)^ + (10.5 - 10.5)^ + (9.5 - 10.5)^ + (8.5 - 10.5)^] 





4 

- X 10: 

5 



(C.64) 



The process in taking the plateau is too delicate to delegate to a computer, and we usually judge the 
plateau by hand. 
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